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A^' ^^OTATIO N 


The present issue of the scientific works consists of two parts. The 
results of the flow around wings, the determination of the optimal forms, and 
the interaction of the v;ake v;ith the cocurrent flow at supersonic and hypersonic 
speeds of the free-stream flow are given in the fii’st part. 

The second part comprises methods of numerical and analytical calc’lLation of 
one-dimensional unsteady and two-^iimensional steady motions of fliel-gas mi:ctures 
;d.th exothermic reactions. 
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fore;ord 


The present collected papers are dedicated to the aerodynamics of the super- 
and hiypersonic ilov^s and theories of th-- supersonic motion of a gas with v/aves of 
detonation and combustion. In ub.ese trends, investigations are performed and 
seminares are made in the Institute of Mechanics during many years under the 
leadership of the associate member of tnc- Acadeiiiy of Sciences of the USSR, 

G»G. Chernov. During the last time, a number of new results were obtained, 
compris:Lng the contents of the given issue of the scieiitific vrorks. In the 
first part, problems of the flow around triangular and V-shaped wings with 
gas flow of tilgh supersonic velocity ai-e investigated- (A.L. Gonor, N.A. Ostapenko 
and V.I. Lapygin) « Further, in the paper of L.V. Gogish and G.Yu» Stepanov, an 
integral method of calculation of the interaction of turbulent wake 'with a flow 
in a canal or a jet is suggested. Problems of the optimum forms of a body at 
hypersonic velocities, in the frame of the approximation law's of resistance, are 
resolved in the v;orks of A.L. Gonor and V.I. Lapygin. The results of an 
experimental investigation of supersonic flow around a right-angled wing, in a 
wide range of elongation and angles of attack, are given in the paper of 
M.P.- Falunin. G.S. Ulyanov studied the effect of the penetrability of a plate 
on the power action of the supersonic floii/. 

The second part of the issue, comprising the methods of calculation of the 
motion of fuel mixtures, is started with the paper of S.A. Medvedev. In this 
paper, the problem of the decoy of the burst is resolved in non-s elf-similar 
arrangement, tal-cing into account the ignition delay and the final reaction rate. 
V.P. Korobeinikov, V.A. Levin and V.V. Markov use similar models of fuel medium 
to solve the probic.s of point explosion v/ith plane, cylindric;;il and spherical 
waves. In the paper of L.I, Zak aJid V.A. Levlii, the gas motion, caused by a 
piston, is investigated assuming that the heat effect of the reaction is small 
and that a linearization, relative to the adiabatic flov/ in front of the piston, 
can be produced. The flow ju-ound a body by a stationary supersonic flow of liiel 
mixture is investigated, in the last three works of the collection. In the first 
one of them, S.M. Gilinskii irivesuigates the growth of the initial disturbances 
£ind the possible nonstationary conditions of combustion at the hupersonic flow 
around the bow of a blvait body. The second paper of S.M. Gilinskii and 
M.L. Khaikin is dedicated to the use of the boundary layer method for the 
composition of an anslj’tical solution of the supersonic flow ar-ound a wedge 
and a cone by a fuel mixture. The third vjork of S.M. Gilinskii "Linearized 
supersonic nonequilibrium flov/ of a fue?_ miicture of gas vs near a wedge" 
comprises an analytical solution of flov/ around a v/edge in the case of a small 
heat- effect of reaction,- obtained, as in the work of L.I. Zak and V.A. Levin, 
on the basis of using linearization relative to the corresponding adiabatic 
flow . 
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FART I 

3 S'JPRR- AND HYPERSONIC GAS FLOW 

■ CALCUUTION OF THE MTROP^’' LAYER ON THE SURFACE 
OF A TRlAI^aUlA.R WING 

By 

A.L. Gonor, N<.A. OsfcapenJco 


On solving the problem of hypersonic flow around a triangular wing [l], a 
singularity in the distribution of the surface speed is involved. For this reason 
the obtained solution must be corrected in the region of the location of these 
singularities, adjoining to the plana of symmetry. A similar situation already 
arised in the general theory of the supersonic conical Tiow, when ^uch singularities 
appeared in the solution on the surface of the body and it was necessary to make a 
specific solution in a region, called th.- entropy layer [2, 3» 6-8], The calculation 
of the entropy layer on a vring ha 3 _mjny _things_ in common \-d.th the mentioned 
inve stigations ; however there are important differences as a result 
of taking the specific features of ^iurfaco solution into account. The 
problem of the flow around a plane triangular wing is analyzed below 
In detail, although this method cem bo triinsposed also for the general 
case . 




1« The formulae of a conical, flow in the coordinate system cf'2'^0_jq? 
(Fig« i) after ^transition to the new variable (9-y satisfying the 

formula v '-P ja ^0 according to reference [l] vdll have the following 

f oimi ; 


cciQ 


coiQ 3f 

JL f- 

U-l 


u^*V~*W 


^■C 


a<f i ■! 


in(fwe^)*2-’^:. ^0,9^0, 


dif 


w9^ = vC&z B . 


O-V 


w aiv ■in ^ <3P fi,2^ 

cosy ..1^7 ? jjcot;- - ■yip Q ■f, 0 ip J 


Here u, v , W — are the corresponding velocity projections on the axes -'jr^PpCf 
relative to the. velocity of the free— streara flow Tj^. ; Pj /■’ ? V are the 
pressure, density and ratio of specific heats. The density and pressure are 
relative to the Quantities and Jl^ n 1 * 030 . ctively. The coordinate 

const determjjaes the stream surface. The projection of the Eulerian equation 
on the direction of the axis cf (le2) is ^n’itten for further direct use. 

In the system (l»l) and..in-tha_bQvmdary_conditions on the shock wave 
after the substitution of 6"-cd* = V=£v the terms, having 

the order of £ are neglected { £ is tne characteristic ratio of the 
densities in front of the shock wave and beiiind it, 1 '> 0). From the 
solution in reference [l], the forrauia for the pressure vdll have the following 
fonn: 


P . sZ/fi/ *6P^ /P, f, 9 ) , 


• ij) 


where : 


i'-n Zcl(8^cc%<p- iinf)- %in~oC . 





Here the bar of 0'^ is caiitedj c>< is the angle of attack. The value of P^, 
and the other parameters of the flow are expressed as r,'' by the solution of 
the integral-differential equation deoemiining the form 'of the shock wave. 

The surface of the wing in the solution of reference [l] is the stream 
surface is the half-angle of the triangular ’.-ring at the top). 

However, in the small .interv'al 0^ cp ^ cp^ the stream surfaces cross the 
surface of the vring (a weak flcv'^ing of makes place). Consequently, it is 
necessary to make a correction fcj- the solution in cermain surroundings of 
the origin of the coordinate. 

Let us introduce the variables A and “ZT instead of u and w; 


tf cos T 


V = iZl s i/7 T 


Formula (l«2) by using (l.l) vrith the nevr variables ;^11 have the form; 


, / dz 

X +: 

J Ir/ 


' COi Bd^ 


flL 

pL'CotQ L 9 'f 


OP 1 

r J 


0^ o ^ 




It is easy to ascertain that, according to (l*3) and (l.4)» the pressure 
gradient will be taken into account only with terms of . Therefore with 

an accuracy cf values of g ^ , \ie will have 


^ i )^0 . 

o * «? * 


(1.5) 


Accoi'cJiiigljj the equality of the formula in the brackets to zero, by using the 
condition on the shock wave, gives the follov.'ing solution wiiich is coincident 
'..'ith [l] ’■ 

( 1 . 6 ) 


where ^ ( cp ^ is determined from the condition on the shock v/ave at 
c^'= cf (Fig. 2) . 
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Fig, 2 


In the surroundings of the shock wave, the second possible solution 
contradicts with the boundary condition. Therefore the surface solution can be 
only formula (l*6), at the same time the transverse component of the velocity W 
has a singularity on the surface of the wing. This shovjs that in the region of 
the flow v/here , , the small corrections, neglected in the derivation of 

formula (1.5) i become comparable \d.th the remaining terms. In particular, the 
pressure gradient appears essential along the wing span. Let us show this, 
introducing for convenience the following designation: 




(1«7) 


Then, putting z [^below, ,the upper index (k) 

will denote the coefficient, standing at , where corresponds to (1.6), 

for we obtain the following relation: ! 
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From equation (l.8), it is shovm that has a singularity at the points 

where becomes zero. Hereby^ the essential term is determined by the 
pressure gradient 

At 0 t the singularity certainly exists} at b( 2) « 0 (plane of symmetry) 

it may be not existing and therefore a special investigation is needed here. Thus, 
in the- solution of [ij (below we shall call it the surface solution), on approaching 
to a certain critical surface a special term appears in the 

formula of , namely? 

This inequality in the solution indicates the presence of a certain term. 

Although this term is negligibly small in -&he surface flow and evenly small in 
all the region, a quantity of a "principal" order appears in the region, vjhere 
nonuniform convergence exists. From the foregoing, it is shovm that such terra 
is a function of B . 

2. The equation (l»5) has two solutions j 


~ =y Tr:0 . - • 

d ■ ' , • ■ 

The first solution has a singularity near 0 • Therefore in the 

investigation of the , ununiform corjvergence of the solution, the analysis of 
the expansion, obtained from the second solution of [?]» is completely reasonable. 

Let us represent T t: (where i- denotes the inside solution) by the form 


t. - 

C 


.2 (l) 

f T. 


( 2 . 1 ) 


Substituting formula (2,l) in equation (l»4)f v;e obtain after a simple 
calculation: 








( 2 . 2 ) 


c 



•• 8 •• 


By means of (l»8), we find that 








TiJ/. ( i) 


~lj- _ ''t.y i 7 

C‘>V«3v'> ^ XT / 


(2.3) 


hence it is clear that| pJ^l must be determined for the 

determination of '' 

Let us assume that these quantities can be determined, and let us consider 
the question of coupling of the inside and surface solutions. To prove that the 
two solutions are coupled, it is necessary to ascertain that they have a common 
region and are asymptotically equivalent in their common region. On the basis of 
the second of these two conditions, it is possible to conclude that the inside 
and surface solutions are not coupled, because these solutions intersect and the 
slope of the curves has a discontinuity at the point of intersection. Thus, 
the first terms of the expansion can not be asynptotically equivalent in any 
common region. 

This might have been expected for tv;o reasons* the scale, of the independent 
variable is the same in every region and the arbitrary constant (or, moreover, 
the arbitrary function of cp ), by which the coupling of solutions might have 
been realized, is absent. 

The problem of coupling can be investigated by several methods [4»7]. 

The method, used below, is closely related to the behavior of the surface 
solution. On. analyz'-lng the surface solution, it was clear that the tertn of 

second order of smallness has a logarithmic singularity at cp C^’') • 

The terms of the next order have still more essential singularity. It is possible 
to show that the ratio of the succassivo terms is an indefinite quantity in the 
surroundings of the singular point* It is possible to eliminate this divergence 
of expansion by the modification method of expansion, as shown in the method of 
PLG [ 5 ], in which the dependent and independent variables are expressed by an 
additional variable and are decomposed in series of £ « This gives additional 
freedom, v;hich can be used to a certain extent for the control of the behavior 
of ec.pansion near the singularity^ 

If we consider equation (I.4) for small values of ‘Y and introduce the 
following approximations: 


U) 


/ham ■ C'p'] 


Cat B - •* 2 ' + '•• , 
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This equation •, indicates the possibility of using the method of PLG: 
Let us rewTite the equation (1*4) in the form: 


dx _ *c)cozB 






Let us introduce the auxiliary function of -P(z) and the variable -z end change the 
previous equations to the equivalent system of the equations: 

■ ■ 

'dr Coie ' 


The function of^(z) is selecjted_Jja_ such a manner that the expansion forT^ 
and Lf begins with the terms and respectively. Thust 

The principal set of equations will take the following form : ' 




^ d2 coid 


(2.4) 


All the variables can be represented in the form of the following series: 


C-Z £ C . , 

n^O 


r = T 


••• , 


(/>» Z + £ i f 
^ ml 1 • 



The functions of c''"^ ' and cos B will be expressed by Z: 

) (• 

^ I^J = C (Zj^-i — 


£92 


.2 0 




Cas9=i-£ • Y + ... 


dz 


(t) 

’f + . . . 



At the same time ^ (fj-0 and ' * because the principal 
terra of the pressure is constant along the vring. Substituting the expansions 
in equations (2./*,); v:e shall have, after integrations 












d2 


ijt SUl[z-i-e[fCf')]cOs[z*<^fff'}] 


C2-^^ 


It is clear from formulae (2.5) and (l«8) that it is enough to know the solution 
of [1] for the determination of -TrC.^^ and in which we are interested. 

According to (2.5) the formulae for and of vail have the following form; 






• V»=2-f^/^ — <p')]Cos[z J-c/f /V;)7 


(2.6) 


The obtained relations enable to extend the uniform convergence of the 
initial solution to a small region in the surroundings of the singularity. 

If we denote ' d- PV/Q<^ then for the new variable Z 
we have: — — . 

■ ' ^ ; V,'; ' , 

The equation continuity (1.1-5) can be written in the form 

+2t/ci^r°»0~ 
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from v;hich v.'e obtain 

here (’) denotes that the value of the variable is taken on the shock wave or 


ow9 , , 


where I for example, 


ixy'-wf 

/2rve^ ; 

After integration, it is possible to find the corrected formula for the line 
of flow 


^ • 
f ' 5 If 

R=u(B'cn f'-el s« * w •'’J • 


In this manner, the convergent solution in the whole outside region, 
including the critical surface and the small region around it, is performed 
uniformly. 

3» To obtain the expansion in the ',:hole region it is necessary to 
find the total inside solution. 

The principal terra for ^ in the inside region has the fom 


\T. =£ ; 0-0 

Consequently, w^ O(s*) . 

Using these relations in the exact equation (l.l) for the variables of the 
inside region, vje obtain: 
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li. 

d tp 


’o'; 



0(i‘) , 


. C3-?0 

2 

It is clear from equation (3»2-2) that the pressure with an accuracy reacVn-n^", t 
across the entropy layer is constant. Consequently, the value of C 2 ,^"^can be 
determined by the press\ire on the critical surface, i.e. by the value of the 
pressure in the surface solution, because and only the term Bf'-’ 

containing the pressure gradient, remains. These considerations can be applied 
also to the quantity C If we investigate the asymptote from of (2.5), 

then on the critical surface, we have; 




So, the surface solution, performed by the method of PLG, is asymptotically 
coupled with the inside solution. 

From equation (3»2--l), vie obtain 




Fxjrther we find that 



ja 


dr 



(3.3) 


where d^ is determined from the condition of coupling vdth the surface solution 
in a certain section cpy • 


• fc/'exp[-2/ ct^T°c/zJ 

Substituting the formula for dr in (2.9), we °ht_a^_ 

/ ‘ fr ^ 'f' 

Further, taking into account (3«l) end (l»3)» will have 


(3.i>) 


( 2 ) 1 


l^dn2cZiUii,'>(e**e*^)-P^^J • 


(3.5) 
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The calculation! performed by the formula (3»5)i showed that the value of~^ 
is negative and the first component in the square bracket is much larger than 
the value of which, in turn, is proportional to the variable c-p . An 

analogous derivation can. be obtained if the method of disturbance given in 
reference [l] is, used for the determination of the form of the shock v;ave« 

Then, for the small values of cp , the equation of the shock wave and pressure 
gradient allow simple analogous representations, from v;hich it follows that 
has a higher order of smallness in comparison tro the first component in the 
square bracket of (3 »5) • 

Now, taking into account (3«4) and (3»5), it is possible to draw a conclusion 
that 9c > 0 and Qi ^ 0 "hen cp^O . 

So, the performed solution, consisting of two asymptotic conjugated derivations, 
is valid in the whole region of the bad flow and totally satisfying the Fery's 
scheme. 
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The graphs of the transverse component of the velocity W at' «* 45®, 

33 n 1,4'i «= 20® and 30° sre show.' in Fig. 3» They ai'e calculated by 

formulae (2.6) and (3.5) (solid line) «nd by the corresponding formula of 
reference [l] (danhed line). The comparison shows that, without taking account 
of the entropy layer near the pla-he of symmetry, the velocity distribution on the 
wing at big angles of attack is appreciably decreased. 
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ABOUT THE MODE OF FLOW AROUKD V-SHAPED '.VITH 

SUPERSONIC LS\DING EDGES 

By 

V.I, Lapygin 

f 


It was shown by experimental investigations [1,2] that on flowing around 
V-shaped vrijigs, a complex system of shock waves is formed in the field of the 
flow. In the present work, a system of shock waves foniied on flowing around a 
V-shaped wing by an ideal gas, is investigated by means of performing numerical 
calculations on the £VM BESM-6, using the method discussed in reference [3]» 

V/e shall consider -the flow around a V-shaped wing at an angle of attack 
(Fig, l). In the conical coordinates y 

the system of equations of gas dynamics (for an ideal gas) in the divergent form 
will be as follows: 









' 


• 





pw 



j=u'y 

a 

JD«U^ 

- 




2 

punf 


pv-vr 


/7+_p^ 

E 

w 


'e% 


. ^ 




xp 




) , X-Cp/c^ 


CD 


here u, v, w are the components of the velooity vector on the axes x,y,z; p andj’ 
are the pressure and density. 

Using the condition that the s'upersonic flow around a V-shaped wing is a conical 
flow, the method of determining %, with the use of the difference scheme proposed in 
reference [5] is used for the numerical integration of system (l). The calculations 
were performed for different values of and , The analysis of the performed 



calculations allowed to reveal several modes of flow around as a function of the 
expansion angle of the V-shaped wing, with qpialitatively different systems of shock 
waves from one shock wave attached ao the leading edge of the wing to a system of 
four waves and more. In particular, the performed calculations shov/ed that there 
are flow modes with a reflection of a strong shock from the symmetry plane of the 
wing. In the plane case, as it is known, a strong reflected shock is not realised. 
The marks used in the text are clear from the consideration of Fig. 1. 



around with a mode close to the outgoing shock wave. The foms of the shock 
waves are given in Fig. 2, and the distribut:ion_of the pressure coefficient 
Cp„ - (_p_ on the wall ■ is given in Fig. 3. 

In the graphs, ‘-r is the distance . from the plane of symmetry to the given point 
on the v;all of the wing, R is the distance from the plane of symmetry to the 
leading edge. 



3 
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C 




The values given belov; in the graphs and tables are dimensionless. These values 
are namely: the component of the velocity vector, relative to the modulus of the 
velocity vector of the undisturbed flow; the pressure, relative to the doubled dynamic 
pressure of the undisturbed flow; the density, relative to the density of the 
undisturbed flow . 

On treating the results of the numerical calculations, the position of the shock 
wave is determined as the region of the sudden variation of the parameters of the flow. 
The results of the calculation of the configuration of the shock waves are represented 
in Fig. 2. ' 

The relations of Cp are given in Fig. 3j and the analysis of the field 

of flow, performed for different angles of 'i at o 3 » 95 » show that with 
decreasing , a transition from a flow around with a convex shock wave, through a 

C mode of a plane shock [4]f to a flow around with a concave shock wave occurs. In the 
case of the existence of the concave shock vrave, an increase of the pressure takes 
place towards the center of the wing. At o 120®, this increase can already be 
identified by the presence of an internal shock v/ave; i.e., the flov; around a 
V-shaped wing occurs with a formation of Mach configuration of the shock waves. 

The question about the foxnnation of rn internal shock wave at once at the transition 
from a mode of a flow around with a plane shock to a flow around with a concave 
shock wave, remains open. As an example, the field of the isobar for X = 120° 
and ox « 150 is given in Fig. 4« The character of the isobar indicates the flow 
around with a Mach configuration of the shock waves. The position of the plane 
shock, formed on the leading edge, is represented by a dashed line in Fig. 4» 

On decreasing the angle X to 40°, the character of the isobars (Fig. 5a) 
indicates a regular reflection of the shock, incident from the leading edge, from 
the plane of symmetry of the wing. It is known that the reflection is theoretically 
possible with both a strong and weak reflected shock, but in the plane case, the 
reflection v?ith a strong shock is not realized. In the present work, to determine 



r 
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the type of the reflected shocks, snd also the character of the reflection (regular 
or Mach), the results, obtained by numerical calculation, are compared with the 
calculation of the syste;:; of the shock waves performed by the fomirula of the oblique 
shock. Such analysis showed that for n /^O® and o 15°, a regular reflection of 
the plane shock, incident from the leading edge, from the plane of symmetry of the 
wing, occurs with a strong reflected shock. In this case, the values of the density 
and pressure, obtained in the numerical calciilation and by the formula of the oblique 
shock, differ by not more than The system of the shock waves in that case is 

given in Fig. 5 a by a dashed line. . " 
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With decreasing the angle T , the system of the shock waves, formed on 
flowing around the wing, is changed, and at '(fa 20°, the field of the isobars 
(Fig. 5 h) indicates the presence of a shock wave, reflected from the plane of 
the wing. The comparison ;d.th the calculation perfonned by the formula of the 
oblique shock, showed that a system of shocks is realized. This system is represented 
in Fig. 5 b by the dashed lines, and namely: the incident shock wave, from the 
]^ading ed^e, is rAflected., from., the, plane of symmetry qf^^^ by the_regular 

way' with a weak wave which. -is regularly reflebtetrfrom the plane of the wing with 
a weak reflected wave. The latter is reflected from the plane of symmetry with the 
formation of a Mach configuration. As for ■ AO'^ the difference of the values of 
the veclocity, density and pressure obtained by the nrunerical- means and by the\formula 
of the oblique shock did not exceed 3^. 

With further decrease of the angle Y ( X = 10°), the number of reflections 
from the plane of the wing and the plane of symmetry is increased. In this case, 
the last reflection from the plane of symmetry is regular with a strong reflected 
shock. All the rest reflections, observed in that case, are regular with Vfeak 
reflected shocks. The system of the shock waves is given in Fig. 5 c. The 
difference of the values of velocity, density, pressure, obtained by the numerical 
means and by the formula of the oblique shock, as for the previous variant, did 
not exceed 2^. It is necessary to mention that in the numerical calculation, the 
region of the flow, designated in Fig. 5 b, c : by the number 2, .is' poorly studied 
and is not practically defined in view of tlie small number of the points of the 
calculated net, falling in it. This case is caused by the small dimensions of the 
region 2, in comparison with the calculated region of the flow. As an example, 
the results of the numerical calculation are given in Table 1, and the results, 
obtained by the formula of the oblique shock are given in Table 2 for V « 10°. 

Table 1 


HoMep c6- 
aacTH Q 

u 

V- 

rtr 

P 


1 

O.SS46 


o 

0,2,23 

3,004 

; 2 






3 

0.Q500 

-0,2143.. 

o 

0,0766 

1,^144 

4 

0,9C23^ 


-0,01994 

0,0678 

1,320 

> 5 

0,9574 

-0,2394 

O 

0,0590 

1,211 

6 

0,9617 

-0,2494 

-0,02183 

0,0522 

1 ,098 


l) Region number 
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Table 2 


HoMep o6- 
aacTH (D . 

r " ■" 

a 


TJf- 

P 

P 

' 1 

0,3995 

-0,00579 

O 

0,222 

3,020 

2 

0,9472 

O 

9 

-0,01759 

0,0865 

1 ,577 

3 

0,9507 

-0,2169 

O 

0,0764 

1,^143 

4 

. . 

O.Sotn 

-0,2283 

-0,01998 

0,0678' 

! ,325 

5 

0,9576 

-0,2403 

O 

0,0599 

1,213 

6 

0,9617 

-0,2494 

-0,02184 

0,0522 

1,099 


1) Region niamber. 


As evident from the comparison of the tables, the coincidence of the results 
of the performed calculations is good, which allows to establish uniquely the 
system of the shock waves formed on flowing a:*ound a wing. The norabers of the 
regions mentioned in the tables, correspond to the numbers mentioned in Fig. 5 c. 
The values of p P jf ^or the region 2 are not given. The values of 

components of the velocity v for the region 1 (Fig. 5 c) in Table 1 are not given, 
because in this region the velocity V smoothly decreases to v = 0 on the X~axis. 
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The relations • P::i P/'^^oe> on 'the axis of symmetry (l) and on the wall of the 
wing (2), obtained in the numerical calculation, for Y, « 10°, = 15°, 

a 29° 30* , /to “ 3*95 ai'e given in Fig. 6« In this figure, the ratio C^/A) 
for the axis of symmetry is the ratio of the distance from the intersection point 
of the planes of the wing to the flovdng point on the axis of symmetry to the value 
of the projection of the leading edge on the plane of symmetry. The dashed lines 
represent the solution obtained by the formula of the oblique shock for the system 
of the shock waves shown in Fig, 5c. The character of the curves confirms the 
existence of the .system of the shock v;aves mentioned above . 

The relation, given in Fig. 6, shows an enough big spreading of the fronts of 
the shock waves at the points of regular reflection. This is related, on one 
hand to the presence of two fronts, the spreading of which is increased approximately 
tv;ice in comparison with a single front, and, on the other hand, to the small 
intensity of the shock waves. 

Calculations of the flow around a V-shaped wing at 140°, 15®, at 

different values of Mach number , were also performed. The results of these 
calculations' are not given here, Ihe analysis of the calculations showed that 
with increasing the Mach nuirber M^ , a transition occurs from a flow around the 
wing with a convex shock wave to a flow around with a concave shock wave and Mach 
configuration, i.e,, the transition is qualitatively as on flowing -around at 

const with decreasing and with values of ■‘X that .are riot too small. 

In this case, at "f “ ^ regular reflection is not observed for any- Mach 

number M^ up to M^^a ^ and a system of shock waves of Mach type is realized. 

In this manner, the perfoimjed analysis shows that at sufficient big angles of"^ 

(in the considered case for ^ « 29° 30* the angle Y .^80-100°), a complex 
system of a large number of shock waves is not realized, and namely these angles 
are of greatest practical importance on constructing aircrafts with high aero- 
dynamic quality. 

The results obtained in the present work agree qualitatively vjell with the 
experimental investigations given in reference [l]. 

The author is grateful to A.L, Conor for ' suggesting the problem and for hds 
attention to the v;ork. 
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QUASI-ONE-DIMENSIONAL THEORY OF THE INTERACTION OF A TURBULENT 
V/AKE V;iTH A SUPEESOMIC FLOW IN A CANAL AND JET 


By 

L.V. Gogish and G.Yu« Stepanov 


Interactions of the turbulent layer of a flow with the outside flow are 
observed for outside flow around aircrafts, flow in canals, • nozzles, diffusers, 
gratings of the blade engines. These interactions deterndne to a considerable 
degree the properties of these apparatus and installations. 

Starting with the work of Krokko and Lis [l], the investigation of the inte3>* 
action of the turbulent layer with the outside flow is performed by the integral 
methods of the theory of the boundary layer. These methods are used together with 
an assumption of the profile of the velocity in uhe tui’bulent layer. They secure 
a satisfactory coincidence of the calculated pressure distribution in the region 
of the interaction with the experimental data. 

Due to the complexity of the phenomena, occurring diu'ing the interaction, the 
most simple of these phenomena are studiexi experimentally and theoretically, such 
as the flow around a step (projection) with u plane supersonic flow or the interaction 
of a shock with a turbulent layer, when the outside How is represented by a simple 
wave (Prandtl-Maier flow) . The respective theoretical solutions in wldch the outside 

O nonviscoUB flow is described. by the PrandtW'laier functions ere actually one 

dimensional, which considerably decreases the calculation difficulties connected 
with their derivation. 

Hovrever, it is found that these simple solutions are insufficient for application. 
Great attention is paid to many practical problms, for examHe to the 
calculation of the interaction of flows arising from the flow around a step 
v;ith a supersonic free Jet (Fig. l 3 .)or in a canal at its sudden widening (Fig. lb). 

In these cases, in the outside nonviscous flow, a complex and previously unknown 
stimicture of shocks and rai'efaction waves exists- Such a structure 'interacts with 
the layer aiH determines its grovfth on the whole length up to the so-called choking 
cross-section (wake throat) . The use of ;xact methods for the calculation of such 
nonviscous ilows is associated v;ith a big v;aste of time and effort and is inadequate 
for an approximate (essentially not fine. dimensional) description of the flow in the 
layer . 

A simple approximate method of the calculation of the flows, represented in 
Fig. 1 a and b, is discussed below. In this method, the average quasi-one-dimensional 
method of the description of complex nonviscous flow &'id the method of the one- 
parametric turbulent layer or wake* are used simultaneously. 


* The wake differs from the more general case of the layer by the absence of 
friction on the wall limiting the flow. 


C 
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To complete the discussion-, the previously obtained solutions of 
the more simple limiting problem of flow around a step by a steady supersonic flovf 
in the presence of a single shock or a rarefaction wave, incident on the neighbouring 
v;ake, and a blowing into the bottom region [2^], are briefly described. ' 

It Was ascertained on the basis of the experimental data that, in the 
separating flows (Fig, l) directly behind the end, an isobaric bottom region is 
arranged; the flow of interaction is strictly studied, starting from a certain sucbitrar^ 
"lacing" cross-section I up to the choking cross-section, the wake throat II, 
behind which the whole flow can be considered supersonic. In the limits of the, 
integral methods of the boundary layer theory appliCcihie here, one-parametric 
family of solutions is obtained. These solxitions describe the flov? of the interaction 
at a finite distance from the end. The condition of the selection of the real 
solution lies in the fact that this solution must be a singular one (must pass by 
the saddle point, corresponding to the physical throat of the v;ake) » 


1 . The One- Parametric Turbulent Wake 


In the general case, the turbulent wake is characterized by the profiles 
of the velocity and temperature. These profiles comprise a finite number of free 
parameters n;^ (i =« 1, 2,..,.), the values' of which determine the character of flow 
in the wake. As a result of the calculation of the interaction, the functions of 
the_ longitudinal coordinate x , the characteristic width of the wake b Cx) , the 
parameters of the profile of velocity n.(oc), the parameters of the^outsiae 
nonviscous flow on the boundary of the wake, the value of the velocity C 
and its inclination to the X axis must be deteimiined in the plane of flow. 
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In the integral methods of caJ.culation, the turbulent wake is defined only 
by the integral values, genera3J.zed by the dimensionless v/idths of the displace- 
ment H* , the loss of momentum and the friction (K = 0,’ 1, 2, ....)• 

For the determination of the unkno;^! I'unctions, different integral relations of 
the boundary layer and other auxiliary equations, and in particular the conditions 
on the axis of symmetry of the wake, are used [2, /+]. The existence of a certain 
arbitrariness in the selection of the number of the parameters and the form of 
equations leads to a certain difference in the results, specially, if the number of 
the equations is small* 

In the practice_of the calculations in reference [4]f a one-parametric profile 
of the velocity u « u (y) in the free turbulent layer and wake v/as taken. 




I'-m/fi ) , = 




, m - 





(l.l) 


The indices " S " and "0” refer to the parameters on the boundary of the 
wake y = b and on its axis <^-=0 , respectively* 

As most simple equations, describing the flow in the one-parametric turbulent 
wake, the equation of momentum and the condition on the axis of symmetry are taken: 



( 1 . 2 ) 


( 1 . 3 ) 


in which j a 0 in the plane case, j = 1 in the axisymraetric case, 

C-U-IU-nYKi ~ M T Is the reduced velocity of the outside flow 

on the boundary or the wake , 5 ' . The turbulent, tangential stress X' is 

expressed by the Prandtl formula. 

, X -const.- 


Taking account of (l.l) 


7r-p. u. arm 


2 2 dj-h) f 


dh) 




z m 


!'C 
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here a = F/d ^the value (-‘^16 )is a free^. parameter (0.012 for the 
neighbouring plane wave). The parameter K = K where 





is the reduced 


width (area) of the displacement^ 


tw reduce.) 


width (area) of the loss of momentum; all are tabulated functions of the 
parametric form ra of the profile of velocity in the wake and of the numbers 
of C or M on its outside boundar-y. 

It is suitable to select either the absolute vddht , or the integral 

width of the displacement, as the characteristic width of the turbulent 

wake. These represent the principal parameter of the interaction of the viscous 
and nonviscous flov/, since the surfaces y « 8^^^ or 5-S’’'Cx) can be considered 
as lower boundaries for the outside nonviscous flow. 

Equations (l«2) and (l.3)» after transformation, take the form: 


SEn H 
3 m 

dm ^ dfn H 
dx 7 d£nC 


2.1 dEnC . . 

dEn 0 * 
dx , 

(1.5) 

dm 

rrt(Z- m) 

1 

U m 

2V 
(d*) 2 

(1.6) 

dx 


dx 

2* 


In the two equations (l.2) and (l«3)» there are three unknown functions of 
x; ra, G and cl- 

In order to close the problem, it is necessary to build an additional system 
of equations for the nonviscous flow. For this purpose, the equation of continuity 
in the form of Krokko-Liz [l] is used for the most simple case of an isoentropic 
plane Prandtl-Maier flow. Tliis equation has the form: 


d 

dx 





(lo7) 
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where ^ - - ’)fc) is the angle of inclination of the nonviscous flow on the 

boundary of the viscous layer “bjVCC) is the Prandtl-Maier fanctiono 

Taking into account that {S- S*) , the equation (l.?) can be 

written in the fonn 


d X 



dir. fpu)s 

dx 




( 1 . 8 ) 


In this case equations (l.5)» (l.6) and (1.8) form a closed system; its solutions 
describe the flow of interaction in the piano turbulent neighboui-ing v/ake. 

Equation (l«8) differs from the knotm Prandtl equation; 




(1.9) 


which is usually applied to perform more accurate calculation of the velocity 
distribution on the outside limit of the boundary layer. It is possible to show 
that equations (l.8) and (l.9) agree with the accurate boundary layer theory. 

From the point of viev; of calcuilation, the difference between the equations 
of Krokko-Luz (l.?) and Prandtl (1.9) is connected only with the selected boundary 
'j j:; o C,x; or , for which the parameters of the nonviscous flow are 

given . 


2 • Flow Behind the Step in Unbounded Flow 


In references [4i 5]» a system of the equations (l.5)t (l*6) (at j O) 
and (l.8) is applied for the investigation of plane flov: in the turbulent neighbouring 
Wake behind the step (Fig. 1 c). Solving this system of equations relative to the 
derivatives of the functions and C C—") we Obtain; 


dm ^ ^ 

dx t ' dx ' a ' 7^ 


( 2 . 1 ) 
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where 


S(l-m) r m(l-m) 



1 - .A 

' i 3 


dtnl4** 


m 


,3EnH** 

^ + / + 

3lnC 




^ mfz~m) dEnU 

(-Cx)rr.^H* /-C‘ ^ dP„C ' A'* ^^^7 






(-a,:cjr 


■H‘ 


fl-m]. 


dtnMl 

~dm J 


The system of equations (2.l) is reduced to one differential equation 
conceming m « m (c) 


^ 

' 4 . ( 2 . 2 ) 


and the two first integrals. 

At the given parameters Coo and equation (2o2) has a singular 

solution, passing by the singular saddle point. The coordinates of this point 
in the plane (ra, c) are determined by the equations ~ - 0 

(Hereby, simultaneously, -A. ^0 J) . As was sho^m in the work [/+], 

the singular solution corresponds to the flow in the undisturbed neighbouring 
wake, and the singular point corresponds to its throat. The nonsingular integral 
ciun/-es of the system (2ol) or the equations (2.2) describe the flows of the 
interaction of the turbulent wake with the supersonic flow in the physical plane 
or in the plane (ra, c) in the existence of isentropic shocks or rarefaction waves, 
disturbing the neighbouring wake. 

The initial conditions for the system (2.1), x = "X , J"- ^^d 

C »• C®, are determined as a result of the gluing in cross section I (Fig. l) of 
the flow of the interaction in the wake with the isobaric flov; of the displacement 
in the bottom region, on the basis of certain integral conditions. The 
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continuity of the absolute ;d.dth d of the layer and the displacements b 
and also the conservation of the mass of the gas in the bottom region were taken 
in the works 5] as such conditions,. In this case, the initial value of 
the pajrametric form m°, the length of the isobaric bottom region ■.>c° and the 
initial width of the neighbouring wa>:e are determined as a function of the, 
beforehand, unkno\m bottom press'ure p® and the given blowing parameter B ..by the 
equations : 



(2.3) 




(2.4) 

'■ G(m,C-) 0 



(2.5) 


In equations (2«3) - (2.5) j B and G are denoted as follows: 


g ,C%osO’) 'x‘] , 


(P^]ba 




where G (m>c) — the relative rate of the returning flow in the wake fo < •> c j 

B — the blowing parameter, S the empirical coefficient of the expansion of 
the isobaric turbulent jet. 

The value of the bottom pressure p® ( and its corresponding velocity C° 
of the outside flow) is found from the condition of the realization of the 
singular solution of equation (2.2), for which the initial condition m® is 
determined from equation (2.3). 

For simplicity of calciiLation, it v;as assumed that there is no initial 
boundary layer in the step. The location of the outside boundary of the isobaric 
mixing layer > 0 ;/ith respect to the boundar)' streamline of the equivalent 

nonviscous flow (Y « O) v;as simply determined, as for the ordinary boundary layer, 
by the formula (for the unbo'onded flow behind the step, from the 

condition of the conservation of momentum ■=% - i in the ^ 

bounded flow, this condition is complicated and for supersonic velocities (^d « d 
it can be disregarded. 



30 - 



If the equation of continuity (l»9) is used in the form of Prandtl, then 
the continuity of the v;idth 'b'' ' of the momentum loss must be used instead of 
the condition of continuity of the absolute width ^ of the layer (taking account 
of the returning flows) © 




( 2 . 6 ) 


Then equation (2.3) vrill have the form: 


I-B 


(2.7) 


Formulae (2.4) and (2.5) will be written in the previous form. 



The value of the free parameter essentially constant turbulence 

depends on the assumed equations, gluing conditions and the profile of the 
velocities. The above-mentioned value of ') “ 0.012, vjas selected in 

work [4] from the comparison of tha calculated and experimental data with respect 
to the bottom pressure P® behind a plane step in the condition of gluing' of ^ 
and with a profile of velocity (l.l) at -t ^ '2^ [ 'ij - 

For plotting the dependence of the relative bottom pressure p°/ behind 
a plane step in an unbounded supersonic flow on the Mach rromber M as of the 
undisturbed flow and the Ijlowing parameter f Cf ^ 

sufficient to use the family of the singular integral curves of the equation (2.2) 
and the family of the solutions of equation (2.3) or the cunres of 
the initial conditions mj_(C, B) . The intersections of these curves determine the 
two-parametric relation or p7P.^( represented in Fig. 2, 




C 
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The calculated curves of Fig* 2 satisfactorily agree with the _experimental 
data in the case of existence of blowing ( pa > o') or suction ^0)^ 

however, for clarity'’, the results of experiments (taking account of th e dispersion 
of the data of the different authors), related only to flows with j^iL~0 
are given in the graph, i*e*, in the case of absence of blov/ing (suction)* The 
curves in Fig. 2 are limited at certain maxima of blowing and suctions* The first 
limitation is related to the accepted assumption of the isobaric flow of the 
displacement in the zone of separation and, in principle, it can be removed if 
the interaction of supersonic and secondary subsonic flov/s is studied* The second 
limitation is more essential and is related to the existence of a limiting line for 
the solutions of the system (2*l) when the denominator in the right sides of the 
equations vanishes J 


A (m, C) ^ 0 


or 


ofn M f" M f II rn fz-m) d£n M 

Which indicates the impossibdJLity of the existence of separating flows,' if the 
flow rate exceeds a certain critical value* 

In Fig. 3j the calculated distributions of the pressures p = pC"2;") on the 
axis of the neighbouring, wd-ce, are compared with the experimental data at different 
Mach numbers of the supersonic flow in front of the step: 


/ - = 1,8(f 

3 - ^ = 2,30 

5 - /M„=3,02 

The numbers correspond to experimental data from the five different sources 
shown in the work [4]» 

The family of the nonsingular integral curves of equations (2*l) or (2,2) 
corresponds to the different flows in the neighbouring wake* These flows are 
disturbed by a shock or a rarefaction wave at finite distances from the beginning 
of the Wake* They are represented by segments of nonsing.ilar integral curves, 
which can be studied directly, or Jointed vjith segments of other singular curves* 



32 - 





p The flow in neighbouring wake, disturbed by a shock wave of finite" intensity 
c^s > ;i has two values. For given boundary conditions of the. Mach number 

in front of the step, the intensity of the shock '^s and the place of its 
location in the vrake two possible flows exist in the neighbouring wake. 

These flows are characterized by the values of the bottom pressure p° and 
which can be arbitrarily called the strong and weak disturbances. This corresponds 
to hysteresis in physical flows of sv.ch a t:.'pe. In the case of existence of ■ 
blovfing or suction in the bottom region, the flow in the neighbouring v/ake is 
two-valued in the whole range of the values of the bottom pressure < p‘‘< 

The flow in a v/ake with blov/ing in the bottom region can be studied as quasi- 
steady, when the wake acts as a source or drain for the subsidiary flow rate, 
and the bottom region becomes a storage for the mass. On this basis, a series of 
change of the parameters of the bottom region was performed in the work [5]» 

This region supposedly characterizes the lov;-frequency oscillations in the wake. 
Such oscillations are really observed, and x.heir amplitude and. frequency in one 
example satisfactorily agreed with the calculated values. 


3* Quasi-One-Dimensional Egautions of the 
Super sonic Jet 


For the approximate quasi-one-dimensional method of calculation of plane 
and axisymmetric thin Jets, having a small curvature, of a nonviscous gas, ell 
the parameters f of the flow ( P ^ f the longitudinal u and transverse \r 
velocity components) are assumed continuous and slightly changing across the 


C 
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cross-section of the Jet| x ■ const* At the boundaries of the jet the parameters 
at ° 'jo j marked by and ^ respectively* The variation 

of all tthe parameters /-/o and also the value of v are considered small of 
the first order* 

Applying the operation of averaging to the analytical functions of the 
parameters of flow, it can be ascertained that the mean value of a function, 
correct up to the squares of the small quantities,' is equal to-the same function 
from the average values of the arguments* In particular, for the mean parameters 
with the mentioned exactness, the usual gas dynamic relations of one-dimentional 
flows hold} at the same time, taking into account the small value of v, it can be 
assumed that / ■^ 

~ fdx = ju-^ w !<u > (i~ 0,i) etc. 


The differential equations of motion after averaging with respect to v are 
reduced to the ordinary ones with respect to -x, vthich is the main purpose of the 
developed approximate method. On the basis of a special investigation and of 
the practice of calculations, it v/as advisable to use the following 

system of equations (the signs of averaging are omitted) : 


e ^uatian 


or continvdty 



equation of momentum in the longitudinal direction 


dx 





(3.1) 


(3.2) 


equation of mbraehtum in the transverse direction 




(3.3) 


equation of the boundaries of the jet 


dyddx-v-Ju , L-.pj (3,4), (3.5) 

Equations (3.3), with the above«mentioned accuracy, can be divided into two 
equations with respect to dv-i /dx. and cLv'Jd-^^ (This physically corresponds 
to the division of the jet into two and to the application of each of them of 
the equation of momentum at the step change Cf/;, ). 
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Performing this division and solving the equations with respect to the 
derivatives, we obtain, after some transformations, a system of six equations. 
This system describes the gas flow in the nonviscous curvilinear jet. 


'dyjdx-.vju , dyjdz-ifju , 

(3.6), (3.7) 


(3,8) 


(3.9) 

dv-„_ r. .1 ^ ° 

dx~ ^o]f n[y/*'-yJ*’) 

(3.10) 


(3.11) 

In the written equations, C = is 

(^ •= U Ic^ } a, is the velocity of somd, ^ 

the reduced velocity, 



T^e 'coefficients ‘=^{, in the case of a circular jet srei ^ .^^£y )l3(y, -t-yt,)- 
^ilz^ilz , o^,rf2y,^yj/zfy,^y,)--2f3^fl2j 

and for a plane jet *^o - - i/x. The system of the six equations (3»6)-(3«H) 

is complete, because it contains eight functions of X, two of which(. ’/ t; and 'rT| 

or i/o and TTJ ) are given. The initial conditions are known in the cross- 
section DC— O' 

The obtained system of equations is used for the approximate construction 
of the boundaries of the nonviscous curvilinear jets 9^ Cx; in the sections, 

where there are no strong shocks inside the jet and the angles of inclination of 
the boundaries to the axis of the jet are small. It is obvious that the greatest 
error in the determination of the boundary. of a jet arises near the edge of the 
nozzle and the points of fall of the shocks, in which the boundary of the Jet 


\indergoes a sharp bend. Therefore, the reduced system can be used for the 
approximate calculation of the plane and circular curvilinear supersonic jets 
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at small degrees of unrated outfLovr. Examples of contours of plane curvilinear 
supersonic jets are shown in Fig. 4* These contours are obtained by means of the 
integration of the system (3.6)-(3«ll) . Examples of contours plotted by the 
method of the characteristic features for the most simple con*figuration of the 
jet| when C^ = are given also in Fig. 4» ( is the reduced velocity 

in the first" section of the jet). 





In order to emphasize the wave (barrel-shaped) configuration of the jet in 
the quasi-one dimensional calculation, it was assumed that the nozzle has a small 
expansion angle 10<^ = O,3-,In the correct construction it was assumed for 
simplicity that ~ 0 . It is obvious that on the average, the boundaries of 

the correct and quasi-one-dimensional jets satisfactorily coincide. In the case 
of the most sharp bend of the boundary (Fig. 4 b) the quasi-one-dimensional jet 
is limited by the section 3-5 where M *» 1. 

As applied to the calculation of the flow of an interaction, we notice that 
only the integral properties of the complex outside flow up to tlie throat of the 
neighbouring wake are of interest. Examples of such properties are the number 
and the dimension of the "barrels’* of the jet, which remain in one-dimensional 
approximation in a "spread" form. Moreover, the angles of inclination of the 
velocity vector of the outside flow on the boundary of the turbulent layer in the 
whole region of the interaction do not exceed 20®. So, the case of Fig. 4 b 
does not belong to the calculation of such a flow. 
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The discussed quasi-one-dimensional method of calculation is distinguished 
by sufficient simplicity! and it is similarly applicable for subsonic and supersonic 
nonviscous flows, both in the thin bounded jets and in the narrov/ canals. For the 
developed theory of interaction, it is particularly important that, by using this 
method, the problem as a whole obtains a closed analytical description that does not 
contain unnecessary informations about the flow and convenient for the programming 
of the nimierical calculations. 


4* Interaction of a Jet with a Wake 


The flow of the interaction of a quasi-one-dimensional jet and one-parametric 
turbulent wake is described by a system of six equations for the nonviscous jet, 
(3»6)-(3»11), and two equations for the wake, (I.5) and (I.6). In the last two 
equations, one should assume that and ~ 

necessary to add an auxiliary fomnila to the system '^f these equations for the 
determination of the length (surface) of the outside and inside boundaries of the 
nonviscous jet. 


This is necessary for the calculation of the width of the jet turbulent layer of 
the mixture. 

The initial conditions of the neighbouring wake are determined, as previously, 
from the conditions of the gluing of the isobaric flow of the mixture in the bottom 
region with the flow of the interaction in the wake. On calculating the turbulent 
layer of the mixture behind the axisymmetric end, one should proceed from those 
approximate relations that exist in the plane flow; the profile of the velocity 
in the layer of the mixture is universal 

• If = m . (4.2) 


the area width of the cross-section of the zone 
of the mixture, 6 mixing coefficient 


Q 0 

is the area (length) of the ejector svirface of the jet. 


where ^ ' Q* ^ 


(4.3) 
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The conditions of continuity of S and in the section of gluing and of 
conservation of mass in the stagnant region are written in the form: 


y‘-'(x}.S\(x}[l-H'(l.C.o]] = S*" , (M) 

, ( 4 ^) 

, ( 4 . 6 ) 


Where 




i (m) 


The relative, integral widths are expressed as functions 'of three arguraentsi 
the parametric form mf the reduced velocity C and the parameter J. • The last 
argument, as shown previously, is equal to 0 or 1 in plane or axisymmetric layer, 
respectively. 


From equations (4.4)“(4«6), it follows that the initial value of the parametric 
form m® is determined as a function of C from an equation analogous to (2,3): 


hH 'O.C',0) 

SCm.CM) ■ U"(t,C’,0Xt-B) ’ 

The coordinate “X. -'x" of the section of gluing is determined from the equation 



8(nH*(mc:i] 




= D . 


(4.8) 


The calculation of the flows of the interaction, which are represented in Fig. 1 a 
and b, consists of two parts. Originally, for a given value of the bottom pressure 
P® (or C®), the flow on the isobaric region behind the step up to the section of 
riuing, determined by equation (4*8), is calculated by integrating the system 
(3.6)-(3*11)» complemented by equation (4»l). Further, by the integration of the 
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;r same system of equations^ the flow in the region of the interactions is 
calculated, , The real value of the bottom pressure P“ or G° and the actual 
flow of the interaction are determined, as in the most simple case in work 
[4], by the condition of passing the, integral curve through the singular 
point (wake throat). The numerical constnaction of the singular solution 
in the neighbourhood of the saddle singular point ^whose position is only 
determinecl iq the ope^r^tion of the calculation jip connected vdth known 
. calculating difficulties, t In .practice^ it seems sufficient to construct 
two nonsingular solutions related tor.the different families. The difference 
of the values_gf the_bgt_tom pressure in them is',' beforehand, less than the 
given value, 



The results of a successive integration (by the method of Rungo-Kutt) 
of the described system of equations for a plane flow at - 3.0, Si ° 1, 
fef " 2, 60 =--0,1, 6,, » 0.1, for different values of the relative 

bottom pressure R,/p^ (P„ is the pressure in the outside medium) are 

represented in Fig, 5» The integral curves in the plane (m, c) are shown in 
Fig, 5 a. These curves characterize the change of the parameters in the viscous 
layer. The contours of the nonviscous jets are shown in Fig, 5 b. The lower 
b-oundary of the jets in the region of the interaction corresponds to the 

width of the displacement of the viscous layer 





The nonsingular integral curves, given in Fig. 5 a, do not correspond to 
any actual flovu All the curves begin on the line of the initial data (in the 
absence of blowing, B ■' O), which is determined by equation (2,3). The right 
branches end on a limiting line which is determined by the condition A(mjC) = 0. 
The left branches are limited satisfactorily conditionally by the separated 
section tJC-- 10, beyond which the solution has physically unreal singularities. 

From the study of Fig. 5 a, it is obvious that, at given initial conditions, 
two singular curves are found. These two singular curves are plotted simply as 
averages between the nonsingular curves, related to the different families, 
and until these nonsingular curves are sufficiently close to each other. The two 
stationary flows of the interaction corresponding to them are represented in Fig. 6. 
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As in the case of the intjaraction of a single shock with a wake, a flow 
with lower bottom pressure, » 0*375 end a flo\i with higher bottom pressure 
m 0»55> can be called weak and strong disturbances, respectivelyo The 
transition from one flow to the other is connected with the essential change 
of the configuration of the region of the returning flows ^d of the location 
of the critical points (m » l) in the wake from x » 5 (at ■ 0.375) to 
X- 8,2 (at 0*55). 

The existence of the double-value in the flow of the interaction in the 
considered, more complex, case also attests the possibility of quasi-st at ionary 
ascillations with a relative amplitude j 


^(prp)l(p° . 


5 • Pseudoshock in a Canal 


As shown in rectilinear long canals, a complex system of weak oblique shocks 
and flow separations, interacting with the turbulent layer near the wall, arises 
instead of a normal shock in such a way that the transition from a supersonic to 
a. subsonic flow in 'the canal and the corresponding loss of the’ total pressure in 
the flow are distinguished by a viscous mechanism, which is commonly called pseudo- 
shock. Following work [7]| we assumed that the effect of the shocks on the flow 
in the core of the canal, outside the layer of mixture, can be neglected, and this 
flow can be considered one-dimensional and isentropic. Vie assume also that the 
friction on the wall can be disregarded. The corresponding sketch of the pseudoshock 
is represented in Fig, 7* 





In section 1, the flow consists of a uniform core and a thin ordinary turbulent 
boundary layer. This layer achieves a separating state under the action of the 
positive pressure gradient. Between sections 1 and 2, a layer of mixture with 
returning flows of fluid grows on the wall. In section 2, the layer of the mixture 
absorbs the isentropic core, after that a balancing of the flow takes place in the canal 
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A second critical point, behind which the asymptotic increase of the pressure to 
the maximum occurs, is located in section 3 when taking account of the friction on 
the walls in a certain section ''4”, the pressure reaches a maximum [8], 

Since the effect of friction on the wall is not taken into account, the pseudoshoc 
can be considered as a flow of interaction of one-dimensional outside flow with a 
layer of mixture (with a wake). 

In the layer of mixture between sections l-4»a one-parametric profile of 
velocities of - ordinary form (l.l) is received. The three unknown functions of the 
problem in the region 1-2 will be the reduced velocity cC'x) the width of the boundary 
layer S Cx) and the parametric fora m (x), and in the region 2-4, where S 
will be the pressure p (x), we have as previously C(x) and m (x) . The original 
system of equations consists of two equations of conservation! the equation of 
continuity isj 


- const 


(5.1) 


and the equation of raonentum isi 

d 




(5.2) 


dx - ' 1 “ 

and also the differential equation on the axis of the viscous layer as: 

( 5 * 3 ) 


^'P . /'ll) 

dz * dx ^ ^ 


y '' 0 

In the case of the isentropic flow, we have in the core of the flow (in the 
section 1-2):. 




For the calculation of the final increase of pressure in the pseudo, shock, 
it is sufficient to use the equations of continuity and conservation of momentum. 
It is obvious that the mean parameters in the final and initial sections of the 
pseudoshock are connected with the relations of the normal shock. 

r 

For the calculation of the pressure distribution in a canal the equations of 
conservation, (5.l) and (5.2), should be used in the differential fora. 
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The original system of equations of the pseudoshock for the region 1-2 
with the isentropic core can be reduced to two ordinary differential equations of 
the first order with respect to m (x) and C (x) : 


dm / : - 3€n H \ 

' 3m / da: ( ^ 1- c^^i(^(c) ' 96n C ^ 

dm mfZ-m) dtnC _ (~ax) 

dx dx 


dinC „ 
:: 0 

(5o5) 

dx ’ 



1 - m 

( 5 . 6 ) 


Equation ( 5 . 5 ) determines the integral curve ra (C) in the plane (m, C) 

^ - 1. d* i* (c)]. SinH/d^nC ^ 

de c /// 


From equation (5*6)# the pressure distribution can be calculated in the physical 
plane C (x) or tt: (x) j 


• Cj" ' ' ’ ■ ■■" 

(-ax)(x - X, ; =/ [i - 


dm m(Z-m) 

n-c^]c 


1 


dC 


( 5 . 8 ) 


Equations (5<>5) and (5»6) are valid up to section 2 of the coupling viscous layers, 
where - 1, or 





( 5 . 9 ) 


For the calculation of the balancing flow, one should proceed differently. 
The dependence of -m on (Co is the reduced velocity on the axis of the canal) 
is more conviniently determined in this case from the transcendental equation, 
which is obtained from equations (5»l) and (5»2): 



l-HY«^.Cj-ir(m,C,} t-cl I 



(5.10) 
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The value of the constant is determined from the data of integration of equation 
(5«7) in section 2. 

The distribution of the velocity C, C 2 :) is obtained from equation (5*6) • 
The distribution of pressure PCxO/p^* is calculated by the equation of 
momentum : 


P*' 


(5.11) 


For the simplification of calculations, the integral parameters of the turbulent 
layer of the mixture were assumed to correspond to the linear profile of the 
velocity / = i-^ in the nonviscous fluid: 

Then from equation (5*7) the following equation is obtained: 


dm _ f y _3 • 

dT ~ ’clj \1 


(5.12) 
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The integral curves of equation (5»12) at » 3 in the plane (m, C) are 
represented in Fig, 8 a« They are valid in the region 1-2 up to the section 
of the coupling of viscous layers at ^ 

The different curves correspond to different conditions in the initial 
section of pseudoshock 1 or to a different width of the displacement of the 
initial boundary layer ♦ This consequently leads to a change of the initial 
parameters of. the flow in the nonvlscous core and of the mean-mass velocity of 
the flow in the section (Table l) * 

The cui*ve 'b «» If corresponding to equation (5t9)f is now determined by the 
equation; 

• (5*13) 


The integral curves in the region 2-3-4 were plotted by the formula (5#l0)i 
taking account of the assumed simplifications. 

The distribution of the pressures P (3^:^ and the parametric forms m (x) 
for three variants of calcixLation are shown in Fig. 8 b, It was assimied that 
iDi^ ■ 0,85 in the final section of the pseudoshock. The obtained parameters 
of the flow in the 'pseudo shock and the compai’ison of its length with the 
experiment in work [8] are given in Table 1 (the linear dimensions are relative 
to the hight hf stated in Fig, ?)• 

Table 1 


EapnajiT paenoTa 
(puc, S) 

1-4 

l ' - 4 ' 

l " - 4 " 

ha|ja;.K!-:-pu nuToioa e ua-io-'iuitOM 
coauiiKit nceBAOCKUHKa: g) ^ 

h 

0,064 

0,212 ' 

0,353 


0,016 

0,053 

0,088 

^ mi 

1,02 

1 ,76 , ■ 

1 ,5 

CE) 

Jl/iHHa ncfi&nocKaHKa: pacneT ' 

22,5 

36 

29 

oitcuepuMeuT 

22,5 

38 

29 


T) Variant of calculation (Fig, 8)j 2) Parameters of the flow in the initial section 
of the pseudoshock j 3) Len^h of the pseudoshock: calculation experiment. 
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It is obvious from the table that there is a totally satisfactory agreement between 
the experiment and the approximate calculation (at one and the same value 
of mi^ - 0,85) . 


6. The Final Remarks 


In conclusion and first of all, one must emphasize that the essentially developed 
one-dimensional theory, in spite of the sufficiently rough initial assumptions and 
with the minimum number of the empirical constants, totally satisfactorily describes 
all the main properties of the investigated separating turbulent supersonic flows, 
including the hysteresis phenomena and the quasi-stationary oscillations. The 
good agreement of the calcvlated and experimental data is remarKeil also in other 
analogous works. This is explained, .cn . one hand, by the existence of certain 
universal properties of the separating flows, slightly depending on the general 
configuration of the flows and enabling to use its simple calculating models, and 
on the other hand, by the smoothing characteristics of the integral methods of 
investigation. These methods guarantee at the same time the conservation of the 
principal properties of the phenomenon. 

The first universal properties of the separation of the turbulent boundary 
layer and the combination of the turbulent wake at definite critical pressure 
ratios P,,/P-i which are, with a high precision, functions only of the Mach 
number of the incident undisturbed flow, were ascertained experimentally in the 
works of I.P. Nekrasov, G.I. Petrov and other authors. The properties of the 
separation, taken phenomenologically, enabled to construct more complex flows by 
the investigation of the local separations and combinations of the turbulent flow 
in the case of achievement of the conditions of the critical increase of pressure 
or the critical angle of turn of the supersonic flow. These conditions obtained 
satisfactory explanation within the limits of the theory of nonviscous fluids, 
taking account of a definite reconstruction of the velocity profile of the 
separating or adjoining vortex: layers and the satisfaction of the integral 
conditions of conservation of mass and momentum. Another explanation was connected 
with the assumption of conservation of the total pressure at the separating line 
of flow. This pressure, in the case of combination, must be more static for the 
magnitude of the increase of pressure in the region of combination. This assumption 
enabled Korsft and Chepraan to determine the bottom rarefaction 'behind the step in 
the plane flow, including the effect of the low blowing, in the case of good agreement 
with the experimental data. However, as the more detailed investigations in works 
[9|3] show, the method of Korst-Ghepman is essentially connected with the local 
properties of the selected profile of velocity and of the flow of the nonviscous 
fluid along the line of flow of a constant rate, which must be considered the 
principal disadvantage of the method. 

Further progress in the theory of the separating flows is due to Krokko and 
Liz [l(] who extended the integral method of the theory of the boundary layer in 
the viscous fluid to these flows, taking into account the separating profiles 
of the: velocity ond the essential interaction of the layer with the out$ide 
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supersonic flow. Unlike the previous approach t this approach enabled to study 
in detail the growth of the layeri to obtain the pressure distribution and the 
location of the points of separation and combination, and to show the effect of 
the initial conditions, disturbances and friction in the flow. 

The use of the integral methods of the theory of Jets and wakes represents 
at present the most promising direction in the investigation of the separating 
flows. 

Among the most important problems for further working out, one should point 
out the extension of the range of the investigation of Mach number M and Reynolds 
number Rej the consideration of the effect of the initial boundary layer and the 
friction on the walls} the study of the flow aroimd more complicates space 
configurations; the introduction of consumption and thermal effects; the changes 
of the design model and working medium, corresponding to the conditions of the 
hypersonic velocities at high altitudes and to the cavitation flows of the 
incompressible fluid. 
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OPTIMUM FORMS OF PLANE AND AXISYMMSTRIC BODIES AT 
HYPERSONIC VELOCITIES 

By 

A*L* Conor 


The solution of a variational problem in a statement | using the resistance 
law of A. Buzeman, was proposed in the works [l-3]» HcJwever, as Kheiz had pointed 
out in work [2], in the more exact statement, the contour of a body of a minimum 
resistance must have a discontinuity of the inclination of the tangent at the end 
point. This is so because, according to the law of A. Buzeman, an infinite negative 
pressure is created at this point, which decreases the resistance in the final value. 
Physically, the pressure can not be negative, and the change of the inclination 
of the tangent at the end point in the supersonic flow must not ahead affect the 
pressure distribution and, consequently, the resistance. The mentioned discrepancy 
with the physics of the supersonic flow requires a new statement of the variational 
problem. 

There are two possibilities. One of them is the use of .the so-called the 
thrust ring and leads to the conception of the absolutely optimum body [2,3]. 

The second possibility involves an additional requirement of limiting the class of 
the bodies so that the pressure is positive everywhere on the contour. With such 
a statement the variational problem was strictly solved in work [4], for bodies of 
a given elongation and in work [6] for thin bodies with different boundary 
conditions. A number of questions, connected with the statement of these problems, 
is investigated in work [5]« The solutions of this problem in case of arbitrary 
isoperimetric conditions and without limitation on the tJiickness of the body are 
cited below. 


• Problems with Arbitrary Isoperimetric Conditions 


The solution of the variational problem about the form 6f a body of a 
minimum drag with a given length and diameter does not cover all the possible 
cases- of interest for practice. 

We shall take in consideration, beside the length of / and diameter, the 
following quantities: in the two-dimensional case— the area A, enveloped by the 
contour of the profile and the moment of inertia M of the contour; in the axisym- 
metric case— the wetted surface S and volume V. Their values are equal to:: 


A = 2j ydx 
f-c 

Xt 




I 


1 


= 'yVx ; 

Xp 

X, 


( 1 . 1 ) 
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Any two values are most frequently given} while the rest are considered free. 

For example JL and V} ¥ and S} t and S etc. The general problem is formulated 
as follows: find the minimum drag Xq in the condition that the pressure on the 
contour is nciv-negative and the values of (l.l) are given. The principal theorem 
of work [4] about the existence of a section of zero pressure on the optimum 
contour is, according to the argument method, not related to the form of the 
boundary conditions} therefore it can be extended without difficulty to the 
general case. 

As a resultjthe problem is reduced to a detection of the minimum drag 



nil 

r j :y“ 1 




(1.1) 


in the condition of fulfilling, on the length x ^ the 

inequality; 


^ J (u u'yd 


(1.3) 


Here ^ - is the dynamic pressure} n - 0 and 1 for the plane and axisym- 

metric flows, respectively. The variable X- is the coordinate in the direction 
of the incident flow; y- in the perpendicular direction. The index 0- refers to 
the initial point, index 1- to the end point (Fig. l). The drag arises only 
from the section OC of the contour. On the arc OF, called the free layer, the 
pressure is equal to zero. 





C. 
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Let us turn to the investigation of a set of curves of zero pressure. 
On this set, the inequality (l<>3) converts to equality and represents a 
derivative from the formula. 



(1.4) 


Hence, it follows in particular that it is impossible to carry out a contour of 
zero pressure from the initial point X » 0, 


Eliminating the integral from the formulae (l»3,) and (l,4), we obtain a 
simple differential equation of the second order for the function i) (x) 




n*f 


y y *xyc y ’ 


(1.5) 


where 'Jc — is the ordinate of the corner point, and the parameter K is found 
from the formula 



( 1 . 6 ) 


As it is evident from (l*6), the geometry of the free layer depends on the 
geometry of the Initial section, called the regular form. The general integral 
of equation (l*5) contains two arbitraiy .constants, which can be so determined 
that the free layer has no inclination in the point 
and passes by the point ^ a result we obtain the equation 


. n*? n*7 

<Ji - y 


n*!/K / \/ 1 n*^t. 


(1.7) 


detennining the geometry of the free layer. 

Applying (l,7) to the corner point, we obtain the relation 

-- — — -y’ */— — ~]/yf'y)"^y” • 


8 ) 


Taking into account (1,7) i we rewrite the conditions (l,l) in a more suitable 
form: 

i/zrr-J 

V07=/ 

■*0 

where J, V,,.,,- are known functions of their arguments. 
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Let us formulate the auxiliary functional 


*° rU il *' 


(1.9) 


Where At -.- constant numbers. The minimum functional of (l.2), in the above 
enumerated conditions! coincides with the minimum of (1.9). Hence the unknown 
solution must satisfy the condition 'b|^=0- Expanding the latter, we obtain: 


Af C 

Sy^Hf (F^~dfy'/dxjSt/dx*SHj Fdx *h[(P ~ yFy')Sx + 
*0 • 

* F y'SyJ f Sd ^ 


1 

1 

c 

yr 


' n*1 (1 


u( n, 

f1~n) 


n) 

■a 


m(K, 


Repeating the usual argument we find that the extremum satisfies the Euler 
equation, which assumes the first integral, 


<3 ^ 

+ r— rr J(/ ^ (i ' » i jij/ ~ C ~ 0 . 

The condition of transversality is represented in the form 


( 1 . 10 ) 


IhJ I- dx * H[c Sx * Fy SyJ -i-SdiO 


( 1 . 11 ) 
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According to (l«6) and (l.lO), the coefficient K depends on the arguments 
(^0 J ‘^o'^sequently 

The form of integral (l*10) is such that the equation of the extremum in the 
general case can be written in quadratvires. Actually! solving (l.lO) with respect 
to the variable ^ and taking into account the parameter y ' f vie obtain 






(lol2) 


Writing (l.l2-l) in the initial point* it is possible to express i/o by 
the parameters in the point of conjugation and in the point If equation 

(1.12-1) is now applied to the point of conjugation and C is eliminated! then the 
following relation will take place 


.i'c .^c > 


(1.13) 


The variation (l.l3) gives the condition 





(1.14) 


The condition of transversality (l.ll) is reduced to an analogous form. 
Investigating it together with (I.14)! is possible to eliminate any one 
variation and! by (1.13)! 'the corresponding variable! for example! Xc • As a 
result we obtain a condition connecting four variations with coefficients that 
are dependent on four variables. Finally! the variation of the equation of the 
free layer (I.8) enables to decrease the number of the independent parameters to 
three. The latters are determined from the equality of the corresponding 
coefficients in the independent variations to zero. The discussed shceme gives 
the possibility to find the values of Xc ^Vc j • Th e 

constants 'X.: | entering in these parameters, are calculated from the relations 

(1.1) , Thiis, the solution of the general problem is closed. 

If the extremum does not pass by the origin of the coordinatesj when the class 
of the bodies with plane nose [5] is engaged in the investigation. The condition 
of the transversality gives an auxiliary relation for the elimination of the 
variable 
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^2, Investigation of the General Solution 


On performing the solution, certain moments, which require special investigation, 
were omitted. First, it was not proved that the inequality (l.3) will be realized 
for the regular form of the body. Second, it nowhere follows that the optimum 
contour contains a free layer only at the end and that it cannot be composed from 
some sections of the free layer and regular form. We shall carry out the proof 
of these situations for the case of a body with a given relative thickness. The 
correctness of the condition (1.3) is most easy to verify by a direct calculation 
of the pressure distribution on the nose of the contour. In the plane flow, the 
regular form is the wedge and the condition (3) is certainly fulfilled. In the 
axisymmetric flow, the pressure distribution is determined by the following formula: 


Cp- pIc} = sin^d- — %ind * -I 


lin^oC 


8 CCS d 




The corresponding graph, represented in Fig. 2, shows that a positive pressure is 
always realized in the regular form. 



Let us now determine the succession of the arcs of the regular form and free 
layer. As already shown, the extremum cannot start from the free layer. Therefore, 
the initial section is the regular form or the arc of the boundary extremum. One 
may ask: does the transition from a regular form to a free layer occur directly at 
the end or we have intermediate arcs of a free layer and regular form by means of 
which the conjugation of the initial and final sections takes place? If the latter 
situation is assumed, a reverse transition from the free layer to the regular form 
will certainly take place at a certain point. Such point is denoted by the letter C 
in Fig. 3; the arcs of the free layer are shown by the dashed lines and the arcs of 
the regular form by the solid lines. 
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To reduce the calculation, we shall limit the investigation to the plane 
problem. The drag coefficient of a body consisting of two arcs of the regular 
form (Fig, 3 ) is determined by the formula; 



Calculating and equating the total variation (2,l) to zero, it is in particular 
received that at the point C, a condition of transversality of the follovriLng form 
must be fulfilled: 


8nf Fdt* H[(F.-ij'Fy‘)Sx*-Fy'SyJ^ ^ 0 




. r r 




On calciilating the first term, it is necessary to take into account that the 
integral of the function F is taken both on the regular form and the free 
layer. As a result, (2,2) will have the form; 




^0 


(2.3) 


On the other hand, the equation of the free layer (l,7) reduces to the 
relation; 


yc-y/, 




C 
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Hence, the variations at the point G are related with the following condition; 




(2.5) 


It is easy to 
inclination of the 
indicts "1" to "C" 
variable X, Then, 


set also a relation including the variation of the angle of the 
tangent at the point C. Therefore, changing beforehand tl 
and ”C" to "Cj”, we differentiate (l»7) with respect to the 
we obtain the formula; 





( 2 . 6 ) 


the variation of which gives the condition: 




(2.7) 


Equations (2.3)', (2.5) and (2,7) form a linear one-dimensional system with respect 
to the variations. Therefore, the determinant from the coefficients must become 
zero. This leads to the relation;: 


yp *-2y y^' 0-v'y2)[y(f^ yp : 


We shall rewrite the latter in the form; 


‘ r ^,7^^ / . (2-8) 

If we now consider (2,6) and (2.8) together, we have at the 

point C, which contradicts the initial inequality In this way, it is 

proved that the free layer can not occupy the intermediate position and be alternated 
with the regular form. 
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^ 3 • The Second Method of the Solution of Variational 
Problems with Inequalities 



C 


In work [6], the above discussed problems are reduced, in the case of thin 
bodies, to the solution of general problems of a variational calculation of the 
Mayer-iagrange type. The principal idea of the solution is borrowed from work 
[7J and includes the insertion of an auxiliary function, enabling to replace the 
inequality, by an equivalent relation. We extend this method to the case of bodies 
of arbitrary thickness. The strong and weak sides of the solution will be seen 
in the course of its performing. 


Let us consider the general problem of minumum wave drag in the case of the 
inequality (1.3) and arbitrary initial and isoperimetric conditions. We shall 
write the formula of the drag in the fonii: 




ft* 



(3ol) 


According to [?], we replace the inequality (1.3) by the relation: 




A 

I 


'■■"dt 


y y 






(3.2) 


where p is a real variable 


Differentiating (3»2) with respect to the variable X and eliminating the 
integral from the both relations, we obtain: 




y'[n(l*y'^)+i/y 




(3.3) 


Introducing the designations: 

y\t ; y"»3 ; = . 


the relation (3*3) will be equivalent to the following system: 


_ fm ti 



0 


i fjit/'-t * 0 

is 0 

f s i'-ms 0 
■* 3 


(.3-4) 
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AS a result I we come to the formulation of the Lagrange problem; in the 
class of the functions t mCv), pCy^ compatible with the conditions 

of the relation (3*4) i the isoperiraetric conditions (l.l) and the conditions on 
the end points | we find such a system of functions, which would be reduced to the 
minimum functional Qf (3«l)» 


We set up the auxiliary functional: 






(Ut^yk 



^0 

* ly(Un]y ^ , 


( 3 . 5 ) 


The vanishing of the function 7 gives the following necessary conditions: 
1. The Euler equations:: 


takes place along the optimum contour. 
2* The conditions of Weierstrass-Erdman: 


'HI I 

are achieved at the comer points. 

3u. Finally - the cond^i tions of transversality are: ^ — 

, her.”. 

Xo’O, <,= /, , yys, , yg'-P • 

Expanding the Euler equations ^3.6;, we obtain the system: 






. f * y " " 

- + ^5 fjf tJj'4 


- =0 


(3.6J 


( 3 . 7 ) 


( 3 . 8 ) 
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(3.9) 


The absence of the independent variable enables to write out the first 
integral: 


( 3 . 10 ) 

Equation (3.9-3) enables to draw the conclusion that in the general case, 
the extremum can be discontLnuous, and it consists of arcs, along which t 



t s * ^ 

y 




or 


p^O 


( 3 . 11 ) 


We shall n ame the first arcs th e regular forms jand the second arcs yie_free liters o 

The conditions of (3.1) in the corner points necessitates that the 
multipliers , X 3 should be continuous and the constant C, 

of integration, should have the same value for all the arcs forming 
the extremum. Moreover, 


(P^o). -(P^o)+ 


( 3 . 12 ) 


The '.form of the conditions (3*8) depends on whether they are recorded in fixed 
points or in points of the natural type. In the last case we have: 


'■0; ; 


(^3 >133 


If the end points are fixed, then the first three conditions of (3.13) disappear. 

We shall illustrate • further solution on the example of a plane contour with a given 
relative thickness. 


The first equation of (3*9) in the case gives the integral 
to S! -giSar fern' and the free layer, according 
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The forms of equation (3 til) and of the boundary condition for the Lagrange 
multiplier enable to suppose that this multiplier is Identically equal to 

zero along the regular form. Such supposition gives the possibility to solve the 
system (3*9) and satisfy all the boundary conditions# Actually, the second equation 
is written in the form: 

^ ^ ifc^ Of''- ' , 

This solution satisfies the system (3 #9) at aivy values of ^3* 

We obtain the equation of the free layer from (3*3) at P - 0. Performing the 
integration, we have: 

The constant • The functions ^ and ^3 along the 

arcs of the free layer are found from the solutions of the second, fourth and fifth 
equations of (3*9) at the conditions (3 #12) and (3»13)» 

The composition of the arcs forming the optimum contour needs a special 
investigation, similar to that performed in ^ 2. 

In conclusion, it is possible to draw the conclusion that both methods give 
an identical solution# Nevertheless, the second method, even in the most simple 
case, leads to the necessity of the investigation of a complex system of 
differential equations, the solution of which is performed in a heuristic way. 
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BODIES OF MAXIMUM AERODYNAJOC QUALITY IN A HYPERSONIC FLOU 

By 

V.I. Lapygin 


In the present work, the form of a conical body having the maximum value of 
aerodynamic quality is determined. It is supposed that the pressure distribution 
is determined, by the Newton formula and that the coefficient of local friction is 
constant f the length and area of the mid-section of the. body are given. However, 
such statement does not give the possibility to perform an analytical solution. 

An essential simplification of analysis was attained by auxiliary limitations for 
the class of the permissible surfaces. Till nov/, the procedure of Strend [8] 
has obtained the maximum spread. Strend solved the first problem about the form of 
a minimum-drag wing that is thin in both the longitudinal and transverse directions. 

Subsequently, the optimization of the mentioned class of surfaces at different 
auxiliary conditions was investigated in the works [l, 9, 10]. In particular, the 
problem of a thin wing of maximum aerodynamic quality was discussed in detail in 
works [l, lO]. The simplicity of solution in such a statement is attained owing 
to the reduction of the three-dimensional problem of the form of the wing to the 
finding of the optimum profile in the longitudinal cross section. On the other hand, 
it is known that the form of the cross section of a wing strongly affects the value 
of aerodyneunic quality [2, 11, 12]. Therefore, a more exact statement of a 
variational problem may give essentially new results for the optimum form of a 
hypersonic wing. In the below- cited solution of the variational problem of the 
form of a conical wing having a maximum quality, the assumption of the thinness of 
wing is removed. 

The solution represented below was discussed initially in work [6]. At the 
moment of its completion, the investigations [l3, 14], in which similar results 
are obtained, became known. However in contrast to work [l4], the results of the 
given work are obtained for a more general statement of the problem, not limited 
by an auxiliary assumption of the form of the lee of the wing. Work [l4] 
involves the same trend. In this work, for the first time, the method of the local 
variations is. used for the variational problem of a wing in the Newtonian 
approximation. 



1. Derivation of the Calculating Formulae. Statement of the 
Variational Problem 


shall consider the hypersonic flow of gas around a conical body. Let the 
axis of the cylindrical coordinate system j be selected along the 
flow, the velocity of which is U (Fig. l). The surface of the conical body in 
the assumed coordinate system is given in the form: 


I 


(l*l> 
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We shall consider that the total drag of the body consists of two parts: 
the wave drag and the frictional drag. We shall neglect the portion of the 
frictional force in the lift force. We shall determine the pressure distribution 
on the body by the Newton formula 


Cp- ~kcQi^(n, U] ^ 


( 1 . 2 ) 


where C p - the coefficient of pressure on the body, relative to the dynamic 
pressure of the incident flow. 


K — the proportionality constant . 
n - the normal to the surface of the body. 

Below, we shall assume that k = 2, ^ = 1. 

We assume that the value of the local coefficient of friction is constant 
along the surface. With the assumed allowances, we obtain the following formulae 
for the coefficients of the aerodynamic forces (in view of the symmetry of the 
problem with respect to the plane ^ « 0, the half of the body, ^ TTj 

is investigated): 


SC. 


-Si — ^ zT l 


wViecf 


ar. I ct, 


fl'v’ 


f 


(1.3) 


(1.4) 




^•is the area of the stern cross section of the body. 


C-c 


'■is the local coefficient of' friction . 
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On investlgatir^thin bodies, these formulae are simplified. Actually, at 
"t introducing a new dependent variable, X H we obtain: 


o 




!/" * 


(1.5) 

( 1 . 6 ) 

(i ^7) 


The coefficient evidently does not enter in the right sides of formulae 
(l»5) and (1.6), which simplifies the investigation. 

The statement of the problem:; to find the function y (cf>) at the given value 
of (lo7)i along which the ratio 


. ( 1 . 8 ) 

takes the maximum value. 

The coupling between ^ and the aerodynamic quality of the body is given 
by the relation 


2 . The Algorithm of the Numerical Solution and the Results 

of the Calculation 


The analytical investigation of the problem encounters considerable difficulties 
connected with the fact that the Exilerian equation in this case is essentially 
nonlinear and has a second order. 

A direct method is used for the solution of the problem. It is known by the 
name "method of local variations" [6, 3]| in which the isoperimetric condition, 

• const, is fulfilled in every step. The convergence of the method in the 
case of absence of the isoperimetric condition is shown in work [5]. Because 
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is maintained constant in the process of calculation, it is possible to hope 
that the used method coincides with the Eulerian equation for this case also; 
the calculations confirm this assumption* 

The calculations were performed on the E7M- BESM-3M. The interval of cp 
Was selected to be 4°i the initial value of the variation of the radius h - 0.01. 


y 



For checking the convergence of the method to the unique solution, the contours 
shown in Fig. 2 are selected as a first approximation. These contours have the 
same area, but they differ in the form and the position of the maximum radius 
We denote the angle corresponding to the maximum radius by . It was found 
that the angle tp'*' at a sufficient big number of interaction (h in the process of 
the calculation is decreased in 8 times), remains as it was in the initial 
approximation, with the exception of the case of the triangular contour. The 
forms of the contours are obtained qualitatively similar. For the explanation of 
this condition , some contours were selected. These contours are close in form 
to the optimum, and have the same area, but they differ in the angle ^ . It 
was found that the dependence of quality on in the vicinity Of the optimum 

angle is very weak (Table l) . 

Table 1 

12° 64° 00° 56° 48° 36° 


"k 0,3904 0,3877 0,4018 0,4015 0,3963 0,3739 


Analogous dependence of ^ on takes place at other values of ' 

The fact of variance of the angle at ■ const is evidently due to the 

insufficient accuracy of calculation (at a moderate number of interaction). 
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One must mention that the velocity of the convergence to the solution 
essentially depends on the selection of the initial approximation. 

In the calculation! the value of the aerodynamic quality is determined 
considerably more quickly} than the contour of a transverse cross-section} 
and} starting from a certain moment} is changed extremely insignificantly. 
Therefore the obtained dependence A (S-^) has a high accuracy; the forms of 
the optimum contours are found with less accuracy. We notice that in the 
process of the calculation} a decrease of the interval of was not performed, 
because the time of the calculation considerably increases in this case, and the 
accuracy of the determination of A is practically not changed. 





The dependence of the optimum aerodynamic quality K at ■ ICT^ on the 
parameter of the volume Z where V— is the volume of the body, and 

S- is the area in the planejis shown in Fig. 3» The maximum values of K for 
V-shaped wings [2] are given by the small circles in the graph. The value of K 
for the V-shaped wing practically coincides with the value of K for the optimum 
body. We notice that these points correspond to wings with a cross-section close 
to a triangle, the base of which is directed to the windward side. 

Let us cite data on the value of the. aerociynamic quality of the bodies, 
investigated in works [ll, 12]. At a given length of a half-cone body with a 
cross-section in the form of a semicircle has a quality K *« 3.6 at - lCr3 
and m 0.0727} with an elliptical cross-section, the maximum increase of the 

quality is 2fo] with a sine cro3S-section~-l2.5^; with a triangular cross-section 
with the lateral sides facing the flov/, the quality approaches the value for 
optimum bodies; for all the remaining cases, the increase is ‘noticeably less. 

If we consider the triangle cross-section, with the base turned to the side of the 
flow, then, as the calculation shows (Fig. 4)i the quality of the optimum body 
will practically coincide with the qu^ity of the body of the mentioned form. 

We notice that the triangular contour, with the base facing the windward side, 
consists of the extrema of the Eulerian equation for the investigated problem. 
However, this contour is not the solution because in the place of conjugation of 
the two arcs of the extrema (half-line and straight line, parallel to the OZ axis), 
the condition of Weierstrass- Erdman is not fulfilled. The forms of the optimum 
transverse contours as a function of are represented in Figs. 5 and 6. 




_ AA _ " 



It is noteworthy that the whole volume of the body is limited from above by 
surfaces representing! within the limits of the calculation errori two planes 
located along the flow. The lower surface of the wing, receiving the whole load, 
has a small convexity to the side of the incident flow. In the vicinity of the 
angle ^^=71 there is a very small elevation, the hight of which decreases with 
increasing the number of iteration. The obtained numerical results confirmed with 
high accuracy the hypothesis that the upper surfaces of the optimum wing are planes 
located in the flow. As a whole, the form of an optimum transverse contour is close 
to a triangle, with the base facing the windward side. Hence, it becomes clear _ 
whey the quality of a wing with an optimum triangular cross-section almost coincides 
with the quality of the optimum. 


L 
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For checking, a calcvdation__was performed wi th the use of formulae (l.3) 
and (1.4) for S S - 2?T . lO*^ ; S = 45T . io‘^ at a = 10“3, without 

an assumption about the thiimeas of the body. The results are obtained identical 
with the abov/e- mentioneti. In particular, the form of the body and the value of the 
aerodynamic quality for the same values of S are practically the same. 


3 • Remark About the Convergence of the Solution 


The problem vuider considerat ion is equ ivalent_tQ -the problem of the absolute 
extremum of the functional J J /<pc{<f 

V ycos * y'si" V’ . f r~i } 2 


Ji » const , 

If the obtained solution will satisiy, with the given accuracy, the Eulerian 
equation for the function 0 and the boundary conditions, then the numerical 
convergence of the method does not arouse doubts. In order to be convinced ty this 
it is necessary to check: 

1. The fulfillment of the boundary conditions at the extremum, 

2. The fulfillment of the Eulerian equation. 

3. The fulfillment of the conditions of the conjugation of extremum arcs in 
the comer point* 

A direct numerical check of these conditions was performed for '5^ -7r 

1. The boundary conditions at the points ^ - O , =TT (the condition of 
the transversality) ~ 0 are fulfilled exactly, because 

yriTj^o. ^ 

2o The fulfillment of the Eulerian equation was checked only in the interval 
^ 0 < because the half-line of ^ is' an extremum. At every 

point of the contour the value of was calculated in the interval 

The results of calculation are represented in Fig. 7» For 0< ^ 

the deviations from the mean value do not exceed 6.5?^; in the vicinity of 
the variance of the values is slightly increased, which is explained by the ^ 
difference of the actual angle from the assumed angle in the calculation y- 

(As shown in Table 1, 56® < *pj- <.Go ) . 


<7 r 
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3» The conditions at the corner point: 

The condition =■ 0 ±g satisfied vd.th an accuracy of 3?^;- the condition 

4=^3 =0 is satisfied with an accuracy of 10^. 

Here AC--J denotes the difference of the values on the left and right 
from the corner point. 

In conclusion, we mention that the accuracy of the method, as clear from the 
result of work [4]» has the order of the interval of . In the performed 
calculations, the interval of ^ was equal to O.O698. Thus, the obtained above 
valuations of the convergence of the solution coincide with the accuracy of the 
method. 
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SUPERSONIC FLOW OF AIR AROUND A RECTANGULAR PLATE 
By 

M.P, Falunin 


A considerable quantity of vrorks is dedicated to the aerodynamics of vdngs 
of different forms at high velocities. An extensive bibliography of this 
problem is comprized in the works [1-4]. 

In the present work, the pressure distribution on the windward side of a 
rectangular plate of infinite elongation is experimentally investigated, a 
study of the spectra of the flow around is also carried out on wide angles of 
attack for a wide range of change of elongation values ( 0 , 025 ^ X ^ 40 ) . 


§ 1. Pressure Distribution on the Windward Side of an 
Infinite Rectangular Plate 


The study of the pressure distribution on a rectangular plate with cutters 
on the ends, decreasing the influence of the end effects, permitted the 
determination of' the diagrams of the pressures along the chord of the plate, and 
also the location of the center of the pressure of the plate. The' pressure was 
measured in three cross-sections in the central part of the plate. These cross- 
sections are located at a distance of 15 mm from each other. In each cross- 
section, 6 drains v/ere located. The distance between the drain points along 
the cross-section was equal to 5 mm; the distance from the edge of the plate 
was equal to 2,5 mm, A sketch of the location of the drain points is represented 
in Fig, 1, 



Fig. 1. 

Key! 

l) Distance, 
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The experiments were carried out at an incident flow velocity correspondipg 
to a Mach number of M = 3,0. The obtained results shown in Figs. 2- 4 indicate 
the following; 


1, The pressure at the drain points located at the same distance from the 
leading edge is the same in all three cross-sections. Therefore in the 
graphs, the pressures are shown at one point for every value of Ti (the distance, 
along the cf^rd from the leading edge relative to the length of the chord) , 

2, At narrow angles of attack ( 30®), the pressure along the chord 

practically does not change and agrees with the theoretical value for a 
rectangular plate of infinite speed at the corresponding angle of attack 
(Fig. 2). 

3, Beginning from o/. 35®, when the leading wave withdraws from the 

leading edge, the pressure at that edge obviously increases, and by increasing 
the angle of attack, it approaches the value of the stagnation pressure behind 
the normal shock at M «= 3» The pressure at the trailing edge slightly increases. 
The pressure distribution has an unsymmetrical form (Fig, 3| I- V), 


4. By increasing the angle of attack, the zone of the increased pressure 
shifts from the leading edge to the middle of the plate, and at o 90 *^ 
the maximum pressure corresponds to h «= 0,5 (Fig. 3| Vl), The pressure near 
the zone of the increased pressure begins to decrease in proportion to the 
displacement of this zone from the leading edge. At the angles of attack 
35®-§i <( 90®, near the trailing edge (at a distance 0,1 of the chord), 

the pressure is less than that at the leading edge at the same distance. As 
the drop between the maximum pressure on the windward side and the pressure 
on the leeward side of the plate is supercritical, the velocity and pressure 
at the edges themselves must be critical. ‘ , 


M • 3,0^ A 


o 3KcnepLiMeHm 
pac7em0 



Fig. 2. 


1) Experimentali2) Calculated 




erimentaL2) Pressure behind the normal shock 
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5. Ovd.ng to the change and displacement of the region of the increased 

pressure, the center of the pressure (in the condition that considers the 
pressure on the leeward side is equal to zero) also shifts from the center 
line of the plate to the side of the increased pressure. The graph of the 
relation Gj is represented in Fig. 4* The data are obtained for the angle 

of attack such that 0° <=^^90®. It is obvious that as long as the 

leading wave is attached to the leading edge, the pressure on all the faces 

is the same, and the center of the pressure is located on the central line of 
the plate (0^ *= 0»5). After the detachment of the shock wave from the leading 
edge and the appearance of the zone of the increased. pressure near it, the 
center of pressure shifts forward by about ^ ~ 0.46; ~ 35- 45®)- 

On further increase of the angle of attack, the center of pressure shifts to 
the middle of the plate, at o>c = 90°, = 0,5. The increase of the angle 

of attack from 90° to 140 ° leads to the displacement of the center of pressure 
to the trailing edge •?:; 0.54) • The attachment of the shock wave to 

the trailing edge of the plate necessitates the decrease of C ^ to 0.5. Thus, 
the relative coordinate of the center of pressure of the forces acting on the 
windward side of the rectangular plate may be changed for the range of the 
angles of attack 0° ^ 180° by about 15- 17?o. 

6. The dependence of the relative coordinate of the o? Imum pir<'ssuy'« 

on the angle of attack is represented in Fig. 5. In Figs. 4 and 5, the curves 
for the angles of attack, 90 °-:^ ex- 180°, are drawn on .the basis of the ^ 

symmetry of the flow around the plate at angles of attack ^ and 180° ...xc . 
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2* Spectra of Supersonic Flow Around Rectangular Plates 
at Lar^e Angles of Attack 


The study of the spectra of supersonic flov^ around rectangular plates 
was carried out within a wide range of change of the angles of attack from 
20® to 90® and of elongation from 0.025 to 40. 

General remarks on the form of the shock wave: The analysis of the 
spectra of flow around and their quantitative treatment point to a sharp 
qualitative difference in the form of the shock wave* for small and big 
elongations. First of all, for such characteristic features of the form 
of the shock wave, it is possible to refer to the following (Fig, 6); 



* All the data shown in the sketch graphs and tables are relative to the 
vertical plane of symmetry. 
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a) The presence of two straight segments of shock wave in the case of plates 
of small elongation and narrow angles of attack; one segment is attached to 
the leading edge, the othei>— at a certain distance from it. In the case of 
a plate with an elongation of the order of one or more, there is only one 
straight segment adjacent to the leading edge, b) At wide angles of attack, 

when the shock wave is totally detached, the shock wave is curvilinear in the 

case of plates of big elongation while it is, for the most part, equidistant 
from the windward side of the plate in the case of plates of small elongation. 

The values of the quantities from which it is possible to judge 

the form of the shock wave and its variation with the angle of attack ■=*^ and 
the elongation ^ & /h are given in Table 1, where ^ is the length 

(spread) of the plate, and h is its height (chord). All the quantities are 
relative to the height of the plate h. The quantities in Table 1 correspond 
to the following geometric characteristic of the shock wave (Fig, 7). 



Fig. 7, 


l) Xn , 2) — the distance from the leading edge of the plate to the shock 

wave along the plate, upstream and in the perpendicular direction, respectively; 
3) 'jf— . the distance from the center of the plate to the shock wave 
perpendicular to the plate direction; 4) X^ , 5) the distance from the 

trailing edge of the plate to the shock wave along the plate, downstream and 
in the perpendicular direction, respectively; 6) , 7) — the coordinates 

of the end of the straight segment of the attached shock wave with respect to 
the leading edge; 8) , 9) — the coordinates of the end of the 
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curvilinear segment of the detached shock wave with respect to the leading 
edge; lO) X*, ll) U* - the point grid reference of the shock wave with 
maximum departure from the surface of the windward side of the plate, with 
respect to the leading edge. 

For illustration, the graphs of the quantities 1- 11 (except 6 and 7) 
are represented in Figs. 8 and 9 for a plate with an elongation A >=> 0.0375. 
The numbering of the curves corresponds to the enumerated quantities. 

Model No. 2 = 0.0375 ; h = 80 mm> M = 3.0 
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Key; 

1) Curves. 


We have to do some explanations in connection with the quantity X,- : 
a) if the curvilinear segment of the shock wave, adjacent to the leading edge, 
converts afterwards to a rectilinear one, which, in its turn, is already 
curved because of the effect of the trailing edge, then denotts the 

distance from the leading edge to the projection of the end of the first 
curvilinear segment on the windward surface of the plate; b) if the shock 
wave is totally curvilinear, then the quantities X and X* coincide; 
c) in the case of a symmetric flow around a plate {cx = 90) , the quantities 
of X ^ and X* coincide. The dash in Table 1 denotes that the corresponding 
quantities could not be measured in the limits of the visual field. This 
refers to the quantities X , X^, , X^, y,, , X*, In individual cases, 

this denotes the absence^ of truthworthy data about the corresponding 
quantities (X^ , 
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Model No. 2; 0.0375 i h = 80 mm; M = 3.0 





Kp: 

l) Curves. 


Analysis of the straight segment of the shock wave;: The straight segment 
adjacent to the leading edge of the plate corresponds to flow around a plate 
with attached shock wave. In this case, in the vertical plane of symmetry, the 
straight segment of the wave proportionally increases with increasing the 
elongation of the plate, until one achieves the elongation of the length 
at which disturbance from the trailing edge reaches the wave more quickly than i 
from the lateral ends. At M e 3.0 and angle of attack ^ « 30°, the 
proportionality is broken for "X > 0.75, and at the angle of attack « 20° — 
for A > 0.125 (see Table 1 and Fig. 10). Llnt-i 1' j :k' in Fig. 10 
correspond to the experimental data for the relative length of projection of 
the straight segment of the shock wave on the plane of the plate. The solid 
lines v/ere calculated with Mach disturbances from the lateral ends of the leading 
edge'. 

The results of the analysis of the spectra of flow around, indicate that 
the rectilinearity of the attached shock wave in the vertical plane of 
symmetry is broken at values of that are smaller than implied from the 

mentioned calculation. 

Dependence of the form and dimensions of the leading shock wave on the 
angle of attack and edge effects; The angle of attack -and the elongation of the 
plate essentially influence all the geometric structure of the shock wave, 
namely the length of the straight segment, tlie magnitude of vdthdrawal, its 
departure from the leading and trailing edges of the plate, etc. 
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Key; 

l) Calculated; 2) Ejqjerimental. 


The corresponding data are given in Table !• In the case of narrow^ angles 
of attack, when the shock wave is attached to the leading edge, plates with an 
elongation exceeding a few units •wHl te l-loy^/ec! ar©tjiid tht j-t-.'-trsl gcirlro?. t\\e 
windward side by a translational flow along the whole height. The shock wave 
has one straight segment, which is curved because of the disturbances running 
from the trailing edge. This bend of the wave is observed d'omstream at a 
considerable departure from the plate. The slope of the straight segment 
corresponds to the calculation. In the case of plates with small elongation 
( A< O.l), the structiire of the shock wave at certain angles of attack is 
different; the translational flow occupies a small segment adjacent to the 
leading edge; the shock wave at the edge is rectilinear, then it is strongly 
curved, converting to the second straight segment, which is almost parallel to 
the face of the plate (Fig. 6a). 
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At angles of attack, occupying an intermediate location between the normal 
flow arovind ( = 90®) and different angles of attack, the sharp difference in 

the form of the shock wave is retained for plates of different elongation 
when the wave is attached to the leading edge. In the case of plates with an 
elongation A 0-5 the shock wave in the plane of symmetry is totally 
curved. The distance from the windward side to the shock wave in the perpendicular 
direction is variable: it is minimum at the leading edge, and maximvun at the 
trailing edge. It is characteristic for plates of small elongation ( A < 0,l) 
at the mentioned angles of attack that, for the most part, the shook wave is 
equidistant from the force; at the same time, the distance between them is 
approximately the same as in the case of symmetric flov/ around ( o<- ■= 90° , 
see Fig, 6b and Table l) , 

At last, the difference in the form of the shock wave for a plate with 
small and big elongation, in the case of flow around perpendicular to the 
windward surface is shown in Fig, 6C, Plates with small elongation have a 
straight segment of shock wave on the whole extension, with the exception of 
the range adjacent to the edges. In connection with the influence of the finite 
dimensions of a plate on the form of the shock wave^we should add that in the 
case of narrow angles of attack, this influence may be extended to considerable 
distances. At angles of attack ^ o 6O-90®, the ends of the plate exert an 
influence on the form of the shock wave at a departure approximately equal to 
4- 3 times the magnitude of the spread (width) of the plate,' respectively. 
Therefore, if ' A > 8 or A < l/8, the middle part of the plate will 

(Exhibit EicTw as a plate of infinite spread. That is to say, ' the relative 
geometric parameters of the shock layer will not depend on the location of 
the fairing segment of that part of the plate. 

Determination of the magnitude of the withdrawal of a shock wave in the 
case of a flow around a plate which is perpendicular to the flow; We shall 
relate the magnitude of the withdrawal of the shock wave to the height (chord) 
or to a side of a square equivalent in area to the plate. If at a fixed height 
of a plate, its elongation is changed, the relation of the magnitude of the 
withdrawal of the wave to the height of the plate ( <5 ) gives a clear idea about 
the dependence of this relation on the elongation. In particular, in the case 
of elongation exceeding 6- 8, the quantity Z approaches a constant value. 

The corresponding data are given in Table 1 and in Fig, 11 (solid line). If 
we relate the magnitude_of withdrawal to the square root of the area of the 
plate, this quantity {S' ) will enable to reveal the equality of withdrawal 
of the shock wave for the quantities A and 1/A • The maximum is attained 
at A = 1 (Fig. 11, the dashed line). In connection_with the great range 
of change of the elongation of plate, the relations <? ( X ) and 2'' ( A) are 
represented in a logarithmic scale on the axis of A , According to the 
data (Table 1, Fig, ll) , in the case of thj change of the elongation of the 
plate from 0,025 to 0,125, the value of iS is linearly changed, and at 
A ^ 8, it remains constant and is equal to cs; 0,72, Thus, in the case 
of the change of elongation by more than two orders, the withdrawal of the shock 
wave changes forty times. The maximum value of S' is one half of 
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1^, 



Fig. 11 



On the basis of the given discussioni it is possible to draw the 
following conclusions. In the case of a supersonic flow around a' rectangular 
plate, the structure of the shock wave strongly depends on the elongation of 
the plate and the angle of attack. In particular, in the case of a perpendi- 
cular flow around a plane plate, the withdrawal of the shock wave may be 
changed by ten times. 

In the case of narrow angles of attack, the influence, of the leading 
and trailing . edges may be transferred to comparatively great distances downstream. 
At angles of attack close to = 90°, this influence is extended along the 

plate to distances exceeding its width by three— four times. Plates with 
elongation A l/S and A ^ 8 are flowed around in their central part 
as infinite plates. 

In the case of an attached leading wave, the center of pressure of the 
forces acting on the windward side of a rectangular plate with infinite 
length is located at the centeral In the case of a detached wave, the 

center of pressure is shifted to the side facing the flow. The magnitude 
of the displacement of the center of the pressure may attain 15- 175^. 

In conclusion, the author expresses sincere thanks to G.S. Ul'yanov, 

A.F, Mosin, A. A. Makshin and L.V. Filyand for their help in carrying and 
working out the experiments. 
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0,720 

1,220 

14,0 

90 

1,220 

0,626 

0,720 

1,220 

16,0 

90 

1,220 

0,625 

0,720 

1,220 

18,0 

90 

1,220 

0,626 

0,720 

1,220 

20,0 

90 

1,220 

0,625 

0,720 

i ; i 2 o 
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SUPERSONIC FLOVJ AROUND PENETRABLE PLATES AT NARROW 

angles of attack 

By 

G,S, Ul'yanov 


The development of the aviation and space technology raises a vdiole series 
of important scientific - technical problems in the stabilization, deceleration, 
descent and landing of different objects. For these purposes, parachutes of 
different constructions and other stopping devices acting at supersonic velocities 
are used to a large extent. The supersonic flow around penetrable bodies has 
a whole series of specific aspects and peculiarities, to which the works [l - 3] 
are dedicated. However, a whole series of problems are still open at the 
present time. For example, in the case of supersonic flow around a penetrable 
plate with a finite thickness under an angle of attack, the necessity arises to 
estimate and point out when one should take into account the tangential component 
of the aerodynamic force. For this purpose, it is necessary to determine 
experimentally the value of the tangential component relative to the normal 
component, its dependence on the degree of penetrability, angle of attack, Mach 
number M, and relative geometric dimensions. This enables to work out in a 
physically substantiated way the equations of conservation of momentum and energy 
for the penetrable surface, 

The present work is dedicated to experimental study of supersonic flow around 
penetrable plates at narrow angles of attack. 


1, Models and Techniques of the Experinents 


Squares plates of 100 x 100 mm with a thickness of 2,5 mm and symmetrically 
sharp leading and trailing edges, were investigated. The perforation of the 
plates was attained by uniform boring of holes of diameter 3 mm along the whole 
surface of the plate. The penetrability of the plates was ^s follov;s;: 0, 8,5, 

16.0, 25,5, , 51.5?^. The ratio of the sum of the areas of the holes Sq to the ' 
total area S is known as the penetrability factor W of the plate S 

(w = (s„/s; 100 M), 

The experiments were carried out in a supersonic wind tunnel at Mach numbers M 
of the flow from 1,5 to 3,0 £ind Reynolds numbers Re ^ = 2,5 ~ 4»10° ( ^ = 0,1 m). 
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The gravimetric investigations of the penetrable plates were carried out 
in the range of angles of attack from 0 to 17®. The plates were fastened in 
the working section of the wind tunnel on two belt suspensions. The aerodynamic 
coefficients and Cy were determined. The correction for the suspension 
device was determined by a separate experiment, in which the belt suspensions 
were blown through without plate. The calculation of interference had not been 
considered. The relative error in measuring the drag force X and lifting force ,'j 
was in the limits of 0,5 - 


2, Results of the Experiments and their Analysis 



We consider that the aerodynamic coefficients G and G ^ are determined, 
with the forces X and related to the total area of the plate, and the 
coefficients with index "P* (Gxp i C vjp )•— to the area taking account of the 
degree of penetrability. The coefficients C r\ and for the normal and 
tangential components of the aerodynamic force were determined by recalculation 
and were related to the total area of the plate. 

The relations GxpCt^) and at the ^ach number M o 2,5 and at' 

different values of the penetrability pairameter W are represented in Figs 1 
and 2, It is obvious from the given graphs that the resistance coefficients 

) nonlindarly increase with increasing the angle of attack, A similar 
pattern of monotonic nonlinear increase of the resistance coefficients with 
increase of the angle of attack was obsejrved in the range of Mach number 
M = 1,5— 3.0, 





C 
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Fig. 2. 


The aerodynamic coefficients of the lifting force ^ mp(C^ ) ^he 

Where ^ is in radiants and the coefficients Kc^ and.Ki^p are given 
in Table 1. ' 


Table 1 


C^-.K,^(W,tA).U 


M ' 1,5 

2,0 

W'/. O 8,6 16,0 25,5 51,5 

O 8,5 16,0 2^5 51,5 

Kf^ 2,522 2,064 1,662 1,204 0,401 

2,064 1,662 1,318 1,032 0,344 

K. 2,522 2,236 1,949 1,605 0,917 

2,064 1,834 1,548 1,376 0,745 

2,5 

3,0 

W*/ O 8,5 1^0 2^5 51,5 

O 8,5 IQO 2^5 51,5 

1,834 1,490 1,146 0,860 0,287 

1,433 1,146 0,86 0,688 0,229 

K, 1,834 1,605 1,376 1,204 0,631 

Up \ 

1,433 1,261 0,974 0,8599 0,516 
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The graphs of the aerodynamic quality K j- (o^) for plates 

vd-th different penetrability parameters W and Hach number M t= 2.5 are given 
in Fig. 3. It is obvious from these graphs that the quality of the penetrable 
plates strongly depends on the degree of penetrability. For example, the plate 
vd.th a penetrability W «= has a maximum quality approximately eight times 

less than that of the solid plate. The maximum aerodynamic quality for 
penetrable plates at a Mach number M = 2.5 is observed approximately at angles 
of attack .= 6- One should also notice that K slowly depends 

on the degree of penetrability. Similar relationships in the behavior of the 
aerodynamic quality are observed at other Mach numbers M. The values of 
K Loc) for Mach numbers M = 1.5— 3*0 and plate penetrabilities w = 0- 51/® 
are given in Table 2. 





Table 2 


w % 

0 

8,5 

16 

25,5 

51,5 


M = l ,5 

4.65 

3,55 

2,95 

2,05 

0,55 


2,0 

4,55 

3,45 

2,60 

1,80 

0,55 

2,5 

4,45 

3,05 

2;25 

1,70 

0,55 


3,0 

3,70 

2,80 

2,10 

1,55 

0,55 
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Fig. U 


The relations Cx (w) and G>j (w) at different values of the paoraraeter 
and Mach number M = 2.0 are shown in Figs. 4 and 5» The increase of the degree 
of penetrability of the plate at a fixed angle of attack strongly decreases 
the lifting force coefficient and slightly increases the resistance coefficient. 
The influence of the degree of penetrability on the resistance coefficient 
at narrow angles of attack is more strong than in the case of wide angles. For 
examplei for M = 1.5 and angle of attack = 5® » the value of C^p increases 
five times with the increase of the penetrability factor from 0 to 51 . 5 ?^f anci 
at = 15®— only by two times. One should also notice that the coefficient 
(w) is a linear function of W. As a consequence of linearity, the resistance 
coefficient Gx (W) can be calculated by the formula: 
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In the mentioned formula, the penetrability factor of the plate is taken 
in the form of the ratio ^ . The values of the coefficients 

and are given in Table 3. 



With increasing the Mach number M, the aerodynamic coefficients and Gy 
of the penetrable plates, at constant angle of attack and degree of penetrability, 
decrease in a linear fashion (Figs. 6 and ?)• 
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The relations C;, (M) and Gu (M) for different values of the penetrability 
parameter W and at an angle of attack -= 10® are shown in Figs. 6 ^d 7. 

For other angles of attack, the observed relationships are similar. The 
quantity G^ (M) can be calculated by the formula* 




Table k 



W/. 

0 

8,5 

16,0 

25,5 

51,5 

c/ » 5 


0,017 

0,040 

0,033 

0,027 

0,013 

V 

0,280 

0,240 

0,200 

0,160 

0,0X) 

lO 


0,127 

0,097 

0,093 

0,070 

0,008 


0,630 

0,480 

0,430 

0,290 

.0,080 

15' 


0,180 

0,147 

0,120 

0,003 

0,033 


0,910 

0,740 

0,605 

0,440 

0,145 


The polars Gy = f (G^ ), for penetrable plates at a Mach number M = 1.5, 
and Cv D = f (C ) , at a Mach number M = 3, are shown in Figs. 8 and 9. One^ 
should notice sol^e peculiarities in the behavior of the integral carachteristxcs 
of Cn and C.-r for the penetrable plates. In the case of the penetrable 
plates, with increasing the degree of penetrability and angle of attack, the 
coefficient of the tangential component of the aerodynamic force G.-c increases 
and approaches in Value the coefficient of the normal component 0,^ , and even 
exceeds it. For example, for the penetrability W = 51.5/», ^ 15 and Mach 

number M.= 3.0, is approximately twice . The coefficient C,^ ne^ly 

does not depend on the Mach number M and it is 9 linear function penetra- 

bility (W) of the plate. It can be approximately calculated by the formula: 


where K _ «= 0.32 for 

c o 


oc. 


15®-, 0.23- 10®; 0.19- 5®; 0.08- 0®. 




iinilM BMBffll 

g/^ BwS BBa 

.RSiSBBHBSManI 


flBmSBBlI 

BBli^Bifll , 

BBBgpBBBl 
flss^iissr 
KifilSBiflfll 
nr^BBR^i 


IB 
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The relations C;^ («x) and (c^) for different values of the penetrability 
parameter W at a Mach numlaer M = 2,0 are shown for illustration in Figs, 10 and 11. 
The values of the coefficients G and 0 for the mentdoned ranges of the Mach 
numbers M, angles of attack and penetrability of the plhtos, are given in Table 5» 
The obtained experimental data show that in the case of supersonic flow around 
penetrable plates under an angle of attack, it is impossible to neglect the value 
of the coefficient G •j;' , in comparison with G , The appearance of the tangential 
component of the aerodynamic force is connected with some physical peculiarities 
which take place in the case of flow around a penetrable plate under an angle .of 
attack, 



A1*2.0 ; W% a B,s IS 25,5 5/5 

o A O • 


The essential pecularities are the following : on the passage of gas through 
the holes of a plate, the velocity of the flow changes in value and direction 
and as a result, an auxiliary force arises. The increase of the thickness of 
the penetrable plate leads to an appreciable reconstruction of the flow and to a 
change of the pressure distribution on the walls of the hole^. 
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PAHT II 

MOVEMENT OF GAS WITH EXOTHEI^nC REACTIONS 

formation of plane detonation wave at the decay of discontinuity 

IN FUEL GAS 
By 

S.A. Medvedev 


The problem of the formation of a detonation wave in an Inflammable medium 
with a finite reaction rate and instantaneous energy release in a finite or semi- 
finite volume of reacting gas, is solved. 

The problem of the decay of the initial discontinuity in a fuel gas [l] is 
one of the problems whose solution is necessary for the answer of the question 
about the possibility of formation of detonation waves from. physically practicable 
initial conditions. The self-similar problem of the decay of arbitrary disconti- 
nuity, in the case when a heat supply takes place in infinitely thin fronts of the 
combustion or detonation waves, was solved in the work [2], In the present work 
we use a model of an inflammable medium in which the ignition decay is taken into 
account, and the heat release is described by relaxation equation in which the 
characteristic time of the heat release may depend on the gas-dynamical parameters. 
The finite difference numerical method is used for solving the system of equations 
of the gas dynamics and chemical kinetics. This method enables the direct 
calculation of the discontinuities and their interactions. The calculations, 
which were performed for the cases when the density of the energy supply in the 
case of instantaneous heat release (’’explosion*') is comparable with the 
density of the energy release in a gas with a finite reaction rate, showed that 
the intensity of the initial shock wave generated as a result of the decay of the 
discontinuity is less than the intensity of the leading front of the detonation 
wave in Chapman - Jouguet system. They showed also that the detonation wave in 
the undisturbed gas is formed as a result of the fusion of the shock wave with 
the heat-release zone formed on the contact surface, with the initial shock wave. 
From the obtained numerical solutions for the cases ' when the specific energy of 
the explosive gas is less or equal to the specific energy of a gas with a finite 
reaction rate and when the dependence of the rate of the heat release on the 
pressure is taken into account, it follows that the shock wave generated behind 
the initial shock, with a reaction zone behind it on fusion with the initial shock, 
is more intensive for small energies of explosion than a similar wave at high 
specific energy of explosion. The velocities and pressures behind the fronts 
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It consists of simple waves, progressive flows and crossed simple waves. 

The trajectory of the shock wave consists of rectilinear and curvilinear 
parts. If the rarefaction wave did not succeed to prevent the ignition 
on the contact surface, then during the expansion of the gas in the combustion 
process, compression waves (or shock waves) should be formed. These waves 
propagate on both sides from the place of ignition through the products of 
the explosion and through the heated gas of the initial shock wave. The 
detonation wave, generally speaking, may appear before the compression wave 
(or the shock wave of the reaction [8]) reaches the initial shock wave (and 
in such case may not be fully compressed) . Investigation of the flow, after 
co-flowing of the wave of reaction occurs with the initial shock wave, is 
of special interest. 

The model of the medium, which is used for solving the problem of the 
formation of detonation during the decay of discontinuity in a fuel gas, is 
set in the following syst em o f equations i 



dv ^ podP 
dt p 

^^LdivV = 0 
dt P 


Jr 

-- -aP e KT 
dt 


l(c-i)-{o. c^i; 

l,C>l} . 


Here V - the concentration of the fictitious component determining the time 
.of induction (O C l)f p —the concentration of the product of the reaction 
determining the heat supply Q. p(t) to the mass unit (04- 4- i)i Q*~ 
calorific power of the unit mass of the gas, P - the characteristic time of 
recombination, m and — constants characterizing the order of the reaction, 
a — factor, E— activation energy. The remaining signs are standard ones. 

As follows from system (l), the transfer process is neglected. 

Let us calculate the velocity of the propagation of the primary shock 
Wave Mg and consider some properties of the solution of the stationary 
progressive wave type [7]. 


As we assumed that the density and pressure before the explosion on the 
right and left from the membrane are equal, the Mach number of the shock wave 


is found from the following equation i 


^0 (p 1) L(p 0[rfr+l)g * ? 1 


Id. 


q.Q/a] 




’ 1^0 


( 2 ) 
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o 


Dg - the velocity of the initial shock wave| 

a^ - the velocity of sound in front of the initial shock wavoo 

The relations M 0 obtained by solving this equationi are 

drawn in Fig. 2. 




Let us now calculate the velocity of the detonation wave, which is the 
solution of equation (l). Let a plane shock wave with Mach number Mj 
propagate in the fuel gas. After the transformation of equations (l) to a 
coordinate system moving with the shock wave, for a stationary one-dimensional 
flow behind the shock, the following equation is valid:- 


dX 



This equation combines the velocity v* relative to the shocl^ with the concentration 
of the product 'X . By integrating this equation and satisfying the conditions 
for the shock, the following relation is received: 


(7* 





From this equ^ion, it follows that in the case of a complete heat release, a 
steady state solution behind the shock wave may be present only for Mach number 
M^ N My, where M^i is the Mach number in Ch*-J. system, and is determined by 
the formula: ^ 
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Thus, the distribution of the parameters inside the detonation wave 
depends only on the specific form of kinetics, whereas the velocity of 
propagation of the stationary wave of Ch»-J» depends only on -y and Q, 

In particular, the velocity of propagation of the stationary double-front 
wave in the Chapman-Jouguet system at given values of y and Q coincides with 
the velocity of the wave of Gh.-Jo with a distributed heat release. The 
stationary progressive wave with any number of reactions has a similar propertyo 

The relation ■My(YfQ)is drawn in Fig, 2 together vfith the 
relation (ry , Q) , 

2, To obtain a quantitative information about the formation of detonation 
and about the gas movement after sufficiently big intervals of time, it 
is necessary to integrate numerically system (l) at the corresponding initial 
and boundary conditions. The one-dimensional equations reduced to the 
dimensionless form, in Lagrangian coordinates, are written in the form:- 


dx 

at 


= a 


3u dP ’ fj 
at 5x 




axp 

dx 


d$ - „ ddxJdXa 
at at 


H ''pj P, pi' 'P>i P 


, 


Pdx 






Initial conditions:; 


L^x,<oo CiQ^x,)-P{O,x^)-0 

O.;cp<oo u(0,x,)-0 J( 0 ,x.) = jpr) . ^'"^0 • 


The boundary conditions j u(t, O) ■ 0. 
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The dimensionleaa quantities are expressed in dimensional quantities 
(maked by dashw^i^^ the form: 

, 7T73r, ^ ' 

P-P'/fP, , S-S7a,‘ , 9‘Q'/a; , , 

L =L'fi/a, , r ppP'B . 

Here, a., J’j } ?x i - are the velocity of sound, density, pressure and 
temperature in the undisturbed gas on the right, respectively. 


m~1 

B^aP^ exp (-^J 

The quantities behind the primary shock wave are racked by index & . The 
velocity of this wave is calculated by formula (2) for Oo the left, equal 
to on the right. The relation Qo/Qi* K is ilso a characteristic parameter 
of 'the problem. Thus, the formulated one-dimensional nonstat-idnary problem 
contains the following dimensionless determining parameters: 


][ , m ,d. , E ^ Q , K , L , r , 

The one-dimensional variant of the conservative obvious numerical scheme 
of the second order of accuracy with direct calculation of the discontinuities 
was used in the calculations after it was fitted to problem with distributed 
heat supply. This one-dimensional variant vms worked out in the works [9 and 10] 
for the solution of the elastic-plastic problems and problems of gas dynamics. 

The untaown quantities £ira calculated by the formulae: 
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The integral lower indices are given to the boundaries of the particlesi 
the half-integral— to the particles. The upper indices designate the numbers 
of the temporary layers. The quantities x and u are determined for the 
bovindaries of the particles} the quantities 6 t C, p - inside the particles. 

The velocity is calculated at half-integral instants. 

The interval of time was selected from the requirements of the stability 
of calculation. 

3. The calculations were carried out at E • 10} Q ■ 24, F»0.3,ra-1, 

L - 3 for 0 and l and for K ■ 0.8, K ■ 1.0, K » 1.5. The value of 
y was taken equal to 1,4. From 100 to 150 particles were placed in the section 
with length L. The coefficient of the artificial viscosity d-* was equal to 
five. All the particles in the calculating range were originally the same. 

The slightly smoothed profiles of the velocities for different instants in 
the case of 1, K - 1 wore given in Fig, 3<> At t - 0.146,- the profile 
consists of a compression shock on the right, linear part in the rarefaction wave 
on the left, and constant velocity between the shock and the rarefaction wave. 

At t - 0,767 this simple distribution is disarranged because heat release started 
on the contact discontinuity (marked by circles) and it led to the generation of 
an internal shock wave (the mark * indicates where the time of induction elapsed) 
in the gas behind the initial shock wave, whereas in the strongly heated products 
of combustion, where the density is small and the velocity of sound is high, a 
compression wave propagates. As a result of the compression wave, which is 
propagated in the gas behind the initial shock wave until the second shock wave 
arises, the velocity level between the shock waves is higher than the value 
Ug a 1,77 obtained by the exact formulae of the decay of the discontinuity. 

On the velocity profile at t - 1.236, two shocks of velocity corresponding to 
two shock waves are recorded. The second shock wave is the leading front of 
the formed detonation wave and behind it there is a heat release zone. The 
particle, in which 0.5 Q was separated out is marked by a pointer, and the 
particle, in which 0.9 Q. was separated out — a pointer with two lines. The 
heat supply between the contact surface and the wave front leads to the 
intensification of the compression wave going from the zone of the heat supply 
to the products of the explosion. Therefore, the velocity of the gas flowing 
from the zone of explosion decreases and a gap appears on the velocity profile. 

The curve of t ■ 1*439 represents the distribution of the velocities after the 
occurrence of interaction between the initial shock wave and the shock wave with 
the following zone of heat supply behind it, and the shock wave with the zone of 
release behind it went out in the stationary gas. We notice that the velocity 
behind the shock exceeds the velocity behind the shock front of the wave of Ch.-Jo, 
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. and the products of the explosion by then slowed down still more strongly© 

At t = 1.762 and t » 2.048f the velocity behind the shock is practically 
equal to the velocity behind the front of the wave of Ch,-J,, and a portion 
of the products of detonation and explosion has a negative, that is to say, 
the gas flows in blie csfposite direction© The steepness of the parts of 
negative slope in the profiles of the velocity increases with increasing the 
time. Taking into accoimt the relatively high temperature of the products 
of the explosion and the established artificial viscosity in the scheme, it 
is possible to consider that the retonation shock wave was formed at the 
moment t - 2.048, 


C 
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A fragment of the graph in Fig, 3j where all the boundaries of the units 
are drawn, is represented without smoothing in Fig, 4 , This fragment gives a 
concept about the magnitude and character of the oscillations of numerical nature 
for the investigated wave propagated in the Ch,*-3’* system. 





The x-t diagrams for the shock waves and the leading edges of the detonation 
waves (points or points connected by lines), of the contact surface 
(triangles) and of the trajectory r (t) of the retonation wave (rhombs) are 
shown in Figs. 5“9» The dependence of the pressure and velocity behind the 
front of the shock and detonation waves on time, Po(t) and Up(t), are also 
drawn in these figures and are marked by crosses and circles, respectively. 

In the figures, there are also the reference V£ilues of the velocity Uy and 
pressure behind the detonation wave in Ch,-J, system, and the velocity 
U 9 and pressure Pq behind the initial shock wave, which are obtained from 
the exact relation. The initial slope of the trajectory of the contact surface 
is drawn by dot~dash line; the exact slope of the trajectory of the initial 
shock wave - by double dots-dash line; the slope of the trajectory of the 
detonation wave in the stationary gas in the Ch.-J, system— solid line. 

Let us return to Fig. 5« The aupper succession of the black points represents 
the trajectory of the initial shock wave. On the lower series of the points and 
triangles, the finned pointer shows the place, where for the first time on the 
contact surface, the time of induction elapsed and the combustion started. Since 
that moment, the contact surface begins to lag behind the trajectory of the 
ignition, on vfhich the time of induction elapses; however, the shock wave with 
the zone of heat release behind it, does not still arise for some time. The 
formation and acceleration of the shock wave with heat supply begins there, 
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where the graphs of the pressures and velocities behind the initial wave and 
on the trajectory of the ignition bifurcate. Till that instant, the greater 
part of the yielded energy is carried away from the zone of heat supply by 
the compression waves going in the direction of the light products of the 
explosion, which are stopped at that time. After the generation of the internal 
shock wave with a sufficient intensity in the compressed gas behind the initial 
shock wave, the time of induction of the particles quickly elapses on passing 
through this internal wave and the movement of that wave gets a self-accelerating 
character. The trajectory of the shock wave with heat release behind it is bent 
upwards and intersects with the trajectory of the initial wave? after that only 
the shock with heat release behind it (detonation wave) remains. Here, the graphs 


C 
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of u and p behind the initial shock are cut off. In the graphs for the 
shock wave with heat supply, it is obvious that the pressures and velocities 
are greater behind the front of that wave than behind the- front of the wave 
of Chapman-Jouguet in the stationary gas. However, they rapidly decrease to 
the values corresponding to the parameters behind the front of the wave of 
Ch.-J. as soon as the detonation wave goes out in the stationary gas. 


U,*S.S6 



Fig. 6. 

The trajectory of the shock and detonation waves, of the contact surface 
and retonation wave, and also the graphs (t) and Uo(t) for the case of 
- 0, K « 1 are drawn in Fig. 6. It is obvious from the given data that 
although the internal shock with zone of heat supply behind it is generated 
at £he same time as in the case of » 1, K « 1, it is accelerated afterwards 
more slowly than in the case of «><»» 1, K ■ 1, and at the fusion with the 
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initial shock, the developed detonation wave does not form. In the graphs Po(^) 
and Uo(t) , there are no bursts of pressure and velocity, and subsequently a 
relatively slow acceleration of the shock wave with zone of heat supply behind it 
is observed. In the investigated section of time, the velocity of the shock and 
the velocity and pressure behind it are still far from the corresponding values 
for the wave in the Ch,~J. system o 

We notice that even in that case, when the resulting wave is relatively slow, 
the retonation compression wave did not succeed yet to reflect from the wall, 
and thus it is possible to neglect the influence of the wall ^n the flow in the 
neighbourhood of the shock with heat release behind it. 
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The x-t diagrams and the graphs Pp(t) and U^Ct) for the case of -1, 

K = 0,8 are given in Fig, As was expected, the beginning of the heat release 
on the contact surface (finned pointer) is protracted in comparison with the 
case of « 1, K = 1 , The beginning of the generation of the internal shock 
with zone of heat release behind it is accordingly drawn back. In consequence, the 
generated internal wave with heat supply behind it may be for a sufficiently 
long time accelerated in the gas behind the initial shock wave before it Hows 
with ito At that time, a detonation wave generates with a leading front that is 
more intense than that in the case of <=< » 1, K « 1, It is obvious from the 
comparison of Figs* 5 and 7 that the waves of pressure and velocity generated 
behind the front exceed the corresponding values for the case of ‘^■l,K-l, 
As in the case of ■ 1, K « 1,. after the exit in the stationaiy gas, the 
velocity of the detonation wave rapidly decreases down to the velocity of the 
wave of Ch»-Jo, and the velocities and pressures behind the front of the wave 
approach to the corresponding values behind the front of the wave of Ch,-J. 

The X- t diagrams and the graphs of Po(t) and Uo(t) fortf< ■ 1, K ■ 1,5 
are given in Fig, 8, In this case, the isaitial shock wave is the most intense 
and the internal shock with heat supply generates more quickly than at K ■ 0,8 
and K » 1, However, its intensity by the moment of the fusion with the initial 
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shock wave, as in the case of ■ 0, K ■ 1, is insignificant, and the 
detonation wave goes out in the undisturbed gas with a leading front 
is more slow than that in the system of Chapman-Jouguet. Subsequently, the 
velocity of the detonation wave of ■ 1, K ■ l*5j increases more rapidly 
than that for the case of 0, K « 1# We notice also a different form 
of the trajectories of the contact discontinuities for the two cases. 



4. The relief of the velocity in x, t coordinates for the case of =1, 

K « 0.8 is given in Fig, 9» It is clear how the internal detonation wave ojiises, 
and how also the retonation compression wave, from which the retonation shock 
wave is subsequently formed, arises. Finally, it is clear how the detonation 
wave transits to the system of Ch.-J, after going out in the stationary gas. 
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This pattern may serve as one of the illustrations of the process of the 
formation of detonation on the ignition of the gas at the closed end of the 
tube and as an example of the autonomous flow which arises in the 
unconstrained gas in the absence of external influences owing to the released 
energy in the chemical reaction which is characteristic for the movement 
occurring in the self-sustaining detonation waves© 
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EXPLOSION IN FUEL MIXTURE OF GASES 
By 

VoP* Korobe^iikovi V«A. Levin, V.V. Markov 


Let in an infinite mass of stationary gas, in which chemical reactions 
may proceed, instantaneous energy release E^ occurs at a point along a plane or 
straight line. Let us consider that the gas is ideal, nonviscous and non-heat- 
conducting. The density' J% and pressure Po in the stationary gas are constants. 
One-dimensional motions of the gas with different form’s of symmetry are 
investigated. 

The strong shock wave generated as a result of the energy release initiates 
the chemical reactions in the fuel mixture of gases. On calculating the 
combustion process, a model is assumed, which takes into account the time 
delay of ignition and the subsequent simultaneous proceeding of the direct and 
reverse reactions. 


The equations describing the proceeding of the chemical reactions are taken 
in the form of Arrhenius relations. Tiie reaction determining the period of 
induction is described by the equation in [l, 8, 9], 


qf C 


(1) 


and the reaction with heat release - by the equation in reference [2]j 

( 2 ) 




dfi 


n\. 


it - 

P f e 


RT 


where C - virtual concentration, p - mass portion of the unburned gas, 

Q - calorific power of the unit mass of the fuel mixture, E - activation j. 
energy of the induction period, Ej^ - activation energy of the direct reaction, 

R - universal gas constant, K > 0, k^ > 0, k > 0, m, , n, , , m^, , -^4 

—are constants. The quantity C « 1 is on the front of the shock wave. /The 
vanishing of C denotes the end of the period of induction and beginning the 
reaction with heat release. . / 

The motion of the fUel gas generated as a result of the explosion, will be 
described by equations (l) and (2) together with the equations of motiom 



*(z*'*pv),^0 


( 3 ) 
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w) 

(?) 


Here j-i f t P - , v-1,2, 3 -'for motion with planci 

cylindrical and spherical waves, respectively. The parameters of the gas 
must also satisfy the boundary conditions in the center of the explosion and 
on the shock wave. The velocity vanishes in the center of the explosion, and 
the following relations are fulfilled for the shock waver 


p J}\ ~=-(d~v )'% * 

r-' Po 2 ‘ ' r-f 


[ 6 ) .. 


where D - the velocity of the shock wave, and the indices 0 and 1 are related to 
the values of the parameters before and after the shock wave, respectively. 

The system of equations (l) - (5) is linear, and the problem in the given 
statement is not self-similar. Therefore, it is possible to carry out the complete 
investigation only by the method of numerical integration of the equations in 
partial derivatives. 


Similar to work [2], we shall study the motion of the gas at an instant near 
to the initial motion. 

It is possible to obtain an approximate analytical formula for the time 
of induction. In order to do thatj we write the thermodynamic functions 
in adiabatically expanding particle behind the shock wave in the form given 
in work [4]j 




(7) 


where the index * - is related to the value of the function at the moment of 
passing by the particle of shock wave, which is marked by t^ , a and 6 - 
are constants. 
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An analytical solution of the problem about a strong point explosion was 
used for the determination of a and &. The values of a and 6 were determined 
from the condition of equality of the first derivatives with respect to V, 
from the exact solution and formulae (7)» *^he values of a and B for certain 
U and 7 * are given in the tables 


V 

1 

2 

3 

r\ 

a 

8 

a 

8 

a ■ 

8 

1.1 

0,085 

0,476 

0,141 

0,703 

0,168 

0,838 

1.2 

0,163 

0,455 

0,264 

0,661 

0,314 

0,785 • 

1,3 

0,261 

0,435 

0,374 

0,624 

0,442 

0,750 

1.4 

0,333 

0,417 

0,472 

0,580 

0,556 

0,684 

3,0 

1,000 

0,250 

1,250 

0,313 

1,400 

0,350 


Substituting formula (?) in equation (l) and integrating it, we obtains 


f »/- 




['■(f) J ■ 


( 8 ) 


where 


A = 28 fifn*Q 




/ n 

V- p* 


If to is the instant at which C vanishes, the difference t© “■ t* gives the 
» value of the period of induction T, as a function of t, . From (8) we 
obtains; 


' -l] 




(9) 


Previously, for example in work [7]| the following formula taken from work [l] 
was used for the induction times 


f; 


ind 


Betp(C/RTl 


ClO) 
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Below, for the investigation of the initial stage of the explosion, the next 
formula is used, following work [3]j 


, ( 11 ) 

where f ' are constants. Both these formulae give the value of the 
period of induction on combustion of stoichiometric mixture of hydrogen with 
oxygen, initiated by a shock wave. 

Fig, 1 is presented for the comparison of the graph of the dimensionless 
values of 'X^ ), calculated by formula (9) (curve l) with those obtained 

by integrating equation (l), taking account of (ll) along the streamline and 
using the exact solution of the problem about a strong point explosion 
(curve 2), for the case of Y - 1.3 and ^ - 3« It is obvious that the 

approximate formula is in good agreement With the calculation. 



Fig, 1 


The performed analysis showed that already a;t the instant close to 
the initial, the time of the ignition delay abruptly increases, which leads to 
the separation of the ignition zone from the shock wave. 





Fig, CL Fig, ^ 


The results of the calculations showing the influence of the different 
determining parameters of the problem on the withdrawal of the ignition front 
from the shock wave front are represented in figs, 2~l^, The value of the time, 
at which the shock wave passes through the particle, relative to the characteristic 
time of induction is plotted on the abscissa. The dimensionless distance from the 
shock wave, at which the period of induction is completed in that particle, 
is plotted on the ordinate. 

Fig, 2 illustrates the influence of the adiabatic ‘index y* on the decay 
of the detonation front for the case of ■ 3j ■ 10 erg, (E/q) m cr' m 10, 

The dependence of this process on the form of the symmetry for y "1,31 
Eo “10 erg/cm^"^ and •• 10 is represented in Fig, 3* The influence 
of the magnitude of the energy E of explosion and the activation energy E of 
the induction period when "f » 1,3 and V m ^ ±b represented ^in Fig, The 

curves 1, 2, 3 are for (E^ » 10 erg, cr' = 13)i (E o • 10 erg, o-' » 10), 
and (Eq s 10^^ erg, o lO), respectively. 

Two curves are shown for comparison in Fig, 5* ^e, represented by solid 
line, is obtained in the present work by the formula for the time of induction (ll), 
and the other, dotted, — in work [5]| on the basis of formula (lO) , In the 
initial stage of the explosion, it is possible to obtain a distribution of all 
the unknown functions on the cooi^inate at certain instants. Such distribution 
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Is necessary for the continuation of the solution at still later lnstantS| 
at which it is already impossible to neglect the influence of the energy 
released by the reaction and of the pressure in the undisturbed gas* The 
calcvilation of the initial distribution shows that the value of C monoto- 
rically decreases from 1 to 0* 
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The value of /B rapidly decreases at once after the ignition fronts 
and then it increases to the equilibrium value in the center of the explosioh^ 
which is calculated by the formulai 


5 - 

Kf ' ‘ 

The mode of the curve of ia represented in Fig. 6. 
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The calculation of the later stages of the explosloni taking account of 
the energy released by the chemical reaction and the conterpressure, was carried 
out for the case of the cylindrical symmetry using the differential method of 
S,K« Godunov [6]* Direct application of this method near the center of the 
explosion is difficult because of the high temperatures. Therefore we are 
obliged here to carry out the calculation by another method. The central 
interval, representing the standard calculating \mit, is introduced in the 
investigation; one of the boundaries of this unit coincides with the center 
of the explosion. On this boundary, besides the condition of the equality of 
gas velocity to zero, some other conditions are set out. The value of is 
supposed to be equal to the value calculated by (12), density- to zero, 
pressure— to the pressureffirom the lower temporary layer. The validity oftthe 
first two conditions is secured by the fact that the equality of the density 
to zero and the infinity of the temperature in the center of the explosion 
remain at all the time of the gas motion, because this takes place in the case 
of ordinary explosion with counterpressure in view of the conservation of 
entropy in the zona of continuous flow [?]• In our case, because of the 
chemical reactions, the entropy may only increase (it is assumed that the 
reactions as a whole are exothermic). The dimensions of the neighbourhood of 
the center of the explosion were automaticaly selected so that a small amount 
of mass and energy was included in it. The ignition front was included in 
the boundary of the calculated unit. Its position was determined so that the 
values of the modulus C extrapolated to it fi*om the neighbouring units was 
smaller than the small positive nwber ^ ._Enuations (1-5) were solved for 
the dimensionless functions j P= ? . a.nc| the ' ' 

variables: t =^ ^ ‘ The parameters,^by wM.ch the numerical 

calculation were carried out ijy the differential method, were as follows: 


_ 10 

■» 10 erg/cm 

fO , , 

Q “ 7 . 10 cm/sec 


- lO'^’^/cm^ 

6 .2 

“ 10 dW/fciii “ 1 atm 

t*n _/» 3 


a^--Ejq-Z.5 
i - - 0 

nil ^ - 2 

Pi* n^~ 2 n-1 


10 



121 - 




Control of calculation was carried out by the laws of conservation of 
mass and energy. The values of the relative errors ^ and £ ^ were 
calculated: f . 


/ -£il£o 
f f o ’ 

1 o ^ T'^ f ' — the total energy 

of the gas being in motion after the deduction of the internal energy of the 
original state. 


f^__the dimensionless energy of the explosion. 






where 



— the mass of the disturbed gaS| 


the mass of the undisturbed gas in the volume limited by the shock wave. 

The calculation was carried_out from t • 10 to t - 500. It was discovered 
that begirming from the Instant t « 19» the ignition front merges with the 
streamline and then it follows together with it. Oscillations of the ignition 
front were observed beginning from the instant t ■ 50. Some residts of the 
calculation are shown in Figs. 6-7. It is possible to judge about the calculation 
accuracy of the gas dynamical functions from the following data. The initial 
error (t ■ 10) r ~ The error at the instant t ■ 19| 

(fp - - 3*3^J ■ - 0 . 069 s. For still later instantsi the error did not 

increase. 
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Thus the motion of the gaS| in which exothermic chemical reactions, 
generated as a result of the instantaneous release of energy E© in a point, 
can proceed on a straight line or a plane can be divide into two stages. 

The first stage is that where the quantity of the energy released in the 
region bounded by the front of the flame is small compared to the energy E© 
of the explosion. In this stage, the flow is described by the formulae giving 
the solution of the problem of a strong point e^qslosion. The chemical reactions 
occur at the background of this flow. At the seme time, in the case of 
explosion, the generated supercompresaed wave of detonation decays on the 
ordinary shock wave and i^tion front. The lncrease,with time, of the distance 
between the shock wave and the ignition flronts essentially depends on the form 
of the syimuetry, the energy Eo of the explosion, the activation energy E for 
the induction period, and the adiabatic index y . . 

The second stage differs by the fact that it is necessary to take, 
into account the energy released by the chemical reaction. The numerical 
investigation of the flow in this stage for the values" of the parameters 
mentioned above (l3) showed that in spite of the energy supply, the ignition 
Ij’oiit continues to lag behind the shock wave and flows together with the 
trajectory of a determined particle of the gas, i.e. in the present case, the 
detonation combustion is not regenerated. 
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ONE-DIMEMSIONAL NONSTATIONARY MOTIONS OF A FUEL MIXTURE OF 
GASES IN CASE OF SMALL HEAT EFFECTS OF THE CHEMICAL REACTIONS 

• By ■ 

Ii.I* Zak and V«A> Levin 


The one-dlnenslonal unstatlonary motions of a fUel mixture of gasesi 
accompanied by the presence of shock wavesy and n involving exothermic chemical 
reaction mechanisms are investigated. As a result of the proceeding of the 
chemical reactions in the gas flow behind the shock wavOf heat is released. 

This heat affects the motion of the shock wave and the gas flow behind it. The 
quantity of the heat energy reserved in a unit mass of the iUel gas mixture is 
supposed to be a sufficiently small value. In an analogous statement t the 
propagation of a shock wave in an Ideally- dissociable gas was investigated in 
work [5]j in which the problem of a piston moving with a constant velocity was 
solved. 


Movement of a Piston in a Fuel Mixture of Gases 


Let at the initial Instant a piston starts to move in a stationary gaseous 
fuel medium, A shock wave is fonned in front of the piston. This wave initiates 
the chemical reactions occurring with a release of heat. We shall consider the 
case when the heat effect of the chemical reactions exerts a small influence on 
the gas motion. This assumption will be achieved in the case of impoverished fuel 
mixtures. If Q denotes the quantity of the heat energy released by complete 
combustion of a xinit mass of the fuel, then, in the present case, the following 
relation must be fulfilledt: 


fli 





O 

Hero - the density of the burning gas, •f' - the total density, a^ - soxuxi 
velocity behind the shock wave. The superscript "0" corresponds to the parameters 
of the gas in front of the shock wave. 


It is necessary to mention that the unequality (l) will be also fulfilled 
in the case of piston movement with a sufficiently high velocity and for 
high-grade fuel mixtures. 
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Let U8 vfrite the eyetem of the equations describing the motion of the gas 
in the formj 


3p_ 

at 


3u 

*u — — =0 


dX 







dxJ 









Here h| P| / - the enthalpyi pressure and velocity of the medium, respectively, 
yQ- the mass concentration of the uiureacted gas, Q - the total heat release in 
the unit mass of the gas and assumed to be of a small value in the abovennentioned 
sense;; 't - the adiabatic index* 


Let us write the equation simulating the proceeding of the chemical reactions 
in the formj 


dj^ 

d t 




(3) 


For the investigation of the gas motion, it is convenient to convert from the 
Euler variables jc, t to the Lagrange variables t, T" , where T- the instant at 
which the given particle is set in motion when the shock wave passes through it* 
With these variables, the system (2,3) has the form*: 


3x 


=" . f S --fuM 


^ '■ 'Bt at ' at f ■at ^ at 


^4 . 


at 


(4) 

(5) 


Here D(t)— the velocity of the shock wave front* The parameters of the flow 
must satisfy the conditions on the shock wave at t ■ t 


2D/7'o‘’V~ 


IlL 

J'*n 2 )“/ ’ •* 

fff - -•* 


* — r P ^ 




j 


-f 


(.O 


and the condition on the piston at X ■ 0 
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Let the velocity of the movement of the piston slightly differs from the 
constant value Ug | i«e*| 

In connection vdth this, we shall seek the solution of the system (4i 5) in the 

formj ; ^ z 7 — z 

u--Uj}*6u,]^ p-.pJUeP,) ^ 

. ■ (7) 


The values with the subscript ”0” correspond to the common gas-dynamical flow 
from a piston moving with a constant velocity! and are determined by the 


formulae: .... . 

) . p -..cLp\J-p‘s , 


. f*' /, 2 fL V’ 

■ C.2) 

Substituting formula (?) in equati<m (A) and using the assumption - f 

we obtain for the values u, Pf f . 

r V 1 


i , Po^odt ' 3i dz ^ ’ 





At the same time the chemical reactions will proceed in the' given field of the 
gee-dynamical parameters! andi instead of equation (5) i we have the equationi; 


3t 




(1.0) 


For ttie parameters of the gas on the shook wave at t "C i we obtain! 
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p - - 


< ri 




On the piston at TT • 0, the foUovdng condition ncill be fulfilled 


u,(t,o)^V,(t) . 

The equation deacribing the proceeding of the chemical reaction la integrated 
irrespective of system (9)f and its solution, taking account of the boundary 
condition on the shock wave| vdll be ■ F(t-Tr), ' The specific form of the 
function F(t-r) depends on the constants ro^ and entering in .equation (5)* 

Excluding the parameters and jo in equations (9)| we obtain the 
follovdng equation for the quantity Uj 

, ,, 

3t’~ fpi)M;*2 Un‘-lj JtJr • (12) 

The general solution of this equation will be» 





where 


cj^- 








For pressure Pi we obtain 




Using the boundary conditions on the shock wave, (ll) we find thatt 


Hi)' 






(j-1 

IZ 




i> ^ , Cjtl 

- 1 ) 2 <^ ^ 


(15) 


C 
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Using the relations (l5)f it is easy to obtain the relation between the functions 

and t 

' ^ . 2A1e • (l6) 

This relation shows how the disturbance of the pressure changes when it is ^ « 

reflected from the shock wave. The quantity .^uC^^o+vJ 

is called the reElection coefficient of small dlsturbancasf and it concides 
with the corresponding coefficient on reflection of disturbances from the shock 
wave in the common gas dynamics [1-3] • 

Substituting formulae (15) in formulae (l3) and ( 14)1 wo express the 
increment of the gas velocity and pressure by the increment of the shpck wave 
velocity:: 



/j 1 *c-i/'aj + /) \ * t) „ fut *t \ c .1 


fot*r\ r / , 7 ? 


(17) 


For finding the function D , ( C>)| we use the condition for the piston:: 
u(t,0) ■ (t). From thia-condition, we find the functional equation which 

is satisfied by the increment of shock waye velocity:- 




( 18 ) 


Thus, the problem of a piston movement with a nearly constant velocity in the 
fuel mixture of gases is reduced to the solution of a functional equation. 

The solution of this equation^ as shown in f 2 T written in the 

form:: l- J . 





( 18 ) 


C 
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The solution consists of two absolutely convergent aeries. The first sum 
ensures the change of the velocity of the shock wave owing to the heat release 
in the gas flow, and the second— owing to the change of the velocity of 
movement of the piston (as in common gas dynamics) • Let us consider in more 
detail that component of the increment of the velocity of the shock wave which 
is caused by the proceeding of the chemical reaction. For that, we put 
Vi(^)-Oin ( 19 ) • In t at case, it is easy to obtain the limiting increment 
of the velocity from ( 18 ) , with transition to the limit at ^ 00 

/ 

, , Zu{ni-t) 

( 20 ) 


Where F («=* ) ■ ^ — the equilibrixim value of the concentration which is the 
solution of the algebraic equation 1 


«>(?. .s-JO-aP-kM - ( 21 ) 


Thus, the limiting increment of the velocity of the shock wave is determined 
by the equilbrium value of the concentration of the component participating in 
the reaction. The knowledge of the function F (Ei) is necessary for the 
detailed distribution of the gas-dynamical parameters with respect to the 
coordinate and time. For example, at ra^ - m^j ■ 1, the solution of equation (lO) 
has the forrai 


and thereby the form of the function F ( ^ ) is detarmlned. For that case, 
the limiting value of the increment of the shock wave velocity as a function 
of the Mach number of the shock wave of the main flow Me, at « 1.3 is 
represented in Fig. 1. _The change of the increment of the -shock wave 
velocity with the time t ■ (kj + k^^) t for some values of Mo is shown 
in Fig. 2 . It is obvious from the ''graph that the velocity monotomically 
tends to its limiting value. The change of the pressure on the piston with 
time is represented in Fig. 3* That change is nontnotoulc. There is an 
obviously expressed maximum of pressure, which was also discovered on the 
body, on calculating the supersonic flow of a fuel mixture of gases around 
a cone [ 4 ]. • 

Let a shock tube in the form of a canal of a constant cross-section be 
filled with two stationary gases separated by a diaphragm. On the left fl*om 
the diaphragm, there is an inert gas with a high pressure, the so-called 
"driver gas", and the right part of the tube is filled with a fuel mixture 
of gases with a high pressure. After the rupture of the diaphragm, a shock 
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Wave will propagate in the fuel mixture and ignite it* Heat will be released 
in the gas flow behind the shock wavei and the motion of the shock wave as well 
as of the contact surface becomes nonuniform* A Rieraannian central wave of 
rarefaction propagates with high pressure in the gas* In contrast to the 
classical problem of the decay of the arbitrary discontinuity, this wave will 
not be self-similar* 



Movement of a fuel mixture of gases In a shock tube* 


C 
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The flow in the rarefaction wave is described by the formulae* 




a * or. 


(23) 


Here V' - tunknown function determined by the continuHy conditions of the 
velocity and pressure on the contact surface* The subscript ^ corresponds to the 
state of the stationary gas with high pressure* 


Let us use the general solution (l7)f obtained previously, for the description 
of the movement of the fuel mixture of gases* The following condition must bo 
fulfilled on the contact surface* 


2f, 

c fi L± 

2 ■' 


= PoO*^Pt(t)) , 


(24) 




Reforming condition (24) and expanding in a series through orders of the small 
parameter £ , we find the relation between the increment of the pressxxre and 
the increment of the velocity of the contact discontinuity in the form* 






(26) 


Here the value of Mo is a function of the decrease of the pressure and 
temperature and it is determined by the relation* 


p. r 


' f . 


hll 





Ll 

r^f- • 


(27) 


Substituting the corresponding formulae for (t) and u^. (t) in the condition 
(26), we obtain the functional equation for the unknown function Dj (^)* 
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where :: 


S2 -cj 


\--- 


^r'^: - n 


J_ ] a" 

,i2= 


SI* 


h = 


^j^o-fr’) 

¥^o-0 

^r'^c-fr-}) 




SI *cj 


.1 \ 

r*< ” /^J 


( 28 ) 


Using the equations (l3) and (14) i we find the relation between the Ihnction 
f^ and f;^ on the contact discontinxxityi 


(29) 


Fron this relationy the physical meaning of the quantity Aj becomes clear* 

It is the coefficient of reflection of the compression wave going backward from 
the shock wavey from the contact surface of the wave* 

The solution of the ^ctional eqaution (28) will bej ‘ 


F’,((h 


2cJA (^o ' 

2Ml*uff^l*l) 


l(nS[Ff{^‘‘"^)]' . (30) 

rt»fl ^ 


For the limiting increment of the velocity of the shock wavey we obtains 



Zh\*u(t\l*l) 


Lh. 



(31) 
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The final velocity of the shock HavOf both in this case and in the case of the 
movement of a pistoOf is determined by the equilibrium concentration of the 
component participating in the reaction. The change of the velocity of the 
shock wave with time for the values ■ 1.66} y ■ 1.3 and for some Mach 

numbers Mo at the temperature drop m is represented in Fig, The 

velocity of the shock wave monotonically tends to its limiting value. The 
change of the velocity of the contact discontinxxity and the pressure on it with 
time are shown in Figs. 5 and 6* The velocity of the contact discontinuity 
nonmonotonically decreases. It has a velocity minimum» The character of the 
change of. the preas\ire is the same as on the piston, i.e* the presstxre nonraonoto- 
nically increases, tending to its limiting value. 


C 






fig. 3 
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The knowledge of the increment of the shock wave velocity in the experiment 
enables to determine some functional relations between the constants determining 
the proceeding of the chemical reactions* ThuS| in the case of a single reaction 
of type ( 3)1 it is possible to determine the affinity constant from the relation 
( 31 ) • One should also mention that the discussed method of the construction of 
solution can be generalized for the cases of many chemical reactions proceeding 
both simultaneously and successively. 
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HYPERSONIC NONSTATIONARY FLOW OF A FUEL MIXTURE OF 
GASES IN THE NEIGHBOURHOOD OF THE CRITICAL LINE OF 
A BLUNTED BODY 

By 

SJ4. Gilinskii 


On shooting a body vd.th high velocitias in hydrogen-air and hydrogen-oxygen 
mixtxiresi soma characteristic modes of flow around were ob served r stationary— with 
plain shock and thermal fronts; nonstationary— of pulsating character with a 
strictly periodic structure of the combustion region in the form of circular waves; 
modes with a complex sti*ucture, with the formation of internal exploisons inter- 
acting with the leading shock wave and bending it# The theoretical investigation 
of the stationary modes was carried out previously^ and the main results are 
discussed in fls'J . 

In the present work| the pulsating nonstationary modes of combustion 
generated at hypersonic flow around the leading part of blunted bodies are 
investigated. The investigation is carried out numerically. The flow is 
studied in the neighbourhood of the critical line. The method of "termination 
of series" is used and the solution is searched for in the form of a first 
approximation. The integration of the two^iimonsional equations in partial 
derivatives is carried out by the finite difference method with the use of the 
characteristic relations for discontinuities and for the contour of the body* 

The initial stationary flow is determined from the solution of the boundary- 
value problem for the system of ordinary differential equations by the method of 
iteration. 

The growth of finite disturbances artificially Inserted in the boundary and 
initial conditionsj and also of the small disturbances generated because of 
nximarical errors is studied with time. It is ascertained that the inserted 
disturbances lead to the formation of the oscillations of the loading shock 
wave and the parameters of the gas behind it. The amplitude and frequency of 
oscillations depend on the magnitude of the activation energy. 

The detected modes are near to self-oscillating ones with constant 
amplitude and frequency. The region of stability, with respect to the small 
and finite disturbance, for the two-dimensional flow is significantly larger 
than that in the case of the propagation of one-dimensional detonation. 

§ i. Let a supersonic flow of a fhel mixture of gases with the parameters! 
velocity , prossxxre Poo and temperature T«> t runs against an axisymraetric 
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blunted body* We eball assume that ignition does not occur in the incident 
floW| and that an exothermic irreversible reaction A-*>B proceeds behind the 
shock wave because of the increase of temperature and pressure according to 
Arrhenius law: ^ 


dS 

¥ 

di 



e 


■i/nr 


= a>/'/7,Ty . 


(l.l) 


Here, as usualy ^ the concentration of the unreacted gas. O'* order of reaction, 
activation energy, L— reaction rate constant, gas constant. 

The system of gas dynamic equations together with the kinetic one (l*l) 
for the axisymmetric flow can be written in Eulerian variables in the following 
form:: 


■ « 
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du V du 
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= -Lfi'^p"''^expf-£ji/p} , 
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where Z ? ® — the polar coordinates with the origin in the center of the 
curvature of the body in the critical point at zero tlmej u, r- the projections 
of the velocity vector iT in the direction Z and & as shown in Pig* 1, 
h, T - the pressure, density, enthalpy and tomporatiure, respectively^ 

Y •“ adiabatic index, Q - energy release of the reaction* 
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We Interpolate the dimenaionlesa variablsa (with index l), with the help 
of the following relations; 


^ ^ ‘"7^ ^ . 
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' L 

To: 


(1.3) 


Here the subscript « « refers to the parameter of the Incident flow, and the 
superscript •’0" denotes that the value of the parameter is considered at zero 
time, Vmax the maximum velocity of the gas in front of the shock wave. 


The fora of equations (1*2) remains in the diraensioxiless variables (below 
we shall onut the index ’•!")• The boundary conditions on the surface of the 
body are given - conditions of non-flowing and in the incident flow for , 
Poo , Too and x • The shock wave is considered a discontinuity, where 
Hugoniot relations are fulfilled* 


Thus, the solution of the problem in the dimensloiiLess variables is 
determined by five main parameters M^o t y » ^ t ^ parameters 

characterizing the disturbance* The heat release q of the reaction at a fixed 
Mach number Moo characterizes the degree of the supercompreasion of the 
stationary detonation wave* It is coupled with the ordinary used contraction 
coefficient^^ by the relation! 
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where Mu- the minimum Mach number M of the propagation of the detonation wave 
(Mach nxmber M of Chj^jman-Jouguet) • 


The integration of the three-dimensional equations (l*2) can bo reduced to a 
successive integration of two-dimensional equations if we use the so-called 
method of ^termination of series” suggested by Van-Dike* This method was used 
in a number of works for the calculation of the flow of viscous [l] and ideal 
gases [2]* 

/ 

Let us represent the unknown functions in the neighbourhood of the critical 
line in the form of aeries with various powers of pin 9 • Taking account of 
the eveness of the functions we shall haver 


u(t^z^B) = ’U^cosQ *u^co^B sin^B * ••• 

= ... 


) 



The first ’♦temdnation" (the terra 'terroination'* is universally recognized) 
includes all the first terras of ejcpansion (l»5)i with the exception of the 
expansion for pressure in which the first two terras are taken; the second 
"termination"- the first two expansion terras of (l.5) plus three- in the 
expansion for pressure^ etc. 

The first "termination" is substituted in the system of equations (l«2), 
ai^d in each equstion the coefficients with powers of sin 8 are equated to zero* 
A system of equations concerning the variables t and Z is obtained for the 
coefficients of expansions (l*5)* system contains seven equations with six 

unknown functions* Follcndng work [l]f we shall omU the equation) obtained by 
equating the coefficients with sin'* 8 in the energy equation to zero* 

The final system of the equations of the first "termination" in the 
variables t and ^ 1st 

can be written in the following forrai 
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( 1 . 8 ) 


Here) ’Js= and — are the equations of the 

conto\xr of the shock wave and the body respectively. 
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Tho system of equations (l»6) has the following characteristics i 
ft) along the ^trajoctoryT 



(1.7) 
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dt ' dt 
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dt ’ ’ (l»9) 
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b) along tho characteristics C and C 
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c) along the characteristic j 
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The boundary conditiona on the contour of the body take the form* 


(1.15) 

For the calculation of the functions behind the shock waT6| it is necessary 
to substitute expansions (1*5) in the Hugoniot relational using the scheme of 
"termination” described above. Analogous expansion should be vrritten for the 
contour of the shock wave and the bodyt 

* ... . 

* ... , 

We shall be limited by the first two terms in expansions (l.l6) and shall 
assume that the coefficient with sin ^9 does not depend on time. This is 
equivalent to the assiunption that the shock wave (and body) in any point in the 
neighbourhood of the critical lino moves in the direction of the normal with the 
same velocity | and that it is near to a sphere. 

The relations coupling the parameters of the gas in frppt of and behind 
the shock wave (at ^ » l) are as follows! 
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For the solution of the system of equation (l.6) with the boundary 
conditions (l.l5) and (l.l7)i it is necessary to know all the functions at 
zero time. 
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The main calculations were carried out from the distributions obtained 
by solving the stationary equations (lo6)o The system of the stationary 
equations for the system (l*6) is solved with respect to the derivatives of 
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The solution of the boundary value problem for the system of the ordinary 
differential equations (1*18) is obtained by the iteration method '• The Cauchy 
problem is repeatedly solved from the initial data behind the shock wave before 
the bodyo The unknown parameter £g is determined by the Newton method from 
the condition of non~flowing (l»15)» The indeterminate form in the point ^ ■ 0* 
for some equations of (l*l8) is treated by extrapolation of the functions from 
a certain distance from the body. The integration of equations (iflS) is carried 
out according to the, clear Eulerian scheme. 

The system of the nonstationary equations (l»6) with the boundary 
conditions (l*15) and (l#l6) is solved by the finite-difference method within 
the interval O < ^ < 1 and with using the characteristic relations 
(l*7-l*14) for the detennination of the functions at the boundary pointso 
The space interval (0*l) is divided to (N + l) nodal points at equal interceptS| 
and the values of the functions for the following time-layer in the internal 
points are determined by the clear scaling circuit of Lax with a second order 
of accviracy (Fig. 2 a). At first, the values of the functions in the two 
half~integral points 5 and 6 are calcxxlated for the layer (j and then, 
by the obtained values of the functions in the half-integral points, all the 
functions in point 4 on the layer (j + l) are obtained. 

The calculation of in the half-integral and integral nodal points 
are performed by the solution of the ordinary differential equation (I.I 4 ) 
along each time layer, using the clear Eulerian scheme. 
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The withdrawal E of the shock wave from the body and its relative 
velocity £ on the following layer is found by the Newton iteration method* 

The trajectory of the shock wave within one time interval Is approximated by 
a parabolai and the value_of £ is the solution of the characteristic equation 
along the characteristic C drawn backwards from the boundary nodal point of the 
upper layer (Fig. 2 b). The scheme of the second order of accuracy is investigated 
and the value of the functions in the intersection point of the characteristic 
with the J- m layer is calculated by equations (l.ll~l.l2). The parameters in 
the point 9 are determined by linear or quadratic interpolation from the values 
of the functions in the nearest nodal points. 

On the surface of the body (Fig. 2 B) , the calculation is performed as 
follows. The pressure is originally calculated .on the layer (j + l) in the 
point l6 by integrating equation (1.14) at the end of the interval (O.l). 

The pressxire Pj and the velocity components are calcvLLated from equation (l.6), 
taking account of the bovindary condition (1.15) from the clear formulae j and the 
density is calcnilated by the iteration method from the equation that is a 
consequence of the relations along the trajectory and the C "*■ characteristic 
drawn backwards from point l6. The values of the functions in point 13 are 
determined by linear or quadratic interpolation. 
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After the calcxilation of the density^ the concentration is determined 
from relation (lo8) along the trajectory by the clear formula* 


^2- The algorithmi described in § 1, was programed on EVM-BESM-3M* 

A standard net with a number of nodal points N ■ 40 in the space variable E, 
was used in most of the calculations* Besides the standard net| a net with 
numbers of nodal points N «• 80 and 120 was used to estimate the accuracy of the 
individual variants* 

The main calculations were carried out at Moo ■ 5t T ** ^*4t <1 ** 0*3 and 
0*4* The activation energy E and the reaction rate L were simultaneously 
varied so that the relative width of the nonequilibrium zone of the stationary 
solution comprised about half of the width of the whole shock layer* The 
value of the withdrawal of the shock wave is close, to that observed in the 
experiments in vrork [7] and obtained in the calculations of stationary- flow 
around bodies [6], 

For the purpose of studying the influence of the hach number Moa and 
comparison with the known investigations of "one-dimensional flowsi the 
following values of the parameters were also considered! Moo ■ 7*6} 
q « Oo615 ■ 1*6), y » 1*2, B ■ O*5»l*25o Those conditions were partially 

considered on work [8], where the growth of the disturbances .with time in the 
one-dimensional flow in front of the piston was investigated* The values of 
£ 0*5 and 0*6 correspond to the steady and unsteady modes of flows according 

to the linear theory [9]* 

Two forms of disturbances of stationary flow were investigated* A strong 
disturbance was originally inserted in the Initial distribution* It led to the 
assignment of constant values of the functions of Pi I » P;^ I ->^1 and of a 

linear distribution for the velocity Ui with respect to ^ * The boundary 
conditions on the shock wave and on the body were satisfied, and the value of 
the withdrawal «* 0*2 is arbitrarily given* 

The calculation was carried on from these initial data to the setting of 
the stationary flow* The results of such calculations are represented in 
Fig* 3* Hero, for convenience, the dimensionless time f ■ t x (l*62)'^ was 
Introduced* This is related to the number of the time intei'vals, through which 
the results were delivered to the printing of the EVM* 

In the case of a perfect gas, and for small activation energies (5 0*5)i 

steadiness rapidly occurs. Slowly decaying oscillations with a nearly constant 
frequency are observed with the increase of the activation energy* Irrespective 
of the reaction rate, at E ■ 1*0, the generated oscillations due to the strong 
initial disturbance OTadually decay, after that the flow is described by the 
stationary solution U*l8)o 
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The change of the velocity of the incident flow by the lawi 


^ J" ^ . , i’Lon%t 

■ ■ ( 2 . 1 ) 

at constant density and preasure was with another form of the given 
disturbance* 

The growth of the diatrubance (l9) with time for some values of the 
activation energies and for different Mach numbers M is shown 'in Figs* 4 
and 5 » 

In the first casci at ■ 5 and for an activation energy E - 1*0| the 

oscillations always decay with time, so that this range of the values of E can be 
considered as belonging to the region of steadiness with respect to the finite 
disturbance* For an activation energy E « 1 , 5 , the oscillations intensify with 
time* 
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An analogous result was obtained at 7,6. When the activation energy E 
was 1 . 0 | the oscillations decayed, and at E ^ 1,25 they intensified with 
time* It should be noted that the palculations in the last case were carried out 
on a comparatively short time interval* This is connected with the fact that for 
great values of the activation energy, lying in the region of unsteadiness, the 
numerical errors increase owing to the increase of the gradients of the functions* 
For the continuation of calculations, a decrease of the interval of time and 
space from a certain moment was required# 

We investigated also the growth of very small disturbances connected with 
the increase of the numerical errors* For that, the calculation was performed 
from the stationary solution of equations (l*l8) without introducing any 
artificial disturbances* 

For values of the activation energy from the steady region with respect to 
the finite distui'bance, very weak oscillations are observed near the stationary 
solution* These oscillations practically do not increase by widely changing E 
from 0 to 1 * The deviation does not exceed ~ 0*3?5* For values of E from the 
unsteady region, a relatively rapid formation of considerable oscillations of the 
leading shock wave and the whole flow behind it is observed* The amplitude of the 
oscillations of the shock wave (Fig* 6) grows slowly with time, and the fi'equency 
remains almost constant* The qualitative character of .the cxurvo does not change 
by using a more fine net* The change of the profile of pressure from the shock 
wave to the body at "2 ■ 1— £ with time, shown in the same figure, illustrates 
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In order to exclude the Influence of the errors of the numerical method 
on qualitative conclusions about the growth or decay of the disturbances, a 
groat consideration in the process of performance of work was given to the control 
of the calculation accuracy. 

For the values of the parameters lying in the region of steadiness, besides 
the comparison of the results obtained with different space nets, the comparison 
with the stationary solution of equations (l*l8) served as a criterion for the 
accuracy of the numerical scheme. The small fluctuations of the parameters of 
the nonstationary but near— stationary solution insignificantly depend on the 
value of E ( 0.3yi) for sufficiently large instants. The comparison of the 

results obtained with the standard space net and with a not with half Interval 
gives approximately the same evaluation (see for example. Fig. 5). 

In some cases, the calculation with the standard net led to a physically 
imreal result. Thus, Fig. 3 represents the result of the calculation at E - 1.0 
and L - 0.14.10^ , obtained with N » 40 (solid line) and N ■ 80 (dotted line). 

In the first case, the shock wave raonotonically moves away from the body, in 
the second— the oscillations decay with time. 

In the unsteady region at E ■ 1*5 (M oo ■ 5)i the nonstationary solution, 
on the average, has fluctuations near the stationary. At the same time for the 
standard net, the deviation of the maxima of the fluctuations,' as a rule, is 
greater than the deviations of the minima* The relative error, as follows from 
the comparison of the results with N « 40, 00 and 120, is the greatest in th 
neighbourhood of the maxima. On decreasing the space internal, the curves are 
shifted to the side of smaller values of T .In this case one should evaluate 
the maximum error from the difference of the corresponding maxima. This 
evaluation gives a value of 5»7^ for the first three fluctuations. 

The numerical errors arise in the first place from the inaccuracy of the 
calculation of the kinetic equations (1.6) and (l.lO), where the product of 
the big (L) and the small (exponent) multipiliers is involved in the velocity 
of the process. In the case of the big gradients,^these errors rapidly increase, 
which loads to the physically unreal result. 

In order to exclude the influence of these errors, an additional series of 
calculations with a still more simple kinetic model was carried out. It is 
assumed that the mixture bums Instantaneously in the detonation front, and the 
heat release q of the reaction depends on the velocity of the flow in a system 
of coordinates connected with this front. 





- 150 - 


Such an enthalpy function changes the e:q>licit kinetic process in the gas 
(distributed heat supply, ignition delay, and others) to an implicit one and 
seirves as if it is an average characteristic of that process. The form of that 
function is determined by the parameters of the incident flow Vo« , Poo y 
and the dim^sions of the body. 

According to theory and experiment, the form of that function is 
qualitatively determined as represented in the left upper angle of Fig. 7* 

Here, a narrow range of the Mach number M is only shown, where the flow is 
reconstnicted from an adiabatic one without heat release to a detonation one 
with a complete combustion of the mixture in a relatively narrow zone. 





By carrying out concrete calculationa of the function cp (M), it is 
possible to approximate the analytic relation: 


The curves of the drop-out of the shock wave and the Mach number M of the 
incident flow in a coordinate system connected with the shock wave are plotted 
in Fig. 7 for the following values of the parameters: M ,>3 ■ 5.1, q 
- 1.4, K » 10, » 5 corresponds to q^a.^ " at cf - 1? 

The stationary solution with such parameters corresponds to a flow with a 
detonation wave, in which the mixture completely bums with the heat release! 
q ■ - i - q ^3^ <■ 0.3. Thus, such solution will be near to the limiting one 

at L for the case represented in Fig. 7f the difference here is in 

Mach number M 


OO 
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It is obvious that ovdng to tha numerical errors at "X, 0*759 the solution 
of the system of equations {l«6) becomes essentially nonstationary* It is 
interesting that in this case 9 self-oscillations with constant fi*equency and 
amplitude of the leading detonation wave take place* The mechanism of the 
maintainance of such oscillations is similar to that described in the experimental 
work [3]» 

Owing to some reasons 9 the stationary flow is led out from the equilibrium 
state at which| for example, the shock wave begin to move in the direction of the 
body* At that time, the thermodynamic parameters, pressure and temperature behind 
the shock wave, begin to fall in view of the decrease of Its intensity* Accord- 
ingly and because of the high activation energy of the fuel mixtures, the heat 
release sharply decreases and the shock wave acceleratingly moves to the body* 

Then its retardation takes place because of the compression of the gas in the shock 
layer "the reflection". As a result;, a peak of pressure and temperature 
is formedj the mixture again "ignites** and the shock wave maintained by ignition 
moves away from the body* Farther, because of the radial expansion of the flow, 
it begins to attenuate and it stops at a certain distance* After that, the process 
described above is renewed* 

Ail these moments of the process are represented in Figs* 7 and 8* The three- 
dimensional diagram of the pressure as a function of the two variables t and ^ 
is plotted in Fig* 8* 



Fig* 8 
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The described mechanism of the pulsation of the leading shock wave with the 
periodic generation of the waves of the burnt gas, which after that are carried 
away downstream, can serve in some cases as explanation of the presence of a 
strictly periodic wave structure of the combustion front* The latter in the 
supersonic part degenerates in surface, separating the burnt and unbumt parts 
of the gas. 

In conclusion, we mention that the values of the activation energy of a 
mixture, corresponding to the neutral oscillations in the investigated two- 
dimensional flow, are lying far in the region of unstability for one-dimensional 
flows in front of a piston. This follows from comparison with the results of the 
linear and nonlinear stability theory both for the strongly-compressed [8, 9] and 
weakly-compressed [4, 5] detonation waves. This effect is explained by the 
stabilizing influence of the spreading of gas along the surface. In this case, 
disturbances, do not succeed to develop arxl are carried away downstream. 
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USE OF THE BOUNDARY LAYER METHOD FOR SOLVING THE PROBLEMS 
OF THE MOTION OF GAS MIXTURES WITH EXOTHERMIC REACTIONS 


S.M* Glllnsklli and M«L« Khai^cln 


A great number of works [l^ 5-3] Is dedicated to the analTtlcal investigation 
of the motion of nonequilibrium reaction or relaxation raediao 

In the present work, a hypersonic nonequiiibrium flow of a mixture of gases 
near a wedge or a- cone, and also in front of a moving piston is investigated by 
the boundary layer method [2]. The obtained solutions can bo simply used for 
the evaluation of the influence of the oscillatory relaxation or the nonequiiibrium 
dissociation* 


§ lo Hypersonic nonequiiibrium flow of a fuel mixture of gases 
around a wedge or cone 



Let the heat supply to the gas takes place in a combustion wave of a finite 
width and lot this process be arbitrarily described by one irreversible reaction 
[3]t 


, m m-1 



( 1 . 1 ) 


w^ere /3 - the relative concentration of the original reacting components of 
the mixture, n>- the order of reaction, p - pressure, T — temperature, 

E - activation energy, R - the imiversil gas constant* 


We shall be limited to the investigation of plane and axisymmetric flows* 
Following work [2] for the variables of the boundary layer, where ^ •— the 
coordinate along the contour and •- the flow function, wo shall write the 
gas-dynamical equation system ' 


du 


dv- dp 




9x 


dx' 
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/ 


U i-t Sp 

= -Z. 
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■■ 


U ip 


m m-f , Bp \ £ 
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( 1 . 2 ) 



155 


Here, ^ — coordinat.e, orthogonal x; u,\r«- the components of the velocity vector 
in the' directions xS !J‘ 1 - radius of the curvature of the flow lino, 

V ■ 1,2, for plane and axisymmetric flows, M— molecular weight of the mixture, 

Q — quantity of heat supplied to the unit mass of gas at complete combustion* 

For the axisymmetric flows, — the distance from the axis of symmetry to 
the considered point of the shock layer, ? — • the distance from the axis of 
symmetry to the point of the contoxir (^‘ig* l)* The flow function is interpolated 
by the relations: 


ay' _ 1 

3^ ' I 



(1.3) 

(1.4) 


The solution of (1.3), (l.A) ie written in the form of series of powers of 
the small parameter £ s 


y'cfy;*,.. 




O t 






(1.5) 



The variables of the boundary layer j diagramatic pattern of 
flow near the boc^- 
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After the performance of the standard operations for the determination 
of the first ter my of the series j ve get a system of equations i: 


_ 1 ’ 



fo . j.'i-' 3Po 
R 5 * * 


( 1 . 6 ) 


and for the determination of the functions ^ y ’Pi haver 



(l<»7) 


The representation of the concentration^ in the form of (l»5) ia a 
consequence of the assumption that the limiting flow at £ 0 is a non—' 

equilibrium one, and this enables to take into account, as a first approximation, 
the effect of the nonequilibriura state. 

The arbitrary functions shoxild be determined from the boundary conditions 
on the given contour 
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( 1 . 8 ) 


and on the shock wave, where all the functions are expanded in power series 
of the form of (l»5)* 


The relations for the shock wave in the considered case coincide with 
those given in work [2] for the ordinary adiabatic flow, and it is necessary 
only to add the conditions of concentration (the index S will bo attached to 
the parameters directly behind the shock wave)* 

- 0 ~Zt 

* 0(6^) . (1,9) 

The systems of equations (l*6) and (l*7) are integrated in quadratures, 
but unlike the case of the adi^atic flow this procedure is more complicated, 
and the solution is more cumbersome* 


We introduce the dimensionless variables x, ^ and as follows t 


, v'--L-\,pry , 
(l u^pJ "'iLnU) , 


(i.io) 


In this form, the function x* (sp) on the shock wave becomes ^ 

For a wedge and cone, the integration of the last two equations of (l*6) 
along the line « const gives: 


where 





( 1 . 11 ) 


B- 




— , c 1 + 


Po, 




( 1 . 12 ) 


In the most simple case , at m ■ 1, for which all the results will be 
reduced for shortness, the integration of (l*ll) gives (see [4]). 
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= e°Ei(i; 'a}-Ei(^ 

CO i;" 

tni; *1. ~I , , t;~a>0 . 


( 1 . 13 ) 


Formula (l*ll) gives an unclear dependence of the concentration 
on the coordinates of the point (x» ) in the sKock layer to zeroj and only 

in the case when the activation energy E equals aero ($ ■ O), it is possible 
to solve this relation with respect to ^o' 

When the activation energy differs from zerof it is necessary to select, 
as Independent variables, the value of the concentration /Pg inside the shock 
layer and the value of the concentration >3® on the contour in such a way 
that the coupling of the old and new variables will be given by the equations i: 


(l.U) 


Further, it is necessary to express all the unknown functions in the 
variables p© j i^i an obvious forra, end then, for example, plot graphically 
the dependence on x. and y_j , It appears that all the unknown functions of the 
first and second approximations are expressed in a combination of exponential 
integrals, transcendental and rational algebraic functions of and 


The calculation of those functions reduces to a partial differentiation, 
or. integration along the coordinate lines ^ ■ const, y-* ■ const. 


FinalIyc_,we_obt ain for a wedge in zerO"a pproximatiQn> 
= y coU , P^--P^V , 


sin 2c<, ‘ 





ff ■ 
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C ^ 


K 




if 




a C 
c- 1 


where y (x) determines the form of the shock wave. 


( 1.153 
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In the first approximation, we shall haver 

casU 


/'f 9f2c-f-A^>jf/>/’-G/r/-A"Aj7 ' 
*/’ />oj * 0^]} -p^ » 

;S^ -- exp[a/ff-^Jc)][‘Pfi;)-^(' ^i)J , 

» • 


(l*l6) 


(1.17) 


■<*> "^Ol 


It is obvious from the formula that the concentration along the flow line 
monotonically decrea3es_from one to zero ( Fig, 2). 



0 to 10 iO 1,0 sc to X 

Distribution of the concentration y3 (x) along the generating lino 
of the body 

first approximation 

second approximation taken into account 
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In the first approximation, we shall haver 


A -A^^p[^/0-A/c)][‘P(-^)-^(-^J] , 




wHa-t*e nr -,r 

•*•(■ ‘^)'-[i/i"J,JfJ]["^if°-‘:h f ]■ 

* 

^ • [C •K*{^r)]b'°Bi(cx'K)-Eih)] 

i; >0 


cos 


td.i 


'%‘-^os Poj J 


,!■-< 


(1.17) 


It is obvious from the formula that the concentration J3 along the flow line 
monotonic ally decreases-l^o P one to zero (Fig» 2) » To 



Distribution of the concentration p (x) along .the generating line 
of the body 

— first approximation 


second approximation taken into account 
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Formulae (l»15) and (l*l6) are essentially simplified and become evident 
with respect to ‘X , when the activation energy is equal to zero. For 
shortness, we shall not give these formulae here, since analogous formulae will 
be written out below. 


It should be noticed that the pressure and velocity in the zero approximation 
are calculated by the Newton formula* Talcing account of the first approximation, 
the velocity u is constant along the flow lino, but it changes across the shock 
layer. The pressure changes nonmonotonic ally and it has a maximum. On increasing 
the activation energy this maximum increases and shifts to the equilibrium region 
of flow. The curves ■ f (x) along the wall of the wedge are plotted in 

Fig. 3 for the values of B ■ 0 and 0*7 1 C » 1.3 j; 





( 1 . 18 ) 


Pressure distribution along the side of the wedge for two values 
of the dimensionless parameter 6 

Taking account of the first approximation, the density and concentration 
monotonic ally decrease along the line ■ const* For the determination of the 

function 'io CPo j PoO havei 





A' 
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1 -^ A ^ M -*<^‘)1 ' 

wV|erti X(^)^°^c-I^A[a’e°£i(i;-a)-2ELfi:}*' ^ ]* -^[^ Ei(i!-a)‘Ei(i^)] * 

\ ' • ' ■ 


K-[x~ e"" £i(t;^-a)*£i{K^)] , ^ */ — ~r 
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We present the function § in an integrand of the imetationary integral "i 
in the form of a power series in the neighbourhood of the point a* Performing a 
term-by-term integration of this series we shall hayej 


^ ^-z I M ~r 1 , 


where ‘J is calculated by the recurrent relational 


(1.21) 


' K-t If’f f 

^ f /*/ 2f^*o) tt ^ 


The component of the velocity is calc ulated as follows i 


J esc. . _ />/^ - 


( 1 . 22 ) 


( 1 . 23 ) 


The position and inclination of the shock way® to the x direction are 
determined as followsi: - 
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y. 




X tin 2d 


yl ' Si'' ^ - — 

® OO 9 V 

•To j ^ tt 


( 1 . 2 /,) 

(1.25) 


In the zero approximation, the difference of the aolution for a wedge 
and cone appears only in the form of the shock wave and the value of the 
velocity and the velocity u.e }; the pressure, density and concentration 
in these solutions coincide o 

It follows from system (l.l6) that in the f irst a pproximation for a 
reaction of thel first order, the parameters Pt,Pi< are 

determined for a cone from the same differential equations recorded along 
the flow lino for a wedge, with an accuracy up to the coefficients which depend 
on and are determined in the shock wave. Therefore, for those 

functions we shall not give here the tedious formulae in the case of 6^0* 

Wo shall write out the more simple formulae for the case of B - 0, 
when all the functions are evidently expressed in x ^ , 


1 1 


zero approximation 
u^-V^cotd , ^ 




u. =-• 


\/« 


2coioi 


first approximation 
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Jk iinoi , 
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'*^1 , ’/t 


sin^al f 1 


The bohavicmr of the gaa-dynamical parameters in the shock layer near 
a cone I taking account of both approximations! is qualitatively similar to the 
behaviour of the corresponding parameters near a wedge* 

As an illustration! a comparison of the position of the shock waves for a 
wedge ani cone with the same expansion angle 2*»< ■ 60® at /^oe>- ! q ■ 0*5 

is drawn in Fig, 4* The shock waves for frozen and equilibrium flows are shown 
by the dotted and dash-dotted lines! respectively* 

The pressiire distribution along the wall of the wedge and the generatrix 
of the cone for these conditions is represented in Fig, 5, It is obvious that 
at the same activation energy the pressure on the wedge tends more rapidly to 
its equilibrium value than on the cone. 



Form of the shock wave for a wedge ■ l) and a cone ( v ■ 2) with 
expansion angle m 60®! Mpo=«> i <1 ■ 0.5, 



1 ' . ■ 

1 Pressure coefficients for a wedge and a cone at 2 - 60«| 

, i ' ,q»p.5o 

In the else of S - 0, the pressure on the wedge and cone is caculatod 
by the form ula e i 


r 1 P'P- 
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* (1.29) 

«hloh =c„.=rt to tho l«o» ror^ao for th. adiabatic flow around a wod*. 
OT cono, ir wa put in then q - 0 (aea [2]). 
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§ 2, Nonegulllbrium flow of a fual mixture of gases, cauaed 
by a piston moving with a. high supersonic velocity 


to § 1, it is possible to investigate the problem of the 
motion of a mixture of gases, generated when a piston moves with a high 
supersonic/ velocity. The original system of equations in ILagrangian 
variables has the formj 


iL - ^ 

d m ' 


dP u n e / £ ) 

‘‘Pi'-pTel . 

where — the fjiitial distance of a particle of the gas, by which the 
shock wave did not pass yet, to the plane of the axis of symmetry} 

Q - the quantity of heat released on complete combustion of a unit 
mass of the gas, 

0 

-P- density, F* - pressure of the unburned gas. 


(2.1J 



We shall introduce the parameter £ - , following £2j v/e shall 

seek for the solution of (2.1) by the method of the bound^y layer in the form 
of the following series J 


(2-23 
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The boundary values for the zero terms are obtaljied by the analoguous 
expansion of the boundary values of the unknown functions. For the zero and 
first terms of series (2.2), we obtain certain systems of equations. 

For the zero terms:: 



C2.33 


where K - K.£ . Frbm (2 .3)1 it follows that « Ro(fc)* As R. (t), we 

shall take the law of motion of the shock wave. Then the equality R^ ■ 0 
is fulfilled on the shock wave. Lot us turn from the Lagrangian variable m 
to the Lagrangian variable r using the relation. 


m = 




From this relation, it is obvious that la the time of passage of the shock 
v:ave by the particle with the coordinate m. The system of equations for the 
first terms has the following form: 


aK, 1 

dm ’ 


at^ 
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Lot us write the system for the zero and firstterma in a diraonaionlesa 
formi: 


4 










where T— a certain characteristic time* After simple convorsionsi the system 
for the "zero” terms will have the following fomn 








If- e,p{- g j 




( 2 . 6 ) 


Here R -a — the law of the movement of the piston* The third and forth equations 
of (2*5) are an integral-differential system of equations with respect to the 
\inknown functions fio and R^, o For solving this system, the approximate method 
of solution was used, which is algorithmically realizod_on the computer as follower 
the iteration method was used, in which the value of R was taken as the zero 
approximation for So} was_ determined from the third equation, and then the 

following approxin^tion for Ro was determined from the forth equation# The K 
approximation for Re, is expressed by the formula 






Here R ti the law of the movement of the piston. The third and forth equations 
of (2.5) are an inte^al-differential system of equations with respect to the 
unknown functions fto and o For solving this system, the approximate method 
of solution was used, which is algorithraica-ll^ realized_on the computer as follows*: 
the iteration method was used, in which the value of R ,v was taken as the zero 
approximation for was_determined from the third equation, and then the 

following approximation for Ro was determined from the forth equation. The K 
approximation for Re, is expressed by the formula 



C 
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The iteration process was cut off on achievement of the condition: 




where S — a certain constant* 


After that I '^the third equation is solved with the found i which is 

taken as the law of motion of the shock wave, i,e« is found, and then 

the remaining terms of the zero approximation are found* 


The system of equations of the first approximation has the following formal; 
1 ‘ f. 

p, JK ■ KafCm • pr. / . 


o X 




i' ifl ) i , />_ £<:' £1. y,f£, 2i I 
H\ 3t‘ A ' p./fjp^ p./fJrfJh/ pfii 


=r 




For simplicity, calculations were carried out for the case or piston movement 
with a constant velocity* 



Hg* 6 

D’ependence of the pressure on the piston upon the dimensionl 
time t. V ■=. 1, 2, 3— 'Pln^e, cylindrical and spherical cases, 
respectively. 
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Law of motion of a plane piston. 

/' 

I 





D.p«.denoe of th. on the pieton upon 
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The graph of the distribution of the pressure on the piston at different 
values is shown in Fig» 6. Tne graph of the law of motion of the shock wave at 
different V values — in Fig, 7, The distribution of the pressure at a fixed 
instant on the"^ coordinate — in Fig, 8, The heat distribution on the piston 
at different values— in Fig, 9 , 


The problem was solved for the following numerical v^es of the parameters. 

I C-ri:Q ; ^^=6 ; Q * 30 ; B = T = fO ® . 

For the| investigated case (Ro w Dt)| from (2.3) we have»: 


+ C/rj , whet''; ) 

• ?o-PJ(c(x)-2q^,) , 


Substituting in the forth equation of (2,3) t we iobtain.ah equation for 
the determination of 


dt 




The boundary condition for the shock wave is as follows i; 
at t-.t 


Since Rq ■ Dt, then C (t:) is a constant quantity. It is obvious from the 
equation that ^0=^/3^, therefore all the values of the zero approximation 

are functions of t - t; ♦ In the zero form is expressed by the following 
relation: 


f ? 

expi /• 

/ f 


% 


(2.6) 
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Jtb a fixed X and t — » co the integral on the right side is dropped, 
therefore the integral on the left is also dropped© At m ^ 1 (this case 
will be investigated) this is possible when at a fixed 't and t — >•<» 
we have , i.e. any particle bums for an infinite time. From the 

fact that po m f (t -Z), ±t follows that for any particle the process of 
combustion will equally proceed if only we take t ■ "C as the start of the 
accouiTi; of time for the particle. As the integral on the right side of (2.6) 
is not taken in an evident fonii, then for finding we take advantage of 
the approximation method, using the following consideration. As )-> 0 

at a fixed value of TT and t — » oa , then it follows that for any value of 
and y there is such a value of t that at t > T + t» it follows that 
/3o <■ c-p , t^ does not depend on TT and ('f ^ • In particiilar, 

wo take r . Therefore, at t > t’ it is possible to neglect 

in comparison with 1. Then for such values of t wo obtain the following 
equation for tne determination of /?£, 



- m n f f Ec' ) 


(2.7) 


here ^ ■ const. Integrating equation (2.7), wo have; 




r 

•fn 


hare *2-0 known function of t*. At the interval 4 *7^ t't 

Wo approximate ^ to a paroblat 

a, b, Z, t’ are unknowns. We designate 

“■ S.?\ 




Then the four conditions for the determination of the four unknown constants 
take the formi 


2jt + , 


3j 


li-.vA' 
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moreover, V « T , where Zp — a known constant* 

From equ3j,tion (2*4) i it follows that at J-’ ■ 1 the following relation 
is fxafilled: 

The third equation of (2*4) is integrated at once if the following 
equality is taken into account i 


so that 


where 


> f! H 3t 


1 


Let us express from the given relation in /S^ and substitute the known 

”zoro” and "first” terras in the forth equation of (2*4) • We then obtain an 
equation for the determination of 


» 0 




The solution of this equation at 




c.ft, tj = - *p'tf2>‘)(t-zhe(i.zh 4= f-. zq ■ 


'f. 


y>/i-££-y' f dt^Cfc-zqJ/^Jt -(c-!Q]Er’f(ff Jt . 


•'I p 
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where 
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Here the formula for is not given because of its tediousness # With 

respect to we only mention that as at 

then for any fixed °T? and t > oo • The law of motion of the 

piston R^(t) is determined as fol low s: 

i ® 

1 /fim . 

£-»oo " 

Let us consider now some qualitative peculleu'ities or tne found solution. 
We shall show that after the lapse of time Zc from the beginning of 
combustion, the characteristics of flow with an accuracy up to become 

constants equal 'for all the particles:: 



2)^ 



/,<£ ' cit 


t > Z * 

O 




and *0(i^) ■ 


?^^cond * O(i^) ; because 
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^ 2)4*^' 


ufT,t)z— = 7 ) * 

at 2) 


rry^/^’/o,Po.p,) ; 

at 


uCr,t)zT).t£ ^ofc^) 


It should be noticed that since all the functions depend only on the 
argument t then after the combustion of the particles the flow (with 
an accuracy up to ) will be isentropic. The equation of the line in the 
plane R, t, by the attainment of which all the characteristics of flow 
become constants with an accuracy up to i is written in_ the form:: 


lR--R(i,z) 

1 


'RdR(rr-rj 
2Q 
2)' 




3 i ! 

^ Ji. 


Erow which 

1 c 


The given line (we shall call it AB) is parallel, in the plane R, t, to the line 
representing the shock wave (Fig. lO) . In the region lyJjig on the right side 
of Al3 t it, is possible to establish an isentropic flow, assuming that the gas 
is definitely burnt. In order that the solution can be plotted in the region ABC, 
realization of one of the following conditions is necessaryt 
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l) The line AB must be s characteristics as the characteristics in the 
region ABC must not intersectt the line AB, because all the functions are- constants 
in-'ito^ The relationj 


R = ii + a . 

\ 

is realized for the characteristics in the region ABC* Hence, it follows that 

on the line R « D; 

\ B “ u + a auTT ' 

\ ■ ■ 




where a— the velocity of sound. Thus, the velocity of the shock wave, at 
which the line AB is a characteristic, is uniquely determined!- 


■ ( 2 , 8 ) 

In such a case, instead of the interval AC, it is possible to take 
another law for the motion of the piston, and the solution in the region ABC 
will be a simple rarefaction wave} 
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2) The line AB may not be a characteristic. It is possible to plot the 
flow in the region ABC by the use of the characteristics! and as the piston in 
the case of t y Zc moves with a constant velocity with an accuracy up to £ 
then this flow will be a progressive one. In this case the characteristics 
will intersect the line AB. In order that the inclination of these characteristics 
was larger than the inclination of AB, it is necessary that the following 
inequality is 'realized; 


-7,^ c//?l 


or 2 




2 

+ * 


• • e. 


In the present case there is an analogy with the detonations, when the 
mixture instantly bums in the shock wave. From the great number of possibilities 
of the propagation velocities of the detonation wave, given by detonation 
adiabat, the velocity of the propagation of the wave of Ch.-J. is minimum* 
Therefore, it is possible to make an analogy between the velocity of Ch.-J, 
and the velocity 'of the shock wave, which is given by (2,8), Since in the 
case of K —^00 the time of the combustion of the mixture approaches zero, 
then the obtained solution must change to the solution of the problem of the 
motion of a plane' detonation wave, and the velocity determined by (2.8) must 
change to the velocity of the wave of Ch,-J. In the case when the counteiv 
pressure is not taken into acco\mt we haver. 




In our cascj the velocity of the wave is calculated as follows; 


where 



This discrepancy shows that the method used for the solution of the problem is, 
generally speaking, unfit when the velocities of motion of the shock waves are 
near to the velocity of the wave of Ch.-J . For the highly supercompressed 
waves, the heat release enurts a lesser influence on the character of flow* 
Therefore, the solution obtained by the method of the boundary layer [2^ for 
such velocities of propagation of the shock wave will give smaller difference 
from the exact solution with a detonation. The velocity obtained after the 
combustion of a substance can be used for the calculation of the limits of 
applicability of the boimdary layer method in the case of the solution of a 
problem with chemical reactions* Thus , it is sufficient to compare this 
velocity with the velocity formed behind the ordinary detonation wave. The 
velocity behind the detonation wave in the exact solution is; 
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u = 






( 2 . 9 ) 


The value of Uy corresponds to the wave of Ch.-J. 


„ § 0 -^' 

y ' 

In the case under consideration the velocity of the gas after its combustion, 
with an accuracy up to , is determined by the formular 


It is easy to see that in this case whenj 

1 


pii j 


<< 1 


( 2 . 10 ) 


( 2 . 11 ) 


formulae (2,9) and (2.10) coincide with an accuracy up to i.e. the 
velocities behind the detonation wave and behind the region in which the 
mixture bums asymptotically approach each other at -f — ^ Since V 

in problems of the considered typo, then £ , Thus, inequality (2,ll) can be 

considered as the criterion of the applicability of the boundary layer, method 
[2] in solving problems of the considered class. 


In conclusion, it is wortlunoting that the law of propagation of the shock 
wave in the case of spherical and cylindrical symmetry, corresponding to the 
steady mode (R n " const), has the forraj 



Thus, the deviation from the mode of Ch,-J. increases with the increase of 
the value of , 
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LDJEARIZED SUPERSONIC NONEQUILIBRIUM FLOW OF A FUEL 
MIXTURE OF GASES HEAR A WEDGE 
( 

-By 

S.M, Gilinskii 


Lot us consider a supersonic flow of a fuel mixture around a wedge. 

We shall assume that the maximum quantity of heat| which may be released as 
a result of a chemical reaction, is much smaller than the total enthalpy of 
the mixture in the incident flow. In this case, the nonequilibrium flow will 
slightly differ from the adiabatic flow near the wedge, and the problem of 
the flow around can bo solved analytically in a linear approximation. 

If wo consider the profile near to the wedge, the disturbances caused by 
the bend of the profile will interact with the disturbances related to the 
nonequilibrium release of heat in the flow. In the linear approximation, the 
total disturbance is determined by a simple superposition of those disturbances. 
Therefore it is sufficient to study independently the two problems of the 
nonequilibriiun 'flow around the wedge and the adiabatic flow around the profile 
near to the wedge. 

The first problem was investigated in works Cl-3» 6, 7]i the second- in 
works [4, 5i 8| 9]. 

The solution of the problem of a flow of a fuel mixtvire around a; wedge is 
given below, taking into account the proceeding of one irreversible exothermic 
reaction. In the general case the solution is presented in the form of a 
series. For certain particular values of the oinier of the reaction jtiiis aeries 
is summed up, and the solution is expressed by finite algebraic functions. 

The latter enabled to analyze the influence of the different laws ofsupply 
of heat to the gas on the character of the flow behind the shock wave. In 
particular, it is discovered that the solution may have an oscillatory character. 
The amplitude and frequency of these oscillations decrease downstream. 

In the extreme case of hypersonic flow around a thin wedge,, the qualitative 
character of the solution essentially depends on the order m of the reaction. 
Thus at m - 0 and <» , T-* Q — C, the amplitude and 

frequency of these oscillations increase. The amplitude of the oscillations 
approaches a finite limit, and the frequency approaches infinity. In the 
case of m - 2/3, the amplitude of the oscillations approaches zero. 

It is convenient to write the system of the gas-dynamical equations in 

the Cartesian coordinates x i q (Fig* l) ♦ 

1 t 
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Coordlnato system and designations 


Let the parameters with the index "0” correspond to the adiabatic (frozen) 
flow near the wedge# We-r.epr.esentjthe unknown functions in the form: 

y,}] , 

« 

p/*( . ' P- p ^ 


,‘/f) =fc[^ ’ 

where the index ”1” is written for the parameters of the unknown disturbed flow, 
and the functions of the disturbances are designated, without index. 

Vie shall model the kinetices of the chemical reactions by one reactions 




where y3 — the relative concentration of the original reacting components of the 
mixture, m*- the oriier of the reaction, P— pressure, T— temperature, E— acti- 
vation energy, R — the universal gas constant, L-- reaction rate constant. 


- 180 - 




Substituting (l) in the system of the gas-dynamical equations and 
neglecting the terms of second order of smallness with respect to the 
disturbances ^ we obtain the following system of linear equations r 


au aj. o 

P + U + P — ® ^ 




I 


Where 


£ii + if. = 0 

dx 


dv 

dx. 


if 


-- 0 


Z dp dif dp 

M — - f +5 

ax, ay, 


ax. 





% 


(fi]Q 


( 2 ) 


Mo— the ^ Mach number | a— the velocity of the sound of the. undisturbed flow 
behind the shock wave. 


Owing to the assumption of the smallness of the dimensionless parameter q^| 
the system of equations (2) contains a derivative of concentration only 
•j3o which is determined independently from the parameters of the 

disturbance, and it is a function of the undistrubed flowo 


3/ 

^"ax. 


0 ^ 


I m m - 1 




(3) 


The boimdary conditions in the linear theory are sot at the boundary of the 
undisturbed flow, i >e# 


(4) 

at 9^ ^ ^0 

u--Kajf'K) . 


C 
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.Here,. _! y,i}a the equ^ion of the disturbed shock wavei 

I , K p , Kj, —known functions of the Mach number Moo , wedge angle 9 
and the adiabatic index ) which are determined from the linearized system 
of the relations in the shock wave. 

I 

For example r 

We introduce the dimensionless independent variables x, y as follows J 

’■ Lpr' ’ Lpr^ ■ ( 6 ) 


Integrating equation (3) taking account of the boundary condition (5)| we get: 

: ^ 


I ‘ (7) 

The order m of the reaction determines the different laws of supply of heat to 
the gas. In the case of m 1, the width of the zone of the reaction extends 
to infinity and, according to (7)i the profile of the concentration /Pp is a 
continuous function. At m < 1 the reaction ends at a finite distance. The 
profile of the concentration as a function of the value of m may have a 
rupture in the end of the zone of the heat release. This rupture is either a* 
discontinuity or the curvature, or of the leading derivatives (Fig. 2). 



Distribution of the concentration pt> along the line of flowj 
m— oixler of the reaction . 
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For convenience, we determine the concentration by the function 
as follows j 








( 8 ) 


Excluding the density j° and velocity u from the system of equations 
(2), we obtain nonhomogeneous wave equations for the determination of the 
components p and \/' 




' =-? mB, 

ax^ 


Zm- 1 


(9) 


Solving equations (9) and taking account of (2), the functions will be:; 


v = 


P = 


% 






( 10 ) 


For the determination of the unknown functions F and G, we use the boundary 
conditions (4) and (5)» From condition (5) for the shock wave it is possible 
to vrrite, in particular, the following equations: 

p = ^ p = Su 

Kn Z co% ( 9^ * B) ii-n B ^ 


A - 


_ "p _ 


cosfd * 29^) + cos6/tA^ 8) 




2cosf9.*9) sin 9 ^ 


Kj, *29^j * iir.Q/t^}^sln^[6^*B) 


(n) 
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Conditions (4) and (ll) give the following system of functional 
equations _ 

, fix] -G[x)^0 , 

^ (X2) 


Excluding function G from this systemi we obtain one functional equation for 
the determination of F: 



where 






'Os' ^ 

I A 


9 ■ f 


Cu Aji](ct^^9^~ jx^) 7 


f 


1 




) - A^t- 
1 + 


ct^Q, *}i- 


sin (d, • B„) 
sen ^ + 9,] 


oC = c>" :en ///v^^ 


(13) 


The functional equation of the form C 13 ) . was studied in a most detailed 
fashion in work [4]« The parameter 'K is called the coefficient of reflection 
of the low disturbances from the shock wave, and it is equcCL to the ratio of 
the amplitude of the disturbance of the pressure reflected from the shock wave 
(along the characteristic x + /t- y •« const) to the amplitude of the incident 
disturbance of the pressure (along the characteristic x-^ xj b const)* The 
absolute value of the reflection coefficient A does not exceed the unity* 

Its value strongly depends on the adiabatic index and at ■=. oa ^ i A_^i. 
According to (l3)» the coefficient K is positive and does not exceed unity* 


By solving (13) we shall have the following series ([4]) 


Then from (l2) we determine the function G 


(14) 


r/-ci- y 


(15) 



Substituting the functions F and G in the formulae (9) and (lO) 
we shall have* 


oo 

*Act^ B^S^[x-y.^y 

; (16) 


p = 6^ - xj- :\u-td^ Oj - 


■jz 


3 B X 

C 0 I 




t AcIqB^Z HB„[K'(x^0)]0jl<^fx-0)]J) 

“hi'- 

A^ l(ctfs^ . (17) 

It should be noted that formulae (l6, 17) are equivalent to the formulae 
obtained in [ 6 ^ 7 3 the particular case where m ■ 1* 

From equation (2), it follows that 

" ^18) 


where S (y) — a function proportional to the vortex. The function S (y) is 
deterndn^ from the boundary condition on the line ctg — >. and has the 
following formt; 

+ Z A '[&Jq icUjd^ ^ u) * (tj(ci^ - ii))] j , 

(19) 

One may notice that on the svirface of the wedge the function of vorticity 
becomes zero. Across the shock layer this function rapidly increases, so that 
in the neighbourhood of the surface of the wedge there is a thin vortical 
relaxation layer. The existence of such a vortical layer was previously ; 
pointed out iri_ [l, 2, 6, Vl^ /The form of the shock wave is determincsa, 
if the boundary condition (5) is used, for example 

P^Kpf'(^) 

setting the following equation in the argument of the functions appearing in 
the formula for p (l7)f 
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The pressure distribution on the wall of the wedge is of particular interest 
as its knowledge enables to calculate the aerodynamic force of the characteristics* 
Below we shall be limited only to the analysis of this function. The qualitative 
character of the behaviour of the remaining gas “dynamical functions on the 
wall of the wedge will be as follows j 


where 




f h ^ K - 0 • 


BJx] .ig fJ^)l 



0 


at 

at 


( 20 ) 


xfe 




f 

1 - m 
I 

l~ m 


The equilibrium value to which the function p {'x- jO) tends ait x ■ — > ^ 
is equal to: ■ . 

j act^Byi 

~ P° * •Po ct^^d g- ' (2l) 

Formula (20) is simplified in the case of a flow with a very high velocity 
around a thin wedge, i.e., when the following relations are realizodf 


5»/ n«i 


( 22 ) 


Using condition (22), it is possible to determine the limiting 
values of the characteristic parameters of flow 

77[r^~T-!^ 77i^ 

"» 2y . ’ ' Jt2S ’ IfS ’ 


„.8.=s ,.a.-- 2 s = 


and formula (20) is reformed as follows: 


( 23 ) 


- ■■ -• ■ y fy 1 


^ i-.i 


„ °° 
Z r~> 


Whore 


( 24 ) 
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At the edge of the wedge at vc s 0 ^ the floy is approximately 
frozen and X (O) =0 m ^ X (x) — > 1 . 


The calculation of the function (x) can be approximately perfomed| 
being limited to a finite number of terms in the series. The particular casesj 
where the sum of the series is calculated in a finite form are of specific 
interest. Obviously, this can be done if the order ra of the reaction equalst 


if. n . ■ 

In this case, the width of the zone of combustion is finite and the 
distribution of the concentration in this zone is in the form of a polynoraia.!' 
of the power 5. The function % (x) in the points x= 0, 1, 2,...) 

as a function of S has a rupture, either a discontinuity of the cervature or of 
the leading derivatives. 


Performing expansion in the seriesi 




m-iW-2) I* 
Ti 7 ^ 



and performing the summation at the same powers of x , we obtain the following 
formula for the interval of change fi x 4 t '- 


I 


X: 1 I / • 


II 


I 


or taking (23) into account, we get: 




' 2 ~r (usF) * B£^(U9S^) 


Zx 3 (£-0 


( 25 ) 


.] : 


( 26 ) 


If the order of the reaction m ■ 0, then 'V-(x) is a linear function 
increasing from zero to the maximum value at the end of the zone of combustion 




/ + 2S* 




Z(J- 5^} 


] t 2S^ 


( 27 ) 
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In the case of m ■ the curve 7c (x) is a quadratic parabola with a 
maximum within the interval_0_.c % which_exceeds_the unity 




^ U5S‘ 


] 


V = i f 1 I — 1 

‘-^ci 3L2 

The abscissa of ithe raaxirauijj increases raonotonically with the decrease of the 
parameter S| arid the value of the maximum changes nonmonotonically with 
decreasing S. ' 



In the case of ra ■ 2/3, the ciurve *%(x) is a cubic parabolai the 
maximum is attained also within the interval 0 <1 x < ^ 


I 


'2x 




;*5S^ 27 1 


— " 1 


(29) 


1 

It should be mentioned that in the limiting case S « 0| the inclination of the 
quadratic parabola at the end of the interval is equal to unity, and that of 
the cubic one is equal to zero. 


On further increase of the order of the reaction, the maximum will be 
shifted towards equilibrium and its value will decrease. At the same time, 
with increasing the power ^ of the polynomial , the number of these maxima 
and minima increases. The solution in the zone of combustion has an oscillatory 
. character with rapidly decreasing amplitude and frequency of the oscillations. 
Thus, at m - 1, i.e. in the case of exponential distribution of the concentration, 
the first maximum is attained at xxs 3.2 and it equals 1 . 062 } the second 
minimum at x 7:1 12 and equals 0.998,.... 

The existence of weak oscillations for this case was discovered in ^3"}. 


In the case of the determination of .the function % (x) in the interval 
between the discontinuities ^ ^ "Vk necessary to 

equate the first terms n of the series (24) to zoror 


m = 0 


1 (0) 
[ 

«o 

W-/-fZ 7 

i s/i» f 

Y/'K‘xj= f- 

2l"^y ; 

. “""x/ 

s Lux 

UKX J 

fm * V 2 

1 t 


^n*t 

r r K 

^Zn*2 


L 

/-A * I'Kl 

/-x'A 

7 -I 


(30) 


(31) 


.V 2 




,3n*3 


3 ?7 


]. 


(32) 
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The values of the functions “X (x) in the discontinuity points are determined 
as follows: 


I f]= / ^ 4 ~ " 

^ -I (-5 / j lf2S-i UZh J ’ 


^ n/-s/i''7,,2s2j/-,,5g2jr,.2s/ ’ 


(33) 

(34) 
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It is obvious that the values of the functions X (?■} in the discontinuity 
points fluctuate about mity and the decay coefficient of these oscillations 
is equal to the reflection coefficient of the distvirbances from the shock wavef 
the frequency of the oscillations is determined by the parameter K =i 4^ ' 

L ^ S 


The abscis sae of the p pjjats of the extrema- function are equal to: 

-ro) ^ 

^(tiJ . ’ (38) 

f'U) 


U5S- 


and the values of the extrema equal 


f] 


(37) 

(38) 


■ ^roj {-hs^ / /-25)^" ~ 

’ r + 2S^ (” J+ 2S / 

Sf/fSy/-2S) S ?-2S j'' 


(30) 

(40) 


Therefore, as a function of the order of the reaction, ' the functions 
in the intervals between the discontinuities are sements of straight lines, 
quadratic and cubic parabola, etc. The function 'X.(x) fluctuates around a unit 
value with an amplitude and frequency depending on the value of the parameter 
S (or ). 

It is important to notice the different asymptotic character of solution 
at m - 0 and m - and 2/3 when S — > 0. In the first case, the amplitude of 
the oscillations increases and approaches a finite limit; simultaneously the 


189 ~ 


frequency of the oscillations tends to infinity© In the second case, at 
sufficient small values of S, the amplitude decreases to zero when S— » 0| 
and the frequency of the oscillations, as previously, tends to infinity. 

The mechanism of the formation of the oscillations in the two-dimensional 
flow behind the shock wave is illustrated in Fig. 3» Owing to the state of 
nonequilibrium, the disturbances of the pressure come to a point with the 
coordinate x along the characteristic x ->t y » const. At the same time, the 
disturbances formed at the edge of the wedge, propagating along the characteristics 
of both families, repeatedly reflect from the surface of the wedge and the shock 
Wave . 



Diagramatic sketch of the propagation of the disturbances in the’ 
region between the shock wave and the s\irface of the wedge. 

As the equilibrium zone has a finite width, these disturbances accumulate. 
If the coefficient of reflection from the shock wave is negative, the 
compression wave alternates, on reflection, by a rarefaction wave, which leads 
to an interchange of the local compression and rarefaction on moving downstream 
in the near equilibriiun region of the flow. 

The behaviour of function yUP is shown in Figs, k and 5 for certain 
values of the adiabatic index y . The curves in Fig, 4 for m •• 1, 2, 3 are 
obtained numerically. It is obvious that, m 2, the pressure along the wall 

of the wedge increases monotonically from the frozen to the equilibrium values. 
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0 JO k.O , 1.0 8.0 , 



Zi&sJL. 

Pressure distribution along the side of a thin wedge at hypersonic 
flow aroundo The case of smooth solution for in ■ Ij 2,- 3, 



0 10 1.0 3.0 X. 


Pressure distribution along the aide of a thin wedge at hypersonic 
flow aroundo The case of \mevon solutions for m « Oj ^ ^ ^ . 

It is easy to see that in the general case, when the wedge is not thin, 
the solution of (20) can have an oscillatory character, and hero the reflection 
coefficient A and the parameter K will be the determining parameters. The 
dependence of the reflection coeffici=nt \ on .y is very strong. For 
example, in the case of a sufficiently big Mach number M^ f the change of the 
reflection coefficient from a positive to negative value occurs on decreasing 
the wedge angle 0 in a very narrow interval of change of 9 .On further 
decrease of the wedge angle within a sufficiently big interval of the change of 
9 at a fixed , the change of the reflection coefficient is small* 
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In connection with thisj the formation of oscillations in the flow must be 
sensitive for the properties of the fuel mixture ( T i Jn) and the outer 
conditions (Ma^i t O). The latter conclusion is confirmed in the experiments 
on ballistic installations 
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